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SECTION  1.0 


INTRODUCTION 


1.1  Background 

Wind  shear,  particularly  at  the  lower  altitudes  in  the  terminal  area, 
has  been  identified  as  being  hazardous  to  aircraft  operations.  Accurate  and 
reliable  wind  profiles  are  required  for  use  in  fast  time  and  manned  flight 
simulation  studies  aimed  at  fully  defining  and  understanding  the  wind  shear 
hazard.  This  report  describes  wind  speed  profiles  developed  for  the  above 
simulation  to  improve  the  safety  and  reliability  of  operations  in  the 
terminal  area.  A comprehensive  set  of  wind  profiles  and  associated  wind 
shear  characteristics  which  encompass  many  of  the  wind  shear  environments 
potentially  encounterable  by  aircraft  in  the  terminal  area  have  been  modelled, 
subject  to  the  data  available.  For  the  purpose  of  this  effort,  wind  shear  is 
defined  as.  significant  changes  in  wind  speed  and/or  direction  up  to  500  m 
above  the  ground  which  may  adversely  affect  the  approach,  landing,  or  takeoff 
of  an  aircraft.  The  wind  shear  is  mathematically  modelled  and  the  mathe- 
matical scenarios  (environment)  are  presented  in  a format  for  direct  applica- 
tion to  wind  shear  hazard/flight  simulation  studies. 

1.2  General  Description  of  Wind  Shear  Models 

A survey  of  existing  wind  shear  data  and  mathematical  models  which  are 
comprehensive,  of  sufficient  spatial  extent,  and  detailed  enough  for  the 
development  of  mathematical  models  of  wind  shear  is  reported  in  Reference 
[1-1].  This  reference  reviews  the  state  of  the  art  and  describes  the  scale 
and  duration  of  frontal  wind  shear,  neutral  and  stable  boundary  layer  wind 
shear,  and  thunderstorm  wind  shear.  The  physics  of  these  atmospheric  phenom- 
ena with  which  significant  wind  shear  is  associated  is  also  outlined  in  this 
reference  and  is  not  described  in  the  present  report.  Selected  data  sets 
reported  in  the  literature  and  described  in  Reference  [1-1]  have  been  chosen 
for  the  three  wind  shear  conditions  considered,  i.e.,  1)  thunderstorm  wind 
shear,  2)  neutral  and  stable  boundary  layer  wind  shear  and  3)  frontal  wind 
shear.  In  all  three  cases,  the  most  comprehensive  set  of  wind  shear  data 
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reported  in  the  literature  have  been  discretized  to  a grid  system  and  a 
computer  program  lookup  with  interpolation  capabilities  developed.  The 
spatial  variation  of  wind  speed,  both  horizontally  and  vertically,  is  obtained 
from  the  programs  when  called  by  a main  simulation  program.  The  location  of 
data  measurement,  how  the  data  were  obtained,  the  time  of  observation  and 
the  extent  of  the  data  up  to  500  m are  outlined  for  each  of  the  wind  shear 
conditions  considered.  From  these  data  sets,  a description  of  steady  state 
wind  which,  in  the  thunderstorm  and  frontal  cases,  is  effectively  a quasi- 
steady state  wind  having  a 10  sec  averaging  period  are  described.  For  the 
stable  and  neutral  boundary  layer  case,  the  wind  speeds  are  mean  values 
which  are  defined  as  a 10  min  average  or  greater.  For  the  neutral  and 
stable  boundary  layer,  spatial  dependence  of  the  wind  only  in  the  vertical 
direction,  z,  is  considered  and  the  wind  field  in  these  cases  is  assumed 
homogeneous  in  horizontal  extent.  For  the  thunderstorm  and  frontal  wind 
fields,  both  vertical,  z,  and  horizontal,  x,  spatial  dependence  is  considered. 
The  data  available,  however,  represents  wind  fields  in  a vertical  plane 
sliced  through  the  storm  and  spatial  variation  in  the  direction  lateral  to 
the  motion  of  the  storm  is  not  considered.  All  three  velocity  components, 

that  is,  longitudinal  velocity,  W , lateral  velocity,  W , end  vertical 

X y 

velocity,  W^,  are  defined  at  each  point  in  the  vertical  spatial  plane, 
however. 

Since  the  high  frequency  component  of  the  variation  in  the  wind  has 
been  averaged  out  of  the  data  sets  used  to  define  the  mathematical  form  of 
the  wind  shear  models,  turbulence  simulation  techniques  are  provided  from 
which  a turbulent  fluctuation  can  be  generated  and  superimposed  upon  the 
quasi-steady  state  wind  fields.  For  the  thunderstorm  and  frontal  models, 
the  frequency  content  of  the  turbulent  simulation  must  contain  all  fluctua- 
tions greater  than  0.1  Hz.  For  the  neutral  and  stable  boundary  layer  the 
turbulence  model  must  contain  fluctuations  of  all  frequencies. 

Statistical  models  which  will  allow  estimates  of  the  probability  of  a 
given  shear  magnitude  and  the  frequency  of  occurrence  are  developed  to  the 
extent  current  data  permits.  A statistical  description  of  the  wind  fields 
are  required  to  establish  meaningful  magnitudes  of  wind  shear  to  assure  that 
the  values  used  in  flight  hazard  studies  are  realistic  and  might  actually  be 
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encountered  in  the  real  world.  Statistical  data,  however,  is  very  limited 
in  this  regard  and  it  was  necessary  in  most  cases  to  establish  only  crude 
estimates  of  the  probability  of  wind  shear  and  the  risk  of  exceeding  a 
prescribed  value. 

1.3  Organization 

The  mathematical  models  of  wind  speed  developed  for  thunderstorms, 
stable  and  neutral  boundary  layers,  and  the  synoptic  frontal  stortrs  are 
reported  individually  in  Sections  2,  3 and  4,  respectively.  Each  section 
describes  the  nature  of  the  data  set  used  to  develop  the  table  lookup  programs 
and  the  theoretical  turbulent  simulation  models  recommended  for  use  with 
the  quasi-steady  winds.  Statistical  models  are  also  described  and  estimates 
of  probabilities  and  risk  of  exceedance  given.  In  the  appendix,  tables  of 
the  wind  fields  for  easy  engineering  application  are  given.  Also,  numerous 
illustrative  plots  of  typical  wind  speed  profiles  encountered  along  conven- 
tional flight  paths  are  provided  for  visualization  of  the  type  and  magnitude 
of  wind  shear  potentially  encounterable  in  the  atmospheric  conditions  con- 
sidered. Finally,  a careful  listing  of  the  computer  code  and  the  associated 
data  set  are  described.  The  computer  programs  are  given’ in  the  form  of 
subroutines  which  can  be  immediately  coupled  with  a fast  time  computer 
program  for  calculating  flight  through  wind  shear  hazards  or  for  direct 
application  to  flight  simulation  studies.  This  material  has  been  placed  in 
the  appendices  because  of  its  bulkiness  and  because,  although  it  directly 
pertains  to  the  engineering  application  of  this  report,  it  is  not  immediately 
relevant  to  the  text  of  each  section.  However,  the  appendices  have  been 
numbered  corresponding  to  the  section  of  the  report  to  which  they  refer. 
Appendix  1 contains  the  nomenclature  while  Appendices  2,  3 and  4 contain  the 
tables,  illustrations  and  computer  code  for  the  respectively  numbered 
sections.  Finally,  each  section  contains  recommendations  and  guidelines  for 
engineering  applications.  The  range  of  applicability,  engineering  interpre- 
tation, and  examples  of  how  to  use  the  data  are  given. 

The  report  concludes  with  a summary  section  which  summarizes  the  results 
of  the  effort,  presents  conclusions,  and  provides  recommendations  concerning 
future  efforts.  A discussion  of  data  deficiencies  and  the  associated  weak- 
nesses and  strengths  of  the  model  are  also  described  in  the  summary  section. 
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SECTION  2.0 


THUNDERSTORM  WIND  SHEAR 

A model  of  wind  shear  associated  with  thunderstorms  has  been  developed 
from  the  cata  of  Goff  [2-1].  These  data  consist  of  measurements  of  longi- 
tudinal, lateral  and  vertical  wind  speeds  at  various  levels  of  a 500  m 
tower.  Although  the  data  are  measured  along  a vertical  line  in  the  atmo- 
sphere, they  are  projected  horizontally  to  form  a vertical  plane  using  the 
concept  of  frozen  turbulence  or  Taylor's  hypothesis.  Thus,  the  wind  shear 
phenomena  associated  with  thunderstorm  models  can  be  described  as  a two- 
dimension,  spatial  wind  field. 

The  data  of  Goff  [2-1],  after  careful  review  (see  Reference  [2-2]), 
were  selected  as  being  the  best  data  available  to  construct  a quantitative 
model  of  thunderstorm  wind  shear  that  provides  both  vertical  and  horizontal 
wind  shear  values.  As  noted,  the  model  is  restricted  to  two-dimensional 
space.  The  extensive  survey  reported  in  Reference  [2-2]  indicates  that 
there  are  no  three-dimensional  data  sets  available  nor  any  theoretical 
models  associated  with  thunderstorms  which  would  allow  a three-dimensional 
simulation  to  be  carried  out.  Therefore,  it  was  necessary  to  restrict  the 
simulation  of  thunderstorm  wind  shear  to  two-dimensional  wind  fields.  The 
data,  however,  do  include  all  three  velocity  components  of  the  wind  vector 
in  the  plane  swept  out  by  the  500  m tower  as  the  storm  passes. 

2.1  Data  Source 

Thunderstorm  wind  shear  data  are  presented  in  the  form  of  longitudinal 
(W  ),  lateral  (W  ) and  vertical  (W  ) wind  speed  components  in  a vertical 

X y Z 

plane  for  20  thunderstorm  cases.  Data  for  these  20  storms  were  measured 

during  the  months  of  May  through  June  over  the  period  of  1971  through  1973 

with  the  WKY-TV/NSSL  481  m meteorological  tower,  Norman,  Oklahoma.  Time 

histories  of  the  wind  speed  are  converted  to  horizontal  spatial  distributions 

using  Taylor's  hypothesis  (i.e.,  x = Wt).  Ten  second  averaged  values  of 

wind  speed  components  are  provided  in  the  form  of  isotach  maps  for  W , W 

A y 

and  W^,  respectively.  These  data  were  converted  to  a 41  x 11  point  grid 
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system  as  illustrated  in  Figure  2-1.  The  data  for  the  20  thunderstorm  cases 
having  been  discretized  on  the  grid  system  were  stored  on  magnetic  tape  and 
a computer  lookup  routine  developed  for  interpolating  the  W , W and  W wind 
speeds  for  any  location  within  the  x-z  plane.  Tabulated  data  for  W , W 

A y 

and  on  the  41  x 11  point  grid  system  are  given  in  Appendix  2A;  a pictorial 
description  of  the  wind  speed  profiles  and  streamline  patterns  are  given  in 
Appendix  2B,  and  a computer  program  which  stores  the  data  on  disks  and 
carries  out  the  table  lookup  with  the  option  of  superimposing  turbulence  is 
described  in  Appendix  2C. 

The  20  thunderstorm  cases  for  the  purpose  of  this  report  are  assigned 
numbers  1 through  20  which  correspond  to  Goff's  identification  symbols  A 
through  T.  All  tabulated  and  illustrated  data  in  the  appendix  contain  both 
the  numerical  and  the  alphabetical  identification. 

2.2  Data  Processing 

For  each  thunderstorm,  a 10  min  record  was  taken  and  the  data  were 
averaged  over  10  sec  intervals  throughout  the  10  mins.  These  data  were  then 
fitted  to  a regular  array  of  a 60  x 10  grid  by  Goff.  Although  the  actual 
measurements  were  not  uniformly  spaced,  the  data  as  presented  for  public 
distribution  were  uniformly  spaced  on  the  60  x 10  grid.  Each  10  min  record 
consists  of  60  10  sec  averaged  data  sets  representing  a cross  section  through 
the  storm.  Streamline  patterns  as  well  as  isotach  maps  of  W , W and  W 

A y L 

were  also  drawn  from  the  data  and  are  given  in  Reference  [2-1]. 

The  data  presented  in  Reference  [2-1]  were  evenly  spaced  over  10  50  m 
intervals  in  the  z-direction.  Values  assigned  to  each  interval  were  deter- 
mined by  linear  interpolation  from  seven  levels  of  measurements  for  the 
1971,  1972  data  and  from  only  three  levels  of  measurement  for  the  1973  and 
1974  data.  For  the  x-direction,  60  intervals  were  utilized  where  a Galilean 
transformation  from  time  to  space  was  utilized.  The  horizontal  extent  of 
each  cross-section  is  therefore  a variable  and  is  equal  to  the  frontal  speed 
of  the  thunderstorm  times  the  10  min  averaging  period. 

The  coordinate  system  utilized  in  Reference  [2-1]  and  also  in  this 
report  defines  the  x-direction  as  being  measured  positively  in  the  direction 
of  frontal  motion  and  the  y-direction  as  being  measured  perpendicular  to  the 
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frontal  axis  utilizing  the  right-hand  rule,  i.e.,  positive  is  measured 
into  the  plane  of  the  paper.  The  components  of  the  wind  speed  are  posi- 
tive in  the  direction  of  frontal  motion,  i.e.,  positive  from  left  to  right  on 
the  contour  maps  and  the  vertical  wind  speed,  W^,  is  positive  in  the  upward 
di rection. 

The  data  presented  by  Goff  [2-1]  are  the  longitudinal  wind  speed,  W^, 
relative  to  the  motion  of  the  storm.  In  this  report,  however,  the  data  are 
generally  presented  relative  to  the  fixed  frame  of  reference  attached  to  the 
ground. 

In  view  of  the  10  sec  averaging  period  utilized  to  reduce  the  data,  all 
fluctuations  in  wind  speed  of  frequency  higher  than  0.1  Hz  have  been  filtered 
out  of  the  data.  Thus,  the  data  contain  variations  in  wind  speed  having  fre- 
quencies of  0.1  Hz  or  less.  These  low  frequencies  are  expected  to  be  those 
to  which  aircraft  motion  is  most  sensitive.  However,  to  assure  correct 
simulation  where  high  frequency  wind  components  are  significant,  it  is 
recommended  that  a turbulent  fluctuation  be  superimposed  upon  the  quasi - 
steady  wind  fields.  Section  2.4  discusses  and  recommends  a turbulence  model 
for  use  with  the  thunderstorm  gust  front  data. 

The  following  section.  Section  2.3,  describes  the  quasi-steady  wind 
speed  profile  grid  system  and  presents  working  data  for  wind  shear  hazard/ 
flight  simulation  studies.  The  wind  speed  is  referred  to  as  quasi-steady  in 
view  of  the  fact  that  it  is  averaged  over  10  sec  time  intervals  and  thus 
contains  departures  from  the  mean  wind  speed  which  is  generally  averaged  over 
a 10  min  period  or  greater. 

2.3  Quasi-Steady  Wind  Speed  Profile  Grid  System 

The  wind  fields  of  Goff  [2-1]  were  fit  to  a 41  x 11  point  grid  system. 
The  grid  system  as  illustrated  in  Figure  2-1  was  superimposed  on  the  contours 
and  wind  speeds  at  each  grid  point  tabulated  and  punched  onto  computer  cards. 
The  data  were  later  stored  on  magnetic  tape. 

The  grid  system  is  numbered  from  the  left-hand  bottom  corner.  The 
numbers  increase  from  left  to  right  in  the  positive  x-direction  of  the 
original  data  and  from  bottom  to  top  in  the  positive  z-direction  of  the 
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original  data.  The  wind  speeds  are  given  in  units  of  meters  per  sec, 
m s"\  with  being  positive  in  the  direction  of  frontal  motion,  being 
positive  upward  and  being  positive  into  the  plane  of  the  paper. 

The  wind  speed  W is  stored  on  tape  as  the  wind  speed  relative  to  the 
storm  motion.  The  stored  values  of  wind  speed  are  therefore  the  values  an 
observer  would  measure  moving  along  with  the  storm.  To  convert  to  the  earth 
frame  of  reference  the  mean  motion  of  the  storm  must  be  added  to  the  longi- 
tudinal wind  speed.  The  next  section  describes  tables  and  graphical  illustra- 
tions of  the  wind  speed  profiles.  In  these  cases  the  mean  motion  has  been 
added . 

2.3.1  Tables  of  Wind  Speed 

Tables  of  thunderstorm  wind  speed  are  given  in  Appendix  2A.  The  upper 
portion  of  the  table  covers  grid  stations  1 through  21  and  the  lower  table 
covers  grid  stations  21  through  41.  (Note  column  21  is  repeated  for  symmetry 
and  clarity  of  presentation.) 

The  thunderstorm  case  numbers  designation  for  this  report  are  listed  at 
the  top  of  each  table.  The  letters  in  parentheses  and  the  following  series 
number  corresponds  to  the  thunderstorm  designation  given  by  Goff  [2-1]. 

Also  listed  at  the  top  of  the  table  is  the  frontal  speed,  W^,  and  the 
horizontal  length  scales  for  the  given  wind  record.  The  horizontal  extent  of 
each  data  set  varies  because  of  the  data  reduction  procedure.  Hence,  the 
length  of  field,  L,  in  kilometers  and  the  horizontal  grid  spacing.  Ax,  are 
specified  on  each  table.  The  vertical  extent  of  each  field  is  taken  as  500  m 
with  50  m vertical  grid  spacing. 

2.3.2  II  lustrations 

Illustrations  of  the  longitudinal,  lateral,  and  vertical  wind  speeds 
encountered  along  a 3°  glide  slope  through  each  thunderstorm  are  provided  in 
Appendix  2B.  The  glide  slope  is  adjusted  to  terminate  at  the  lower  left-hand 
corner  of  each  streamline  plot  as  illustrated.  The  ordinate  in  the  wind 
speed  profiles  is  height,  z,  nondimensional ized  with  the  length,  H = L tan  3°, 
where  L is  the  length  of  the  wind  field.  Each  profile  is  the  wind  seen  by  an 


airplane  traveling  along  the  flight  path  line  drawn  across  the  streamlines  as 
shown  in  the  upper  figures.  The  flight  paths  are  drawn  to  terminate  in  the 
lower  left-hand  corner  always.  Note  that  the  streamlines  plotted  are 
relative  to  the  speed  of  the  front  which  for  reference  purposes  is  indicated 
with  the  vertical  dashed  line  on  the  horizontal  wind  speed  profile.  The  wind 
speed  profiles  are  relative  to  the  fixed  earth  frame  of  reference. 

2.3.3  Computer  Program  for  Computing  Thunderstorm  Wind  Speed  Profiles 

A computer  program  has  been  developed  which  for  given  input  x,  z 
computes  the  longitudinal,  lateral,  and  vertical  wind  speeds  at  that  position. 
The  six  velocity  gradients  3W  /8x,  9W  /9z,  9W  /9x,  9W  /9z,  9W  /9x,  and  9W  /9z 

X X fc  j"  jf 

are  also  output  by  the  program.  A detailed  description  and  user's  instruc- 
tions for  the  computer  program  are  provided  in  Appendix  2C.  The  computer 
program  also  has  the  option  of  calling  for  turbulence  which  is  super- 
positioned  upon  the  thunderstorm  wind  field.  A non-Gaussian  turbulence  model 
based  on  the  technique  of  Reeves,  et  al . [2-3]  is  employed  for  turbulence 
simulation. 

It  is  intended  that  this  computer  program  can  be  used  as  a subroutine 
for  direct  application  to  fast  time  flight  studies  or  for  fast  time  computer 
simulation  of  aircraft  flight  through  thunderstorms.  As  an  example  of  the 
application  of  the  computer  program,  the  reader  is  referred  to  Reference  [2-4] 
in  which  this  computer  program  has  been  used  to  study  aircraft  dynamics  in 
thunderstorms  and  Reference  [2-5]  where  the  data  were  used  for  manned  flight 
simulator  programming. 

2.4  Turbulence  in  Thunderstorms 

Measurements  of  the  power  spectral  density  function  for  turbulence  in 
thunderstorms  has  been  reported  as  early  as  1962  by  Steiner  and  Rhyne  [2-6]. 
Their  data  was  measured  only  over  the  approximate  reduced  spatial  frequency 
range  of  0.004  to  0.4  rad  m"\  The  theoretical  von  Karman  spectrum  follows 
the  data  in  this  frequency  range  very  well  as  demonstrated  by  Houbolt, 
et  al.  [2-7],  in  Figure  2-2a.  The  Dryden  spectrum,  on  the  other  hand,  does 
not  compare  as  well  with  the  data.  However,  with  the  properly  chosen  length 
scale  a reasonable  fit  is  obtained  (see  Figure  2-2b).  All  reported  data. 
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however,  were  measured  in  the  altitude  range  of  12  to  8 km  and  moreover  they 
do  not  extend  to  low  enough  frequencies  to  illustrate  at  which  point  the  knee 
of  the  turbulence  spectrum  curve  occurs.  The  knee  of  the  power  spectral 
density  function  is  important  in  many  design  applications  (see  Reference 
[2-8])  and  particularly  to  the  empirical  value  of  length  scale,  L,  used  in 
the  analytical  models  of  the  power  spectral  density  function. 

In  Reference  [2-7]  a comparison  of  the  power  spectral  density  function 
of  severe  storms  with  that  of  cumulus  clouds  and  clear  air  turbulence  is 
given  (see  Figure  2-3).  One  can  see  from  Figure  2-3  that  the  turbulent 
spectra  for  severe  storms  behaves  very  similar  to  that  of  cumulus  clouds  and 
clear  air  turbulence  with  the  only  major  difference  being  higher  values  of 
(t)(fi)  which  indicates  higher  turbulence  intensity.  Again,  these  data  are 
measured  at  very  high  altitudes  and  probably  do  not  include  effects  due  to 
the  presence  of  the  ground. 

Reference  [2-7]  computes  two  values  of  turbulence  intensity.  One  is 
computed  directly  from  the  experimental  data  by  integrating  the  power 
spectral  density  function  over  only  the  limits  of  the  frequency  range  in 
which  the  data  have  been  measured,  i.e.. 


o,  = 


(})(n)dn 


(2-1) 


where  u and  Z designate  upper  and  lower,  respectively.  Spatial  frequency 
designated  by  is  related  to  cyclic  frequency  used  in  Section  3-4  by 
n = Uil/2v.  An  alternate  method  of  measuring  the  turbulence  intensity  is  to 
compute  its  value  from  the  autocorrelation  function  evaluated  at  zero. 


(X> 

4>(5l)dfi 


o 


(2-2) 


It  is  found  that  based  on  the  Dryden  model  the  two  measures  of  turbulence 
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Power  Spectral  Density,  <{>  (fi)  , (m/sec)  / Radians/m 


0,001  0.01  0.1  1.0 

Reduced  Frequency,  fi.  Radians/m 


Figure  2-3  Typical  Power  Spectra  of  Vertical  Component  of  Turbulence 

Measured  in  Clear  Air,  Cumulus  Cloud,  and  Thunderstorm  [2-7] 
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intensity  are  related  by 


(2-3) 


where  n and  Q.  represent  the  upper  and  lower  limits,  respectively,  of  the 
u ~ 

range  of  the  turbulent  spectrum  measurements.  Also  appearing  in  this  equa- 
tion is  the  length  scale,  L,  which  is  an  extremely  important  parameter  in 
fitting  the  analytical  expressions  for  the  turbulence  spectra  to  the  data 
(see  Reference  [2-9]). 

Although  the  von  Karman  spectra  tends  to  fit  the  data  better  as 
illustrated  in  Figure  2-2,  the  Dryden  spectrum  is  conventionally  used  because 
of  its  rational  form. 


L a 

w w 

1 

X 

X 

7T  j 

% n + 3(L„ 
y y j_ 

['  + (L„  !!)27 

y 

L [1  + 3(L  fl)^] 
w w w 

z z z 

[>  • (L„ 


(2-4) 


Evaluation  of  the  Dryden  spectrum  in  a turbulence  simulation  routine 
requires  knowledge  of  the  length  scale,  and  of  the  turbulence  intensity, 

*.  Houbolt,  et  al.  [2-7]  recommend  values  of  L = 1036  m and  gives  oi  values 
for  the  vertical  fluctuations  in  the  range  of  4.08  m s"^  to  2.27  m s ^ and 
for  the  lateral  fluctuations  of  4.69  m s"^  to  2.69  m s \ Longitudinal 
values  of  a were  not  measured.  Adjusting  oi  to  o from  Equation  2-3  gives 
values  of  approximately  10.20  > > 4.75  m s"^  and  9.82  > Oy/^  > 5.63  m s . 
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These  values  cannot  be  used  directly,  however,  because  the  ground  is 
expected  to  have  a strong  effect  on  the  turbulence  length  scale  and  intensity. 

No  actual  data  for  L and  a nor  for  how  they  vary  with  height  in  a 
thunderstorm  below  500  m appears  to  be  available  in  the  literature.  Barr, 
et  al . [2-10]  postulate  a decrease  in  length  scale  and  intensity  at  low 
altitudes.  Preliminary  results  from  the  analytical  model  of  Lewellen, 
et  al.  [2-11]  nredict  increasing  values  of  and  near  the  ground  (see 
Figure  2-4).  In  Section  3.4  it  is  recommended  that  near  the  ground  the  ratio 
of  and  be  determined  from  Figure  3-5.  The  relationships  are 

0.,  /o.,  = [0.177  + 2.74  X 10"^  z]"°-^ 
w w 
X z 

and 

Ow  /o^  = [0.583  + 1.39  X 10"^  z]"°'®  (2-5) 

and  for  leek  of  a better  model  this  relationship  is  recommended  herein. 

Also  in  Section  3.4,  an  equation  (3-11)  for  evaluating  oyy^  proposed  by 
Barr,  et  al . [2-10]  is  given.  Evaluation  of  from  this  equation,  however, 
is  very  dubious  for  thunderstorms  and  is  not  recommended.  For  example,  a 
very  large  and  uncommonly  observed  value  of  u*  is  2 m s ^ or  greater  (see 

Section  3.4)  which  from  Equation  3-11  gives  a value  of  = 2.6  m s . This 

value  which  would  occur  only  infrequently  in  a normal  boundary  layer  is  a 
factor  of  1.8  lower  than  the  lowest  value  of  measured  by  Steiner  and 

Rhyne  [2-6]  bearing  in  mind,  of  course,  that  their  data  were  measured  at 

altitudes  between  12  and  8 km. 

It  is  assumed,  however,  that  Steiner  and  Rhyne's  data  will  extend  to 

lower  altitudes  and  the  procedure  recommended  for  evaluating  is  to  select 

-1  . ^ 

a value  in  the  range  of  10.20  > < 4.75  m s . To  predict  ay/^  and  owy. 

Equation  2-5  should  be  used. 

The  length  scale  of  the  turbulence  may  be  selected  from  Figure  2-5  taken 
from  Reference  [2-10]  where 
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Altitude 


RMS  Velocity  Fluctuations,  m 

Figure  2-4  Profile  of  the  RMS  Values  of  the  Vertical  and  Head- 
wind Velocity  Fluctuations  as  Predicted  by  Lewellen, 
et  al.  [2-11] 
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I =1  = 2 ^ 300  m 

w “ I3OO  m,  z 300  m 

X y ^ 


I ^ /z[0J77  + 0.00274  z]"^’^;  z < 300  m 

|300  m;  z > 300  m (2-5) 

Table  2-1  provides  values  of  cyclic  frequencies  for  different  veloc- 
ities of  and  for  spatial  frequency  of  0.4  and  0.004  radians/m, 
respectively.  It  is  observed  from  this  table  that  the  tabulated  values 
exceed  0.1  Hz  only  for  the  upper  limit  of  and  at  velocities  greater  than 
10  m s"\  The  importance  of  this  observation  to  the  simulation  model  of  wind 
shear  being  developed  herein  is  that  the  wind  data  utilized  in  the  math 
model  has  been  integrated  over  10  sec  periods  and,  therefore,  contains 
oscillations  in  the  low  frequency  range  of  less  than  0.1  Hz.  Therefore,  it 
is  proposed  that  a turbulence  simulation  for  the  current  thunderstorm  model 
have  all  frequencies  below  0.1  Hz  filtered  out  using  a high  pass  filter  (see 
Reference  [2-12]).  In  many  ways  this  result  is  advantageous  because  the 
turbulence  spectra  are  better  known  at  high  frequencies  and  can  be  simulated 
very  well;  whereas,  low  frequency  turbulence  spectra  are  not  yet  well 
defined.  The  high  frequency  range  corresponds  to  the  well  established 
inertial  subrange  [2-13]  for  which  (t)(f2)  obeys  the  -5/3  power  law.  A high 
pass  digital  filter  program  was  used  with  the  computer  program  given  in 
Appendix  2C  but  is  not  included  in  the  write-up  because  it  is  a standard 
program  given  in  Reference  [2-14]  and  is  too  long  to  be  included  in  this 
report. 

2.5  Statistical  Model  of  Turbulence 

In  an  attempt  to  make  some  estimate  of  the  most  extreme  wind  shear 
which  could  be  encountered  during  the  expected  life  cycle  of  a thunderstorm, 
the  20  thunderstorm  cases  were  taken  as  a sample  for  statistical  purposes. 
From  this  sample,  wind  speeds  along  the  20  flight  paths  illustrated  in 
Appendix  2B  were  statistically  analyzed.  Table  2-1  shows  the  ensemble 
average  of  the  wind  speeds  along  the  flight  paths  plus  the  one  point 
standard  deviation  of  each  wind  speed  and  the  correlations  betv/een  the  wind 
speed  components.  These  values  are  defined  as  shown  below  for  the 


10 


TABLE  2-1 


ENSEMBLE  AVERAGES,  STANDARD  DEVIATIONS  AND  CORRELATIONS 

FOR  W , W AND  W 
X y z 


z 

<W  > 

<W  > 

<W  > 

0 

X 

y 

Z 

0 

9.0700 

-0.6850 

0.0000 

50 

9. 189b 

-1.0350 

-0.0519 

100 

9.7518 

-1.2028 

-0.1174 

150 

11.4401 

-1 . 3898 

-0.2904 

200 

12.3359 

-1.2510 

-0.2539 

250 

12.5787 

-0.1779 

-0.0397 

300 

11.8227 

0.2164 

-0.0670 

350 

8.8224 

1.1370 

0.6620 

400 

4.3032 

2.1131 

1.2981 

450 

1 . •*b70 

2.7131 

0.8968 

500 

-0.1681 

2:9205 

0.4902 

^0 

X 

y 

z 

0 

4.0341 

3.7439 

0.0000 

50 

3.4693 

3.4720 

0.1404 

100 

3.7829 

3.6212 

0.2641 

150 

3.9754 

3.7328 

0.4325 

200 

4.0718 

3.7115 

0.6148 

250 

3.8818 

4.9743 

1.0989 

300 

3 . 8 2 4 4 

5.1184 

0.8581 

350 

5.4606 

4.7680 

0.9819 

400 

5.7848 

5.8223 

1 . 2641 

450 

5.5734 

1 .il32 

0 . 6 5 b 3 

500 

5.3830 

8.2354 

0.4774 

Zo 

p 

W^w 

p 

W W 

w 

X y 

X z 

y z 

0 

3.4760 

0.0000 

0 .0000 

50 

1 .9b69 

0.0184 

-0.05b7 

100 

1.5427 

-0. 1086 

0.1101 

150 

0.3304 

-0.7416 

-0.0176 

200 

0.1972 

-0.8220 

1.2121 

250 

0.9971 

-0.8448 

3 . 8 6 9 6 

300 

0 . 09*^5 

-1.4322 

2.1130 

350 

2 . b395 

-2.bl72 

1 .2844 

400 

-4.8549 

-2.5340 

1 .9091 

450 

-4.6407 

-0.5253 

-0.5719 

500 

-2.0084 

-0.1117 

-1 .8076 
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longitudinal  wind  speed. 


<W  > = rr  y W 

X N X. 


N 


N 

I(Wx 

i=l  ^i 


(2-7) 


<W  >)' 

X 


N - 1 


1/2 


(2-8) 


and 


Pw  w = ^ 

X y 1=1  1 


<W  >)(W^ 


<W  >) 

y 


(2-9) 


The  relationship  for  the  lateral  and  vertical  wind  speed  and  other  two 
correlations  are  similarly  defined. 

The  results  of  Table  2-1  allow  a trivariate  probability  density  function 
for  the  wind  components  at  each  point  along  the  flight  path  to  be  con- 
structed. A trivariate  distribution  as  described  in  Reference  [2-15]  is 
given  by 


P(W^,  Wy,  W^)  = 


(2-10) 


where 


and 


(2-11) 


K • 1 ■ 1.  2.  3 


The  vertical  bars  designate  the  determinant  of  the  matrix  s^j  where 
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- K 


(2-12) 


From  this  function  the  probability  of  the  wind  vectors  at  each  point  can  be 
estimated,  however,  the  mathematics  involved  is  too  complicated  to  present  in 
detail  in  this  report  and  the  reader  is  referred  to  Reference  [2-15], 

There  is  considerable  encouragement  that  a probability  model  of 
thunderstorms  could  be  constructed  since  the  components  of  [s^j]  and  the 
ensemble  averages  for  the  20  storms  are  reasonably  consistent  functions  of 
the  spatial  coordinate  x,z.  This  is  demonstrated  in  Figure  2-6  v/hich  shows 
the  variation  of  <W^>,  and  Pw^Wy  aloP9  horizontal  lines  through  the  storm 
at  equal  increments  of  At  = 1 min  or  Ax/W^.  The  lines  are  at  the  400,  300 
and  200  m level.  The  data  show  well  defined  variation  with  time  and  con- 
sistent trends  with  altitude.  Equally  consistent  results  are  obtained  with 
the  remaining  statistical  parameters. 

These  well  behaved  statistical  parameters  suggest  that  a study  to 
develop  a complete  statistical  model  of  a thunderstorm  be  carried  out.  The 
results  of  this  study  would  allow  simulation  of  thunderstorms,  similar  to 
turbulence  simulations,  to  be  performed.  Thus,  a simulator  could  be  pro- 
grammed to  simulate  approaches  and  landings  through  a random  selection  of 
thunderstorms.  Moreover,  the  statistical  model  would  allow  an  estimate  of 
the  most  extreme  wind  shear  associated  with  a thunderstorm  and  its  frequency 
of  occurrence. 

Until  such  a model  is  available,  the  best  estimate  of  the  probability 
of  an  extreme  value  in  the  thunderstorm  can  be  made  from  the  tabulated 
values  of  standard  deviation  (Table  2-1)  and  the  assumption  of  a Gaussian 
distribution.  For  example,  one  can  estimate  that  at  the  point  z = 250  m on 
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the  3°  glide  slope  W will  on  the  average  be  12.6  ms,  there  is  however  a 

^ -1 

32%  chance  that  it  will  be  +3.88  ms  of  that  value;  a 5%  chance  tnat  it 
will  be  ±7.76  m s"^  of  that  value;  etc.  This  approach  assumes  that  the  wind 
fluctuations  at  neighboring  points  are  statistically  independent  which  is 
obviously  not  the  case  in  view  of  the  significant  value  of  the  correlation 
coefficients  also  shown  in  Table  2-1. 

Considerably  more  work  is  required  to  produce  a satisfactory  statistical 
model  of  a thunderstorm  wind  field,  however,  preliminary  analyses  conducted 
on  the  sample  of  20  thunderstorms  suggest  a meaningful  model  is  feasible  and 
within  reach. 

2.6  Application  of  the  Thunderstorm  Wind  Shear  Model 

The  user  of  the  thunderstorm  wind  shear  model  simply  selects  which 
thunderstorm  case  is  of  interest  to  his  simulation  problem  and  generates 
wind  speed  profiles  for  given  values  of  x and  z throughout  the  wind  field. 

The  selection  of  the  thunderstorm  case  must  consider  two  factors:  one  is  the 
length  of  horizontal  fetch  over  which  the  simulation  is  to  be  carried  out 
and  the  other  is  the  severity  of  the  storm  which  it  is  desired  to  simulate. 

All  profiles  have  the  same  vertical  extent  which  is  500  m.  However, 
they  have  different  lengths  depending  upon  the  mean  motion  of  the  thunder- 
storm. If  a simulation  is  to  be  carried  out  over  a significantly  long 
distance  in  the  horizontal  direction,  one  must  select  the  storm  which  will 
span  the  distance  of  the  proposed  simulation.  Table  2-2  lists  the  length  of 
each  thunderstorm  based  on  its  case  number.  This  table  can  be  used  to 
select  which  thunderstorm  case  to  use  in  carrying  out  the  simulation  where 
the  length  of  horizontal  extent  is  critical. 

The  severity  of  the  storms  also  differ  for  a number  of  reasons.  One 
is  that  the  maturity  of  the  storms  as  they  passed  over  the  tower  were  in 
different  stages  of  development,  another  is  that  the  center  of  the  thunder- 
storm may  not  have  passed  over  the  tower  and  only  the  fringes  of  the  storm 
were  recorded.  In  selecting  a storm  to  have  a severity  appropriate  for  the 
simulation  to  be  carried  out,  the  user  may  either  inspect  the  graphical 
illustrations  shown  in  Appendix  28  and  thus  pick  what  would  appear  to  be  the 
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TABLE  2-2 


HORIZONTAL  LENGTH  OF  EACH  STORM  RECORD 


Case 

Length  (km) 

Case 

Length  (km) 

1 

4.00 

11 

8.16 

2 

3.28 

12 

3.64 

3 

5.71 

13 

6.35 

4 

7.69 

14 

5.13 

5 

11.43 

15 

7.41 

6 

7.27 

16 

4.44 

7 

7.84 

17 

11.43 

8 

5.56 

18 

16.00 

9 

7.41 

19 

7.14 

10 

8.51 

20 

13.33 

worst  storm  relative  to  the  simulation  he  wishes  to  perform.  Note,  however, 
that  the  wind  profiles  shown  are  based  on  a specified  flight  path  which 
terminates  at  the  bottom  left-hand  corner  of  each  data  set  and  which  has  a 
constant  3°  glide  slope.  Therefore,  the  flight  path  for  which  the  winds  are 
illustrated  may  not  pass  through  the  worst  part  of  the  storm.  In  this  case, 
the  user  may  examine  the  tables  provided  in  Appendix  2A  and  by  examining  the 
overall  wind  fields  select  from  the  tables  the  thunderstorm  case  which  gives 
either  the  largest  longitudinal  wind,  largest  vertical  wind,  or  the  largest 
lateral  wind,  which  ever  may  be  of  interest. 

Having  chosen  the  thunderstorm  of  interest,  one  would  normally  wish  to 
predict  what  the  probability  of  encountering  such  a thunderstorm  would  be 
and  with  what  frequency  would  such  storms  occur.  Statistical  models  which 
would  allow  this  type  of  risk  of  exceedance  estimate  to  be  made  are  not 
available  and  require  further  research  for  their  development.  The  best 
procedure  to  achieve  some  estimate  of  the  probability  of  encountering  such  a 
storm  is  to  utilize  the  information  provided  in  Section  2.5.  Table  2-1  in 
this  section  gives  the  ensemble  averaged  wind  speeds  for  all  20  thunderstorm 
cases  and  the  standard  deviations  of  these  wind  speeds  about  this  mean. 
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Again,  these  results  are  compiled  for  only  those  flight  paths  as  specified 
earlier.  The  user  can,  by  comparing  his  selected  storm  with  the  ensemble 
averaged  value,  determine  the  number  of  standard  deviations  his  storm  departs 
from  the  mean  and  in  this  way  estimate  the  probability  of  the  wind  field 
magnitude.  The  addition  or  subtraction  of  standard  deviations  about  the 
mean  at  each  point  along  the  wind  field  provides  only  a crude  estimate  of 
the  statistics  of  the  thunderstorm  wind  fields,  however.  This  point  is  well 
illustrated  by  inspection  of  the  coherence  between  wind  speed  components 
provided  in  Table  2-1.  There  is  a very  strong  coherence  between  wind  speed 
components  and  therefot^e  each  point  in  the  wind  field  does  not  behave  inde- 
pendently. Therefore,  it  is  incorrect  to  simply  add  standard  deviation  at 
each  point  but  in  lieu  of  a better  approach  this  method  can  be  followed. 

The  foregoing  arguments  clearly  indicate  that  there  is  a need  to  carry 
out  a more  detailed  analysis  of  thunderstorms  such  that  a statistical  model 
which  would  provide  the  extreme  magnitudes  of  wind  speed  and  of  the  wind 
shear  expected  to  occur  in  a thunderstorm  will  be  available.  The  frequency 
of  occurrence  of  the  extreme  is  also  greatly  needed. 

The  turbulence  model  developed  in  conjunction  with  the  thunderstorm 
wind  shear  model  should  definitely  be  used  in  carrying  out  any  simulation 
process.  The  reason  for  this  is  that  the  data  utilized  to  develop  the  wind 
shear  model  predict  at  most  downdrafts  of  3 m s Values  of  downdrafts  as 
high  as  15.5  m s"\  however,  are  reported  in  [2-16  and  2-17],  These  values, 
however,  are  undoubtedly  averaged  over  much  shorter  periods  of  time  than  10 
sec  for  which  the  data  presented  herein  are  averaged.  Neither  reference 
gives  any  information  on  the  averaging  time  utilized  in  arriving  at  the 
quoted  value  of  15.5  m s~\  Two  models  of  thunderstorm  wind  fields  that 
have  been  developed  by  Keenan  [2-17]  are  tabulated  in  Tables  2-3  and  2-4. 

These  wind  fields  were  reconstructed  from  the  flight  data  recorder  of  aircraft 
involved  in  accidents  resulting  from  flight  through  severe  thunderstorms. 
Inspection  of  these  tabulated  results  illustrate  that  much  more  extreme 
downdrafts  or  downbursts  as  defined  in  Reference  [2-16]  occur  in  these  data 
sets  than  in  those  tabulated  in  Appendix  2A. 

There  is  at  this  time  conflicting  data  and  opinions  as  to  the  maximum 
magnitude  of  the  downdraft  that  can  occur  in  a thunderstorm.  Although  the 
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2-3  Thunderstorm  Wind  Field,  BIO,  Similar  to  Philadelphia/Allegheny 
Profile  [2-1^;  HX  - Horizontal  Station  (ft.).  Velocities  Are 
Given  in  Knots  and  Height  in  Feet. 
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10  sec  average  data  of  Goff  [2-1]  will  have  lower  values  than  the  peak 
downdraft  wind  speeds  reported  by  Fujita  [2-16],  the  discrepancy  in  the 
values  cannot  be  completely  contributed  to  averaging  time.  A recent  report 
by  Alexander  [2-18]  gives  a statistical  summary  of  vertical  wind  speed  data 
recorded  at  NASA's  150  m ground  wind  tower  facility,  Kennedy  Space  Center, 
Florida.  One  year  of  continuous  around-the-clock  vertical  wind  speed  measure- 
ments were  processed  to  determine  the  intensity,  frequency,  time  of  occur- 
rence, etc.,  of  the  daily  maximum  vertical  gust.  Both  updrafts  and  downdrafts 
were  studied.  These  values  represent  0.1  sec  averages  and  the  maximum 
vertical  downdraft  recorded  is  9.3  m s'^  although  data  recorded  specifically 
during  the  hurricane  Agnes  indicated  a downdraft  in  excess  of  11.9  m s"^. 

Sinclair  [2-19]  indicates  that  downdrafts  at  1000  m for  an  Oklahoma 
thunderstorm  may  be  considerably  in  excess  of  the  15.5  m s ^ recorded  in 
Reference  [2-16].  Sinclair  has  experienced  and  measured  downdrafts  as  high 
as  28  m s"^  based  on  a 1/25  sec  averaging  period.  Finally,  the  numerical 
model  of  Williams,  et  al.  [2-20]  does  not  predict  wind  speed  downdrafts 
greater  than  10  m s”\ 

Thus,  it  is  evident  that  research  is  needed  to  bring  together  the  data 
currently  available  and  to  resolve  the  magnitude  of  the  maximum  downdraft 
which  can  occur  within  a thunderstorm.  This  v/ould  allow  the  current  simula- 
tion model  to  be  updated  by  superimposing  turbulence  fluctuations  of 
realistic  magnitude  on  quasi-steady  wind  fields.  For  the  time  being,  however, 
to  provide  an  estimate  of  wind  fields  which  would  be  consistent  with  the 
higher  values  of  vertical  wind  speed  reported  in  References  [2-16  and  2-17] 
consider  the  following. 

If  the  wind  shear  model  based  on  the  thunderstorm  data  from  Goff  [2-1] 

incorporates  into  the  turbulence  simulation  fluctuations  which  are  based  on 

4.75  s s 10.2  m s"^  recorded  at  12  to  8 km  as  described  earlier  in  Section 

2.4,  very  high  downdrafts  will  be  computed.  For  example,  taking  the  average 

value  of  ow,  = 7.5  m s"^  and  adding  to  that  the  10  sec  average  wind  speed  of 
2 1 _] 

approximately  2.5  m s , a value  of  10  m s is  obtained  for  one  standard 
deviation  and  a value  of  17.5  m s”^  for  two  standard  deviations.  The 
reported  value  of  15.5  m s“^  mentioned  earlier  is  slightly  less  than  two 
standard  deviations  about  the  10  sec  mean.  Thus,  statistically  it  is  clear 
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that  downdrafts  of  15  m s'^  or  greater  can  readily  occur  in  thunderstorms  if 
the  turbulence  intensities  at  altitudes  of  12  to  8 km  extend  to  the  ground. 
However,  turbulence  intensity  is  normally  attenuated  near  the  ground  and  it 
is  not  confirmed  that  such  high  values  exist  there.  Until  this  is  experi- 
mentally resolved,  it  is  recommended  that  the  model  of  thunderstorm  wind 
shear  based  on  the  extensive  data  from  Goff  [2-1]  should  have  turbulent 
fluctuations  superimposed  with  a standard  deviation  of  = 7.5  m s The 
downdraft  magnitudes  reported  in  accident  investigation  will  then  be  achieved 
in  the  simulation. 

It  is  anticipated  that  the  proposed  turbulence  simulation  model  based  on 
the  work  of  Reeves,  et  al . [2-3]  will  provide  a realistic  turbulence  simu- 
lation to  accompany  the  wind  shear  model  proposed  in  this  report.  However,  a 
better  simulation  of  the  large  downdraft  fluctuations  that  occur  in  a thunder 
storm  can  be  achieved  by  utilizing  a model  of  turbulence  which  includes 
coherence  between  different  levels  of  the  atmosphere  or  between  different 
positions  in  the  storm.  Such  a model  has  been  preliminarily  developed  by 
Perlmutter  and  Frost  [2-21]-,  however,  this  model  requires  further  perfection 
and  the  coherence  function  associated  with  thunderstorms  must  be  developed  to 
permit  its  use. 

An  alternative  explanation  of  why  the  data  of  Goff  does  not  contain 
downbursts  of  the  intensity  reported  in  [2-16  and  2-17]  is  that  the  20 
thunderstorms  investigated  may  not  contain  a "spearhead  echo"  type  storm 
[2-16]  or  the  measured  data  may  not  encompass  the  downdraft  portion  of  the 
storm.  Fujita  [2-16]  defines  a downburst  as  = 3.6  m s ^ at  the  300  m 
level.  The  averaging  time  related  to  this  wind  speed  is  not  specified,  and 
it  is  not  known  whether  this  represents  a peak  gust  or  a value  averaged  over 
some  interval  of  time  which  is  undoubtedly  less  than  10  sec.  The  data  of 
Goff  [2-1]  do  not  indicate  any  values  of  that  equal  or  exceed  the  3.6  m s 

definition  of  a downburst.  There  are  a few  values,  however,  that  do  approach 
the  3.6  m s"^  level;  for  example,  thunderstorm  cases  8 and  11  shown  in  Table 
2-5  which  lists  the  maximum  and  minimum  values  of  recorded  at  the  300  m 

level  in  the  data  utilized  to  construct  the  thunderstorm  wind  shear  model 
given  in  this  report.  The  fact  that  no  downbursts  are  recorded  is  not  sur- 
prising, however,  because  obviously  the  chance  of  a downburst  being  directly 

over  the  tower  the  instant  of  maximum  intensity  is  extremely  small. 
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TABLE  2-5 


MAXIMUM  AND  MINIMUM  VERTICAL  VELOCITIES  AT  300  m 
LEVEL  IN  DATA  OF  GOFF  [2-1] 


Case 

W (m  s'b 

max 

W (m  s"\l 

min 

Case 

W (m  s”^) 

^max 

W (m  s"^) 

min 

1 

1.6 

-2.5 

11 

1 .8 

-3.0 

2 

3.0 

-1.2 

12 

2.0 

-0.9 

3 

2.6 

-1.1 

13 

3.2 

-0.7 

4 

3.0 

0.7 

14 

2.1 

-1 .0 

5 

3.0 

-1.3 

15 

2.0 

-0.5 

6 

3.7 

-1.5 

16 

2.0 

-1 .0 

7 

1 .2 

-2.1 

17 

1.1 

-0.5 

8 

4.0 

-3.0 

18 

1 .0 

-1.0 

9 

4.1 

-0.1 

19 

2.0 

-2.0 

10 

3.0 

-1 .2 

20 

0.2 

-1.1 

Fujita  [2-22]  points  out  the  probability  of  an  airport  being  under  the 
influence  of  a spearhead  echo  is  very  low  and  is  probably  less  than  2 percent 
of  the  thunderstorm  probability.  Moreover,  he  notes  that  the  location  of 
aviation  hazards  for  the  extreme  downburst  is  limited  to  only  a fraction  of 
the  spearhead  echo  area.  With  this  in  mind  it  is  relatively  easy  to  believe 
that  the  data  of  Goff  do  contain  a few  storms  of  the  magnitude  approaching 
those  types  defined  as  spearhead  echo  by  Fujita.  Moreover,  since  several  of 
the  storms  clearly  record  downdrafts  approaching  the  magnitude  of  downbursts 
for  10  sec  averaged  wind  speed,  it  is  envisioned  that  significantly  higher 
downward  gust  for  a shorter,  say,  3 sec  average  are  contained  in  the  original 
data . 

The  data  set  illustrated  in  Tables  2-3  and  2-4  is  being  punched  on  cards 
and  will  be  included  with  thunderstorm  data  now  stored  on  magnetic  tape.  It 
is  believed,  however,  that  the  current  data  set  utilizing  the  appropriate 
turbulent  intensities  will  provide  a valid  simulation  of  thunderstorms  for 
flight/hazard  simulation  studies.  The  larger  selection  of  thunderstorms,  20 
in  number,  enables  a manned  flight  simulation  study  to  provide  the  pilot  with 
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several  different  thunderstorm  situations  that  do  not  duplicate  those  that  he 
has  negotiated  in  previous  tests.  Thus,  he  cannot  learn  a given  storm.  By 
selecting  the  same  sequence  of  thunderstorms,  a second  pilot  involved  in  the 
same  test  program  can  be  exposed  to  the  identical  test  pattern  utilized  for 
the  first  pilot  with  neither  pilot  flying  the  same  storm  more  than  once. 

Finally,  it  is  called  to  the  attention  of  the  user,  that  in  carrying  out 
a simulation  where  avionics  using  ground  wind  speeds  as  inputs  are  being 
studied,  the  storm  must  be  considered  to  be  passing  over  the  airport  or 
anemometer  from  which  the  ground  wind  is  being  determined  at  the  speed  of  the 
gust  front,  W^.  Therefore,  the  length  of  the  thunderstorm  record  used  in  the 
simulation  must  be  sufficiently  long  that  the  location  of  the  assumed 
anemometer  at  the  airport  has  not  moved  out  of  the  range  of  data  set  during 
the  time  taken  by  the  aircraft  to  complete  its  approach. 
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SECTION  3.0 


NEUTRAL  AND  STABLE  BOUNDARY  LAYERS 

The  data  set  utilized  in  formulating  a mathematical  model  for  wind 
shear  or  wind  fields  associated  with  neutral  and  stable  boundary  layers 
relies  mainly  on  the  work  of  Clarke  and  Hess  [3-1].  Although  numerous 
analytical  models  for  boundary  layers,  both  under  neutral  and  stable  condi- 
tions were  available  in  the  literature,  these  were  in  general  based  on  the 
assumption  of  constant  stress  which  is  valid  only  to  the  first  50  or  100  m 
of  the  atmospheric  boundary  layer.  Other  analytical  models  which  did  not 
evoke  this  assumption  and  included  the  influence  of  turbulence  of  the  atmo- 
spheric boundary  layer  were  generally  highly  mathematical  and  required 
numerical  solution  with  a computer.  Therefore,  the  data  set  of  [3-1],  which 
are  presented  in  the  form  of  contour  maps  of  constant  longitudinal  and 
lateral  wind  speeds  as  a function  of  height,  z,  and  the  stability  of  the 
atmosphere  expressed  by  the  stability  parameter,  u,  allowed  tabulation  of 
the  wind  fields  on  a grid  system.  These  tabulated  data  were  then  coupled 
with  a table  lookup  computer  program  to  provide  the  fast  time  wind  speed 
model  for  flight  hazard  simulation. 

None  of  the  data  sets  or  mathematical  models  available  in  the  literature 
for  boundary  layers  (see  Reference  [3-2])  incorporate  horizontal  variation 
of  wind  speed.  Hence,  the  wind  shear  models  for  neutral  and  stable  boundary 
layers  reported  herein  depend  only  on  the  vertical  scale,  z.  It  should  be 
borne  in  mind,  however,  that  terrain  features  indigenous  to  each  landing 
site  can  influence  the  spatial  gradient  of  the  wind  speed  horizontally  and 
also  introduce  a possible  vertical  wind  speed  component.  Fortunately,  most 
airport  Lerrain  is  relatively  flat  and  not  surrounded  by  high  mountains  or 
sharp  cliffs  which  create  flow  disturbances,  both  vertically  and  horizontally 
in  the  wind.  Therefore,  the  wind  shear  model  described  in  the  following 
section  is  expected  to  present  a valid  simulation  of  wind  shear  in  terminal 
areas  of  the  typical  airport.  However,  indiscriminate  application  of  the 
model  is  cautioned  against  until  further  studies  of  wind  fields  over  irregular 
terrain  become  available. 
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3.1  Data  Source 


Boundary  layer  data  are  presented  in  the  form  of  horizontal  and  lateral 
wind  vector  components.  The  boundary  layers  are  described  over  flat  uniform 
homogeneous  terrain  and  therefore  have  no  vertical  component  of  wind  speed 
or  dependence  on  the  x-coordinate.  Of  course,  when  turbulence  fluctuations 
are  superimposed  as  described  in  Section  3.4,  a vertical  component  can 
occur.  The  wind  speed  profile  is  thus  dependent  only  on  the  height,  z,  and 
on  the  stability  of  the  atmosphere  characterized  by  the  stability  parameter, 
u.  The  data  set  used  in  the  boundary  layer  model  is  based  on  the  results  of 
the  extensive  wind  speed  measurements  reported  in  Reference  [3-1].  Also, 
the  influence  of  baroclinicity  is  neglected  and  is  justified  for  this  data 
in  Reference  [3-1]. 

Wind  speeds  and  temperatures  were  measured  over  a very  flat,  smooth 
surface  having  a mean  surface  roughness  parameter,  z^,  of  approximately  3.5 
mm.  Hourly  double  theodolite  observations  of  pilot  balloons  were  taken  at 
five  stations  for  40  days.  Micrometeorological  observations  were  taken  at 
two  of  the  five  stations.  Wind  profiles  were  measured  at  0.5,  1,  2,  4,  3 
and  16  m while  temperature  differences  were  measured  at  one  station  between 
1 and  2 m,  2 and  4 m.  Values  of  u*  were  estimated  by  the  drag  coefficient 
method  and  values  of  surface  heat  flux  by  means  of  temperature  and  wind 
profiles.  Surface  pressures  were  measured  at  each  of  the  five  stations  and 
radiosonde  measurements  were  made  every  three  hours  at  the  central  station. 

In  addition,  three  hourly  surface  pressure  and  temperature  data  from  14 
stations  at  distances  up  to  350  km  away  were  used  to  augment  the  data  obtained 
by  the  research  expedition. 

Geostrophic  winds  were  determined  from  the  data  by  using  second  order 
curve  fitting  procedures  on  the  pressure  data.  Thermal  winds  were  determined 
using  temperatures  from  a surface  network  reported  every  three  hours  and 
from  a temperature  field  measured  each  day  at  1500  hours.  Interpolation 
from  the  three-hour  data  set  were  primarily  used  in  assessing  baroclinicity 
which  was  shown  to  be  small.  A second  set  of  thermal  wind  estimates  based 
on  twice  daily  radiosonde  data  from  a network  of  five  stations  were  also 
made.  Reference  [3-1]  reports  that  surface  geostrophic  winds  were  well 
determined  as  evidenced  by  the  high  correlation  (93  percent)  with  observed 
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winds;  however,  no  claim  of  high  accuracy  in  estimating  thermal  winds  is 
made  in  the  report. 

These  data  thus  reduced  were  reported  by  Clark  and  Hess  [3-1]  in  the 

form  of  contour  maps  of  winds  as  a function  of  dimensionless  height,  z = 

zf/u*  , and  of  the  stability  parameter,  y = KU*/fL'.  The  contours  of  the 

map  are  lines  of  constant  AW  = W (z  = 0.15)  - W (z)  and  AW  = W (z)  - 
/s  ^ ^ y y 

Wy(z  = 0.15),  respectively,  where  W is  the  dimensionless  wind  speed  W/u*. 

The  symbol  f is  the  Coriolis  parameter  (here,  treated  positive  in  both 

hemispheres  with  a right-hand  coordinate  system  implied  in  the  Northern 

hemisphere).  For  use  in  this  report,  the  data  were  converted  to  a 34  x 11 

grid  as  illustrated  in  Figure  3-1.  The  data  were  then  stored  on  magnetic 

tape  and  a computer  look-up  routine  developed  for  interpolating  the  values 

of  W and  W for  given  values  of  z and  y.  The  tabulated  data  are  given  in 
X y 

Appendix  3A,  illustrative  wind  speed  profiles  are  given  in  Appendix  3B,  and 
a computer  program  for  looking  up  and  interpolating  the  data  is  given  in 
Appendix  3C. 

3.2  Presentation  of  Da :a 

The  data  are  presented  in  a right-hand  coordinate  system  with  W^ 
positive  in  the  direction  from  left  to  right  and  W^  positive  into  the  plane 
of  the  paper. 

The  values  of  y range  from  -333.34  to  216.67  in  increments  of  Ay  = 

16.67.  This  unusual  increment  size  was  chosen  for  convenience  in  extracting 
the  data  from  the  contour  plots.  Values  of  z range  from  0.001  to  0.15  in 
increments  of  Az  = 0.0149. 

Because  of  the  similarity  and  scaling  laws  used  in  reduction  of  the 
data,  values  of  wind  below  z = 0.001  are  not  given.  To  establish  the  profile 
from  z = 0.001  to  zero  the  log-linear  wind  speed  profile  was  used  (for  a 
description  of  the  log-linear  wind  speed  profile  see  Reference  [3-2]). 

Values  of  ^^j^(2)  below  z = 0.001  are  determined  from 

W^(z)  = 1 {£n(Ro  z + 1)  + 4.5  zu/k);  y > 0,  0 s z < 0.001  (3-1) 

where  Ro  is  the  Rossby  number  defined  as  Ro  = u^/fz^.  The  variable  z^  is  the 
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empirically  determined  surface  roughness.  Typical  values  of  z are  given  in 
Reference  [3-2].  The  value  of  W is  zero  for  0 < z < 0.001.  Introducing 
z = 0.001  into  the  equation  gives  the  value  of  wind  speed  to  which  all  the 
tabulated  longitudinal  wind  speeds  are  referenced. 

3.3  Mean  Mind  Speeds 

This  section  considers  values  of  the  mean  wind  fields  computed  from  the 
mathematical  model  of  the  neutral  and  the  stable  boundary  layers.  Although 
the  data  include  the  range  of  unstable  condition,  i.e.,  -334.34  < p < 0, 
strong  wind  shear  is  not  normally  associated  with  unstable  boundary  layers 
and  hence  no  values  of  winds  for  this  range  of  y are  given. 

The  wind  data  are  averaged  over  a 10  minute  period  or  longer  and  thus 
represent  mean  values.  A model  of  turbulence  for  the  neutral  and  stable 
boundary  layer  is  given  in  Section  3.4  which  is  superimposed  to  give  the 
random  instantaneous  wind  speed. 

The  grid  system  used  in  storing  the  data  is  numbered  from  the  bottom 
left-hand  corner.  The  numbers  increase  from  left  to  right  in  the  direction 
of  increasing  y and  from  bottom  to  top  in  the  direction  of  increasing  z. 

Thus  column  1 corresponds  to  y = -333.34  and  column  34  corresponds  to 
y = 216.67  where  row  1 corresponds  to  z = 0.001  and  row  11  corresponds  to 
z = 0.15. 

3.3.1  Tables  of  Mind  Speed 

Tables  of  the  longitudinal  and  lateral  wind  speeds  W (z)  and  W (z)  for 

X y 

neutral  and  stable  boundary  layers,  i.e.,  y ^ 0,  are  given  in  Appendix  3A. 
These  values  are  computed  from  the  tabulated  wind  differences  from  the 
relationship 

fl^(z)  = aQ^(z  = 0.001)  - AW^(i)  + W^(z  = 0.001)  (3-2) 

and 

Wy(z)  = AWy(2)  - AQy(2  = 0.001)  (3-3) 
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where  from  Equation  3-1  evaluated  at  z = 0.001, 

W^(z  = 0.001)  = K'^[)ln(0.001  Ro  + 1)  + 0.01125  y]  (3-4) 

3.3.2  Illustrations 

The  longitudinal  and  lateral  wind  speeds  which  would  be  encountered 
along  a 3°  glide  slope  are  plotted  for  various  stability  conditions  and 
Rossby  numbers  in  Appendix  3B.  Height  is  expressed  in  meters,  m,  on  the 
vertical  scale  and  wind  speeds  are  expressed  in  meters  per  second,  m s”\  on 
the  horizontal  axis. 


3.3.3  Computer  Program 

A computer  program  has  been  written  which  computes  with  inputs  of  height, 

2,  and  stability  parameter,  y,  the  longitudinal,  W , lateral,  W , wind  speeds 

X y 

and  the  wind  gradients  9W  /9z  and  9W  /9z.  The  computer  program  also  requires 

X y 

as  input  the  friction  velocity,  u^^,  the  Coriolis  parameter,  f,  and  surface 

roughness,  z . All  velocities  are  input  and  output  in  m s"^  and  lengths  in 

° -1 
meters,  f,  is  introduced  in  s . 

This  computer  code  can  be  used  as  a direct  subroutine  input  to  existing 
computer  programs  for  flight  simulators  or  computer  programs  of  airplane 
motion  in  variable  wind  fields.  The  option  of  superimposing  turbulence  on 
the  mean  wind  speed  is  available  by  appropriate  specification  of  control 
variables.  A complete  description  of  the  computer  program  is  given  in 
Appendix  3C. 


3.4  Turbulence  Model 

A turbulence  model  has  been  developed  for  use  with  the  neutral  and 
stable  boundary  layer  data  which  employs  turbulence  kinetic  energy  spectra 
developed  by  Kaimal  [3-3],  as  shown  in  Figure  3-2.  Kaimal  gives  the  following 
functional  form  for  the  spectra. 


n4>  (n) 
a 


0.164  h/Hq 
1 + 0.164(T)/nQ)^^^ 


(3-5) 
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1.0 


10 


100 


n/Ho 


son  of  Measured  Turbulence  Data  with 
m 3-5  [3-3] 
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2 -1 

(J)  (n)  Power  spectral  density  distribution  in  a direction,  m s” 

^ -1 
n Cyclic  frequency,  s 

0 Turbulence  intensity  or  root-mean-square  fluctuation  in 

a direction,  m s 

n Reduced  frequency  nz/W^(z) 

z Height  above  ground,  m 

W^(z)  Mean  wind  speed  at  the  height  z,  m s"^ 

a Designates  component  of  velocity  fluctuation  (either  w , w 

X y 

or  w^) 


The  influence  of  atmospheric  stability  enters  through  the  variable 
which  is  a characteristic  reduced  frequency  and  is  a function  of  Richardson's 
number  as  shown  in  Figure  3-3.  These  values  of  are  for  the  range 
0.04  < Ri  5 0.20. 

Values  of  for  the  neutral  boundary  layer  where  Ri  = 0 are  recommended 
in  Reference  [3-4]  as 


= 0.01444;  (n^)„  = 0.0265; 
X y 

Thus  for  the  complete  range  of  0 < 
used  in  this  report  are 


- 0.0962  (3-6) 

Ri  < 0.2  the  approximate  relationships  are 


'0.0144 

0 

J 

o 

X 

li 

0.5  Ri 

0.029 

fO.0265 

0 

J 

o 

II 

[l.5  Ri 

0.018 

0.0962 

0 

2.8  Ri 

0.035 

The  relationship 

between  Ri 

Ri  = (zu/i 

<)D  + 

4.5  zy/ic]  ^ 

< Ri  < 0.029 
S Ri  < 0.2 

< Ri  < 0.018 
i Ri  < 0.2 

S Ri  < 0.035 

< Ri  < 0.2 

and  u is 


(3-7) 

(3-8) 

(3-9) 


(3-10) 
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Figure  3-3  Correlation  of  Hq  with  Richardson  Number  for  Stable 
Boundary  Layers 


The  turbulence  spectra  data  of  Kaimal  [3-3]  given  above  appear  to  be  the  most 
comprehensive  presentation  of  atmospheric  spectra  associated  with  the  stable 
boundary  layer  and  are  taken  herein  to  represent  the  state  of  the  art. 


To  utilize  Equation  3-5  a value  of  turbulence  intensity  must  be  deter- 
mined. Figure  3-4,  taken  from  Reference  [3-5],  gives  the  turbulence  intensity 
of  the  vertical  wind  speed  component,  nondimensional ized  with  u*,  plotted 

as  a function  of  nondimensional  altitude,  zp/k.  For  neutral  conditions, 
zy/K  = 0 and  the  ratio  of  vertical  turbulence  intensity,  to  the  friction 

velocity,  u*,  becomes  1.3. 


At  zy/ic  = 1.22,  the  turbulence  intensity  vanishes  as  the  atmospheric 
boundary  layer  becomes  so  stable  that  essentially  laminar  flow  is  achieved. 
In  general,  o^^/u*  decreases  with  increasing  stability. 

A curve  fit  of  Figure  3-4  between  0 < zy/K  < 1.22  is 


a,  fl  .3  - 0.13(zy/K)°'^ 
w 

z 

u*o  [6*49  - 5.32(zy/K) 


0 ^ zy/K  <1.0 
1 .0  < zy/K  <1.22 


(3-11) 


The  value  of  can  thus  be  determined  for  given  values  of  u*,  y and  z. 

No  satisfactory  mathematical  description  of  how  and  vary  with 

large  scale  surface  features  nor  how  they  vary  with  atmospheric  stability  is 
available.  Barr,  et  al.  [3-5]  propose  that  the  ratio  be  treated  as  a 

function  of  altitude  according  to  the  following 


[0.177  + 0.832  z/z.] 
, 1-0 


z < z. 


z > z. 


(3-12) 


z.  = 300  m or  1000  ft.  Reference  [3-5]  also  proposes  that  Oy,^/oy,^  = 

This  relationship  does  not  result  in  satisfactory  agreement  with  Equation 
3-5,  however,  and  the  relationship 


a /o 

Wy 


f[0.583  + 0.417  z/z.]"*^*® 

1.0 


z < z. 
1 


z > z. 


(3-13) 
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Figure  3-4  Vertical  Turbulence  Intensity 


. is  proposed  by  Frost  [3-4].  Equation  3-13  is  developed  identical  to  Equation 

3-12  but  assumes  that  is  less  than  near  the  ground.  The  assumption 

y X 

is  that  for  neutral  conditions  = 2. 6/2. 0/1. 3 which  is  consistent 

with  a number  of  reported  results. 

Values  of  and  Oy^^/crw^  plotted  in  Figure  3-5  as  a function  of 

dimensionless  height  to  facilitate  computation  of  these  values.  Inspection 
of  Equations  3-11,  3-12  and  3-13  show  that  for  the  neutral  atmosphere 


2.6  u*. 


2.0  u*  and  a 
* w. 


1.3  u* 


(3-14) 


Figure  3-6  shows  longitudinal  turbulence  spectra  for  varying  degrees  of 
stability.  Turbulence  in  the  atmosphere  is  observed  to  decrease  with 
increasing  stability.  This  is  most  evident  at  the  low  frequency  range  and  is 
to  be  expected  since  the  thermal  effects  which  create  stable  boundary  layers 
depress  large  scale  turbulent  motion.  These  spectra,  however,  are  reported 
to  be  influenced  by  surface  terrain  features,  particularly  in  the  longi- 
tudinal and  lateral  components.  No  mathematical  models  which  account  for 
terrain  effects  are  available.  The  spectra  having  been  measured  for  rela- 
tively homogeneous  terrain  which  is  characteristic  of  the  majority  of 
airports  should  give  a valid  representation  of  most  terminal  areas.  For  air- 
ports located  near  unusual  terrain  features  such  as  mountains  or  cliffs, 
fluctuations  in  the  longitudinal  and  lateral  components  of  the  wind  may, 
however,  be  higher  than  simulated  by  the  proposed  turbulence  model. 

The  spectra  represented  by  Equation  3-5  do  not  have  a rational  form  and, 
consequently,  are  difficult  to  use  in  turbulence  simulation  schemes  [3-6]. 

To  overcome  this  difficulty,  the  spectral  data  in  Figure  3-2  were  adjusted  to 
fit  a modified  Dryden  spectrum  having  the  functional  form 

Pi.(^.  = (3-15 

1 + 

The  constants  C-j  and  C2  are  adjusted  such  that  the  modified  Dryden  curve  fits 
Kaimal's  data  [3-3].  Values  of  C-j  = 0.1580  and  C2  = 0.0694  were  determined 
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Figure  3-5  Longitudinal  and  Lateral  Turbulence  Intensity 
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Figure  3-6  Influence  of  Stability  on  Turbulence  Energy  Spectrum 


to  give  the  best  curve  fit.  Thus,  the  power  spectral  density  function  chosen 
for  the  turbulence  simulation  model 


0.1580  h/tIq 
1 + 0.0694(n/riQ)^ 


(3-16) 


Figure  3-7  compares  the  modified  Dryden  spectra  with  Equation  3-16.  Although 
this  curve  fit  does  not  provide  a good  representation  of  the  higher  frequency 
spectra,  it  does  fit  the  data  quite  well  in  the  lower  frequency  range  which 
is  expected  to  have  the  most  significant  influence  on  the  flight  of  aircraft. 

This  turbulence  model  utilizing  the  turbulence  spectra  given  by  Equation 
3-16  and  the  z transformation  technique  (see  Neuman  and  Foster  [3-7])  has 
been  developed  and  can  be  coupled  with  the  mean  wind  fields  for  the  stable 
and  neutral  boundary  layers  to  give  a random  fluctuating  field.  The  model, 
however,  is  linear  and  results  in  a Gaussian  distribution  of  the  wind  speeds 
in  the  atmosphere  which  introduces  a small  error  into  the  simulation  (see 
Reeves,  et  al . [3-6]). 

To  illustrate  the  influence  of  turbulent  fluctuations  on  the  velocity 
profile  experienced  by  an  aircraft  during  landing  in  a stable  boundary  layer. 
Figures  3-8  through  3-13  have  been  prepared.  These  figures  illustrate  stable 
boundary  layer  wind  speed  profiles  seen  by  an  aircraft  on  a 3°  glide  slope 
with  turbulence,  computed  by  the  simulation  technique,  superimposed.  The 
inertial  aircraft  velocity  is  64  m s”^  which  corresponds  to  a sink  rate  of 
3.35  m s~^ . The  time  increment  used  in  the  turbulence  simulation  was  taken 
as  0.15  sec  which  results  in  a turbulent  fluctuation  being  superimposed  at  Az 
increments  of  0.35  m along  the  flight  path.  The  figures  show  some  interest- 
ing results  relative  to  the  stability  of  the  boundary  layer.  In  these 
figures,  u*  has  been  held  constant  and  y is  increased  in  sequential  order. 

The  first  velocity  profile  (Figure  3-8)  corresponds  to  the  neutrally  stable 
layer  where  y is  essentially  zero.  One  observes  that  the  aircraft  encounters 
turbulence  from  the  500  m level  to  the  ground  and  that  the  turning  of  the 
boundary  layer  is  reasonably  small  after  100  m.  The  second  figure  (Figure 
3-9)  illustrates  a somewhat  higher  level  of  stability,  y = 25.  The  longi- 
tudinal wind  speed  is  larger  because,  as  noted  earlier,  the  computed 
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Figure  3-7  Comparison  of  the  Modified  Oryden  Spectrum  with  Kaimal's  Emperical  Curve  Fit  [3-3] 
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Figure  3-8  Near  Neutral  Boundary  Layer  with  Turbulence  Superimposed 
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Figure  3-11  Stable  Boundary  Layer  with  Turbulence  Superimposed 
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Boundary  Layer  with  Turbulence  Superimposed 


profiles  are  based  on  the  same  friction  velocity,  u*,  and  surface  roughness, 
Zq.  The  interesting  feature  of  this  figure  is  that  at  large  values  of  z, 
i.e.,  where  z/L'  is  large,  the  mechanical  turbulence  is  damped  out  by  buoyancy 
induced  turbulence,  however,  as  z/L'  becomes  small  due  to  the  aircraft 
approaching  the  ground,  the  situation  is  reached  where  the  atmospheric 
boundary  layer  returns  to  a neutral  condition  and  mechanical  turbulence 
dominates  the  buoyancy  induced  turbulence.  In  Figure  3-9  this  occurs  at 
approximately  125  m.  The  airplane  thus  flies  from  a region  of  rather  quiet 
flow  to  a sudden  and  rapidly  increasing  turbulent  intensity.  The  turbulence 
intensity  becomes  larger  as  y increases,  however,  this  is  not  due  to  the 
effects  of  stability  but  due  to  higher  mean  velocities  resulting  from  a 
constant  u*.  For  consecutive  figures,  as  the  stability  of  the  boundary  layer 
increases  one  observes  that  the  transition  from  the  essentially  laminar  flow 
to  the  highly  turbulent  region  approaches  lower  and  lower  levels.  It  should 
also  be  noted  from  this  sequence  of  figures  that  the  strongest  directional 
shear  occurs  at  intermediate  values  of  p,  whereas,  the  strongest  linear  shear 
occurs  at  the  larger  values  of  y. 

An  alternate  way  of  computing  the  profiles  would  be  to  assign  the  same 
wind  speed  at  a certain  level  and  adjust  u*  accordingly.  Physically  this 
corresponds  to  conditions  encountered  at  an  airport  where  the  wind  speed 
measured  with  an  anemometer  at  a fixed  height,  say,  at  the  10  m level  would 
record  identical  values  but  the  wind  shear  would  be  appreciably  different  due 
to  the  associated  stability  conditions.  Figures  3-14  through  3-19  show  com- 
puted profiles  having  a common  value  of  10  m s"^  wind  speed  at  the  10  m 
level.  The  value  of  u^  used  in  the  computation  is  related  to  the  stability 
parameter,  y,  by 

u*  = 0.58  - 0.0016  y 

The  same  characteristics  of  the  wind  speed  profiles  are  observed,  how- 
ever, the  higher  wind  speeds  now  occur  at  intermediate  values  of  y rather 
than  at  the  higher  values  as  in  the  sequence  of  Figures  3-8  through  3-13.  In 
turn,  one  observes  that  the  directional  shear  is  also  largest  at  intermediate 
values  of  y which  is  in  correspondence  with  the  former  sequence  of  figures. 
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Fiqure  3-14  Near  Neutral  Boundary  Layer  with  Turbulence  Superimposed 
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Figure  3-16  Stable  Boundary  Layer  with  Turbulence  Superimposed 
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Figure  3-17  Stable  Boundary  Layer  with  Turbulence  Superimposed 
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Figure  3-19  Stable  Boundary  Layer  with  Turbulence  Superimposed 


3.5  Statistical  Model  Relative  to  Risk  of  Exceedance 


The  wind  shear  associated  with  the  stable  and  neutral  boundary  layers 
are  dependent  upon  the  variables,  p and  u*,  or  the  mean  wind  speed,  W^,  since 
they  are  directly  related.  A statistical  description  of  these  two  parameters 
will  enable  the  user  of  the  wind  shear  model  to  provide  an  estimate  of  the 
probability  and  frequency  of  exceeding  a prescribed  value  of  wind  shear. 

In  order  to  establish  this  statistical  description  some  estimate  of  the 
probability  of  a given  stability  condition  occurring  simultaneously  with  a 
given  value  of  u*  is  needed.  An  approximate  analysis  to  achieve  this  goal  is 
proposed  as  follows.  From  statistical  theory  the  probability  of  p and  u* 
occurring  simultaneously  is  equal  to  the  probability  of  p given  u*, 

P(p  > Pp/u*),  times  the  probability  of  u*,  P(u*  > u*p),  [3-8].  This  is 
referred  to  as  a conditional  probability. 


P(p  > Pp;  u*  > u*p)  = P(p  > Pp/u*)  P(u*  > u*p)  (3-17) 

The  probability  of  p given  a wind  speed,  W^,  which  can  be  directly 
related  to  u*,  has  been  estimated  by  Barr,  et  al . [3-5].  They  give  the 
probability  of  Richardson  number,  Ri,  for  a given  wind  speed,  W^,  where  both 
and  Ri  are  evaluated  at  6.1  m,  shown  in  Figure  3-20.  Note  that  these 
curves  are  highly  interpolated  and  based  on  only  two  sets  of  data.  The 
Richardson  number,  Ri , is  related  to  p by  the  relationship 


Ri 


p/u* 

</2f  + 4.5  p/u* 
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0 < Ri  < 0.18 


Ri  > 0.18 


and  is  related  to  u*  and  p by 
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Figure  3-20(a)  Cumulative  Percent  Frequency  of  Occurrence  of  Ri  at  Given  Wind  Speeds, 
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Figure  3-20(b)  Cumulative  Percent  Frequency  of  Occurrence  of  Ri  at  Given  Wind  Speeds, 
6 to  12  m s"^ (14  to  28  mph) 


13-14  m s 

(30-32  mph)  ”z=6.1m  (20  ft) 


U) 


x: 

a. 

e £ 


CL 

e 


x: 

CL 

e 


L 


_L.  __ 


o 

o 


o 

CO 


o 

vO 


o 

'4- 


o 

Csi 


J 

o 


m m 
^ ro 
I I 
•<r  cN 
iH  cn 


aDuajjriDDo  jo  iCouanbajj  ^Udojad  SAja^xnainQ 


71 


Figure  3-20(c)  Cumulative  Percent  Frequency  of  Occurrence  of  Ri  at  Given  Wind  Speeds, 
12  to  15  m s”\28  to  32  mph) 


Assigning  values  of  z = 6.1  m,  = 0.01  m,  f = 10"^  s"^,  k = 0.4,  Equations 
3-18  and  3-19  have  been  plotted  in  Figures  3-21  and  3-22  for  convenient 
reference. 

Using  the  above  results  the  designer  may  select  a value  of  y and  a value 
of  u*  which  gives  a wind  shear  condition  in  which  he  is  interested.  The 
probability  of  y given  u*  is  then  obtained  directly  from  the  figure. 

The  probability  of  u*  can  be  determined  from  the  distribution  of  mean 
wind  speeds  given  by  Frost  [3-4]  as  interpreted  from  Justus,  et  al.  [3-9], 
This  reference  shows  that  the  probability  of  a mean  wind  speed  greater  than  a 
prescribed  value  Wp  occurring  during  the  year  is  described  by  a Rayleigh 
distribution,  i.e.. 


P(W  > Wp) 


(3-20) 


Reference  [3-4]  shows  that  this  equation  is  a general  representation  of  wind 
speeds  measured  at  138  airport  sites  throughout  the  United  States  where  an 
average  value  of  the  scale  factor,  c,  is  4 n s"^  ± 0.9  n s~\  Since  u*  is 
related  to  mean  wind  speed  by  Equation  3-19,  the  probability  of  u*  is 
directly  related  to  the  probability  of  U^. 

With  the  information  provided  above,  the  user  of  the  neutral  and  stable 
boundary  layer  wind  shear  model  can  estimate  what  the  probability  of  a given 
wind  shear  occurring  over  an  expected  lifetime  is  by  utilizing  the  prob- 
ability of  y given  u*  from  Figure  3-20  and  multiplying  that  by  the 
probability  of  u*  given  by  Equation  3-19.  To  illustrate  the  use  of  the  above 
model,  consider  the  following  example.  A user  wishes  to  determine  the 
probability  of  y being  equal  to  200  given  that  u,^  is  equal  to  0.5  m s"^ . 

From  Figure  3-21  find  Ri  = 0.164  corresponding  to  y = 200  and  u*  = 0.5  m s~\ 
i.e.,  y/u*  = 400.  From  Figure  3-22  determine  W /u*  = 23  and  therefore  W = 
11.5  m s for  the  given  values  of  u*  and  y.  From  Figure  3-20  the  probability 
of  Ri  ^ to  0.16  at  11.5  m s'^  is  1.0  - P(Ri  ^ 0.16;  W^  = 11.5  m s"^  = 0.12, 
From  Equation  3-20  the  probability  of  W ^ 11.5  m s'^  is  P{W  > 11.5  m s'^)  = 

X X 

0.0003.  The  probability  of  both  y > 200  and  W^  i 11.5  occurring  in  a given 
year  is  from  Equation  3-17  equal  to  0.003%. 
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3-22 


To  establish  the  risk  of  exceeding  a prescribed  wind  shear  in  a given 
period  of  time  we  proceed  as  follows.  The  probability,  p,  that  u will  be 
exceeded  in  any  one  year  is  related  to  the  return  period  T by  the  relation- 
ship Tp  = 1.  The  probability  that  a value  less  than  or  equal  to  p will  occur 
in  one  year  is  q = 1 - p.  In  a period  of  N years  the  probability  Q that  a 
value  less  than  or  equal  to  y will  occur  is 

Q = q^^  (3-21) 

The  probability  or  risk  that  a value  greater  than  y will  occur  at  least  once 
in  a period  of  N years  is 


= 1 - (1  - p)^  (3-22) 

Hence  continuing  the  foregoing  example,  the  probability  that  y = 200  and 
u,,  = 0.5  m s ^ will  occur  at  least  once  in  a period  of  25  years  is  from 
Equation  3-22  equal  to  3.4%. 

To  give  the  reader  some  feel  for  the  nature  of  the  wind  shear  associated 
with  a given  probability  of  occurrence.  Table  3-1  has  been  prepared.  This 
table  gives  the  risk  of  exceedance  associated  with  each  of  the  wind  shear 
conditions  depicted  in  Figures  3-8  through  3-18. 

3.6  Application  of  Wind  Shear  Models  for  Neutral  and  Stable 
Boundary  Layers 

Application  of  the  wind  shear  model  for  the  neutral  and  stable  boundary 
layer  proceeds  by  first  selecting  stability  conditions  for  which  the  simula- 
tion is  to  be  carried  out.  The  parameters  required  are  the  stability 
parameter,  y,  the  friction  velocity,  u*,  and  the  surface  roughness,  z^.  The 
rule-of-thumb  for  the  effect  of  these  individual  parameters  is  that  high 
values  of  u*,  y and  all  result  in  larger  linear  shear.  On  the  other  hand, 

intermediate  values  of  y,  say,  on  the  order  of  50,  give  the  largest  direc- 
tional shear  within  the  lower  500  m of  the  atmospheric  boundary  layer. 
Inspection  of  Figures  3-8  and  3-9  clearly  illustrate  this. 
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TABLE  3-1 


RISK  OF  EXCEEDANCE  ASSOCIATED  WITH  WIND  SPEED  PROFILES 


SHOWN  IN  FIGURES  3-8 

THROUGH 

3-19 

U* 

U 

Ri 

Wx 

f(“x) 

P(Ri/W^) 

P(Ri  & W^) 

Risk  of 
Exceedance 
in  25  Yrs. 

0.50 

200 

0.18 

13.64 

8o85xl0"^ 

0o70 

6.2x10“^ 

0.16% 

0.50 

150 

0.17 

12.24 

8.61x10"^ 

0.74 

6.4x10"^ 

0.16% 

0.50 

100 

0.15 

10.83 

6.54x10"^ 

0.65 

4.3x10"^ 

1.07% 

0.50 

75 

0.14 

10.13 

1.64x10'^ 

0.64 

1.0x10"^ 

2.47% 

0.50 

50 

0.12 

9.43 

3.88x10"^ 

0.58 

2.3x10"^ 

5.59% 

0.50 

25 

0.080 

8.72 

8.61x10"^ 

0.56 

4.8x10"^ 

11.33% 

0.50 

1 

0.005 

8.05 

17.47x10"^ 

0.42 

7.3x10“^ 

16.74% 

0.25 

200 

0.35 

6.32 

8.96x10"^ 

0.72 

6.5x10'^ 

0.002% 

0.34 

150 

0.20 

6.21 

1.09x10“^ 

0.61 

6.6x10"® 

0.02% 

0.42 

IDO 

0.16 

12.11 

1.04x10"^ 

0.68 

7.1x10"® 

0.18% 

0.46 

75 

0.14 

11.41 

2.92x10“^ 

0.64 

1.9x10"^ 

0.47% 

0.50 

50 

0.12 

10.71 

7.72x10"^ 

0.62 

4.8x10"^ 

1.19% 

0.56 

12.5 

0.04 

9.65 

2.95x10"^ 

0.59 

1.7x10"® 

4.16% 

0.58 

1 

0.0042 

9.33 

4.34x10"^ 

0.54 

2.3x10"® 

5.59% 

Having  selected  the  parameters  of  interest,  one  then  determines  the  pro- 
bability that  these  conditions  will  occur  and  what  the  risk  is  of  exceeding 
the  values  that  are  planned  for  the  simulation.  If  it  is  desired  to  simulate 
severe  cases  of  wind  shear,  then  values  of  u and  u*  should  be  chosen  to  have 
large  magnitudes  which  have  correspondingly  lower  risk.  Therefore,  a low 
value  of  risk  is  acceptable.  Simulating  average  daily  conditions,  however, 
allows  a higher  risk  of  exceeding  the  prescribed  stability  conditions  to  be 
accepted.  The  percent  risk  that  the  user  of  the  wind  shear  model  is  willing 
to  accept  is  subject  to  his  own  judgment  based  on  the  engineering  application 
for  which  the  simulation  is  being  carried  out  and  on  the  consequences  asso- 
ciated with  the  test  not  being  of  sufficient  severity. 

One  must  also  bear  in  mind  that  the  risk  of  exceeding  a prescribed  value 
computed  in  Section  3.4  and  tabulated  in  Table  3-1  are  on  a per  airport  per 

76 


r 


year  basis;  consequently,  although  the  probability  of  u*  = 0.5  and  y = 200 
occurring  simultaneously  is  6.2  x 10"^^^  (see  Table  3-1),  for  10,000  airports, 
this  condition  will  occur  at  0.62  of  an  airport  or  at  least  one  airport  per 
year.  In  view  of  the  fact  that  stability  conditions  may  endure  for  one  or 
two  hours  at  an  airport  having  heavy  traffic,  several  pilots  will  be  exposed 
to  these  conditions  every  year. 

Use  of  the  turbulence  simulation  routine  to  be  superimposed  upon  the 
steady  state  winds  requires  no  additional  specification  of  atmospheric  param- 
eters, however,  the  time  step  increment  of  the  random  signal,  DT,  in 
subroutine  STB,  must  be  specified.  In  general,  its  value  is  set  equal  to  the 
time  step  used  in  the  calling  program  carrying  out  the  numerical  integration 
of  the  equations  of  motion  for  the  flight  dynamics  of  the  aircraft.  On  the 
other  hand,  it  can  be  shown  [3-7]  that 

T < 4 X IO'^HqW^/z 

but  since  this  value  varies  with  altitude  an  average  value  on  the  order  of 
0.01  sec  or  less  is  generally  recommended  for  use  with  the  turbulence  simula- 
tion routine. 

In  carrying  out  a realistic  simulation,  it  is  recommended  that  the 
turbulence  be  used  with  the  mean  wind  profile.  It  is  evident  from  Figures 
3-8  through  3-19  that  simulation  of  high  stability  conditions,  that  is,  large 
y,  is  accompanied  with  a transition  from  a rather  smooth  wind  aloft  to 
extremely  severe  turbulence  near  the  200  to  100  m range,  depending  upon  the 
stability.  This  sort  of  turbulence  phenomenon  is  expected  to  be  quite  real- 
istic judging  from  corriments  with  commercial  airline  pilots  and  from 
experimental  data  available  in  the  literature.  Super-positioning  turbulence 
on  the  wind  fields,  however,  does  require  that  a computer  be  available  for 
carrying  out  the  computations  of  the  wind  speeds.  This  does  not  impose  a 
hardship  because  normally,  a computer  is  involved  whether  the  user  is  per- 
forming a manned  flight  simulation  or  programming  a fast  time  computer 
solution. 

Without  turbulence,  however,  one  can  construct  a wind  speed  profile  for 
the  neutral  or  stable  boundary  layer  manually  from  the  tables  given  in 
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Appendix  3A.  As  an  example,  determine  the  wind  speed  profile  for  p = 50, 
u*  = 0.5  and  = 0.05.  To  establish  the  profile  one  first  computes  the 
Rossby  number  Ro  = u^/z^f  = 1 x 10^.  This  value  corresponds  to  Table  3A-1  in 
Appendix  3A.  Selecting  the  column  labeled  y = 50  gives  the  wind  speed 
profiles 


z 

M 

X 

y 

0.15 

63.4 

-32.0 

0.10 

63.4 

-27.2 

0.05 

59.3 

-17.0 

0.005 

29.5 

-4.8 

0.0001 

1 24.1 

0.0 

Dimensional  values 

y\  ^ 

are  obtained  from  z = u*z/f,  = W^u*  and  = W^u*, 

hence 

z 

W 

W 

X 

y 

-1 

-1 

m 

m s 

m s 

750 

31.7 

-16.0 

500 

31.7 

-13.6 

250 

29.6 

-9.5 

25 

14.7 

-2.4 

5 

12.1 

0 

Thus, 

wind  speed  profiles  can  be  established  very  simply. 
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SECTION  4.0 


FRONTAL  WIND  SHEAR 

A model  of  frontal  wind  shear  has  been  developed  in  exactly  the  same 
manner  as  the  thunderstorm  model.  Unfortunately,  the  number  of  data  sets 
available  for  processing  are  only  two,  which  are  shown  in  their  original 
format  in  Figures  4-1  and  4-2.  These  data  were  provided  courtesy  of 
Ms.  Judith  Stokes  and  Mr.  Jean  Lee  from  the  National  Severe  Storms 
Laboratory  in  Norman,  Oklahoma.  The  data  represent  cold  frontal  passage  in 
the  vicinity  of  the  500  m tower  and  the  reduction  of  these  data  sets  is  the 
same  as  used  by  Goff  [4-1]  for  thunderstorm  gust  front  cases  (see  Section 
2.2).  The  data  were  measured  on  January  28,  1977,  labeled  Case  1 in  this 
report,  and  on  December  10,  1976,  labeled  Case  2. 

No  other  data  which  would  provide  the  necessary  detailed  wind  speed  pro- 
files over  at  least  a two-dimensional  plane  in  space  were  found  during  this 
study.  Additional  data  which  is  on  magnetic  tape  at  the  National  Severe 
Storms  Laboratory  have  not  yet  been  processed.  No  warm  front  wind  shear  data 
from  which  a model  for  flight/hazard  simulation  studies  can  be  constructed 
appears  to  be  available. 

The  simulation  model  developed  from  the  limited  data  is  constructed  in 
the  same  manner  as  the  thunderstorm  models.  That  is,  a grid  system  is  super- 
imposed on  the  contours  of  constant  wind  speed  shown  in  Figures  4-1  and  4-2. 
In  these  two  cases  a 113  x 11  grid  is  used.  The  data  were  again  punched  on 
cards  and  stored  on  magnetic  tape.  The  only  difference  between  these  data 
and  those  of  the  thunderstorm  cases  is  the  larger  grid  having  113  columns 
rather  than  41  which  represents  a total  horizontal  distance  of  approximately 
17.2  km  in  Case  1 and  13.2  km  in  Case  2.  This  represents  approximately  30 
minutes  of  recorded  data. 

The  following  section.  Section  4.1,  describes  the  quasi-steady  wind 
speed  profile  working  data  developed  for  the  large  scale  frontal  wind  shear. 
The  wind  speed  is  referred  to  as  quasi -steady  in  view  of  the  fact  that,  as 
with  the  thunderstorm  data,  it  has  been  averaged  over  10  sec  time  intervals 
and  thus  contains  fluctuations  in  the  wind  speed  of  0.1  Hz  or  smaller.  It  is 
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Figure  4-1  Wind  Speed  for  Synoptic  Front  Cold  Air  Out-flow 
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proposed  that  the  same  turbulence  simulation  model  for  thunderstorms  be  used 
for  the  frontal  wind  shear.  Since  absolutely  no  data  was  found  which  would 
give  information  pertaining  to  turbulence  intensities  and  characteristic 
turbulent  length  scales  in  a major  frontal  flow,  the  lower  values  associated 
with  thunderstorms  are  recommended  as  input  to  the  model. 

4.1  Quasi-Steady  Wind  Speed  Profile  Grid  System 

The  wind  fields  shown  in  Figures  4-1  and  4-2  were  fit  to  a 113  x 11 
point  grid  system.  Wind  speeds  at  each  grid  point  were  tabulated  and  punched 
onto  computer  cards.  The  data  were  later  stored  on  magnetic  tape. 

The  grid  system  is  numbered  from  the  left-hand  bottom  corner.  The 
numbers  increase  from  left  to  right  in  the  positive  x-direction  and  from 
bottom  to  top  in  the  positive  z-direction  of  the  original  data.  The  wind 
speeds  are  given  in  units  of  meters  per  second,  m s”\  with  being  positive 
in  the  direction  of  frontal  motion,  being  positive  upward  and  W^  being 
positive  into  the  plane  of  the  paper.  The  frontal  speeds  of  the  two  storms 
are  9.3  m s~^  and  7.1  m s'\  respectively. 

The  wind  speeds,  W , shown  in  Figures  4-1  and  4-2  are  the  wind  speeds 
relative  to  the  storm  motion  which  are  the  values  punched  on  cards  and  stored 
on  magnetic  tape.  To  convert  to  wind  speed  relative  to  the  ground  the  frontal 
motion  must  be  added  to  W^.  The  next  subsections  describe  tables  and 
graphical  illustrations  of  the  wind  speed  profiles. 

4.1.1  Tables  of  Wind  Speed 

Tables  of  synoptic  cold  frontal  passage  wind  speeds  are  given  in 
Appendix  4A.  The  tables  have,  due  to  their  length,  been  split  into  six  parts 
covering  grid  stations  1 through  20,  21  through  40,  etc.  The  tabulated 
values  of  W^  are  values  relative  to  the  storm  motion.  The  frontal  speed  is 
given  at  the  top  of  the  table  for  converting  W^  to  its  value  relative  to 
ground. 

The  frontal  storm  case  numbers'  designation  as  1 and  2 in  this  report 
are  listed  in  the  upper  left-hand  corner. 
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Also,  listed  at  the  top  of  the  table  is  the  horizontal  length  scales  for 
the  given  wind  record.  The  horizontal  extent  of  each  data  set  varies  because 
of  the  data  reduction  procedure.  Hence,  the  length  of  field,  L,  in  kilo- 
meters and  the  horizontal  grid  spacing.  Ax,  are  specified  on  each  table.  The 
vertical  extent  of  each  field  is  taken  as  500  m with  50  m vertical  grid 
spacing. 

4.1.2  Illustrations 

Illustrations  of  the  longitudinal,  lateral,  and  vertical  wind  speeds 
encountered  along  four  3°  glide  slope  through  each  front  are  provided  in 
Appendix  4B.  The  glide  slopes  are  displaced  in  increments  of  4.3  km  for  Case 
1 and  3.3  km  for  Case  2.  Each  profile  is  the  wind  seen  by  a airplane  travel- 
ing along  the  flight  path  drawn  across  the  streamlines  as  illustrated  in  the 
first  figure  of  the  appendix.  Note  that  the  streamlines  plotted  are  relative 
to  the  speed  of  the  front  which  for  reference  purposes  is  indicated  with  the 
vertical  dashed  line  on  the  horizontal  wind  speed  profile.  The  wind  speed 
profiles  are  relative  to  the  fixed  earth  frame  of  reference. 

4.1.3  Computer  Program 

The  card  deck  or  magnetic  tape  on  which  the  frontal  wind  speeds  are 
tabulated  can  be  directly  coupled  with  the  program  for  thunderstorms  given  in 
Appendix  2C  and  requires  and  provides  identical  input/output  statements.  The 
only  modification  needed  is  to  increase  all  dimension  statements,  do-loop 
statements  and  read-in  control  integers  from  41  to  113. 

The  turbulence  option  can  also  be  called  for  if  desired.  The  turbulence 
intensities  and  and  length  scales  and  used  in  the 

program  all  have  values  as  estimated  for  thunderstorms,  however,  in  lieu  of 
better  information  on  turbulence  in  major  fronts,  these  values  can  be  used. 

4.2  Turbulence  Model 

No  data  on  turbulence  in  large  frontal  flows  were  found  in  the  litera- 
ture. It  is  therefore  suggested  that  the  model  proposed  for  the  thunderstorm 
be  employed  along  with  the  lower  values  of  the  input  parameters.  Recall  that 
data  on  turbulence  in  thunderstorms  is  also  very  limited  in  the  region  near 
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the  ground,  and  the  thunderstorm  turbulence  model  is  itself  an  approximation. 
Again,  it  is  pointed  out  that  only  the  high  frequency  components  of  the 
turbulence  (i.e.,  frequencies  greater  than  0.1  Hz)  have  been  lost  in  the  data 
since  the  quasi-steady  wind  field  given  represent  10  sec  averages.  There  is, 
however,  a very  great  need  to  obtain  turbulence  data  in  thunderstorms  and 
large  fronts  in  the  vicinity  of  the  ground  or  terminal  areas  of  an  airport. 

4.3  Statistical  Model 

Due  to  the  paucity  of  detailed  wind  speed  data  for  major  fronts,  no 
attempt  is  made  to  establish  a statistical  model  of  the  wind  fields. 

4.4  Applications  of  Frontal  Wind  Shear  Model 

The  frontal  wind  shear  model  is  applied  exactly  the  same  as  the  thunder- 
storm models  described  in  Section  2.6.  The  user,  however,  has  only  limited 
choice  of  wind  fields  to  use  in  the  simulation  since  only  two  fields  are 
available.  Keenan  [4-2]  has  used  thunderstorm  case  numbers  2 and  3 as  fronts 
since  these  two  thunderstorms  behave  in  a manner  similar  to  fronts.  However, 
they  are  extremely  short  in  length  and  are  not  long  enough  for  a full  length 
simulation  based  on  an  aircraft  approaching  from  500  m along  a 3°  glide 
slope.  Such  an  approach  requires  a horizontal  fetch  of  9.5  km.  Fortunately 
the  frontal  data  that  are  available  for  the  two  cases  of  wind  shear  studied 
extend  over  a sufficiently  long  distance  that  a sequence  of  one  or  more 
landings  or  takeoffs,  can  be  simulated  using  these  data.  Up  to  two  landings 
assuming  a 3 mile  separation  distance  can  be  simulated.  Moreover,  the  wind 
fields  encountered  along  flight  paths  separated  by  small  distances  show  quite 
dissimilar  wind  speed  profiles  as  evidenced  from  the  graphical  display  in 
Appendix  4B. 

Again,  it  is  recommended  that  turbulence  models  be  used  in  conjunction 
with  the  quasi-steady  state  wind  speeds  utilized  in  developing  the  frontal 
wind  shear  model.  It  is  anticipated  that  extreme  wind  speeds  will  be  somewhat 
averaged  out  using  10  sec  averages,  just  as  described  for  the  thunderstorm 
models,  and  therefore  to  achieve  a correct  simulation  high  gusts  should  be 
introduced  in  the  form  of  turbulence.  The  turbulence  intensities  used  in  the 
model,  however,  are  still  open  to  question,  and  it  is  proposed  that  the  lower 


extremes  of  turbulent  intensities  quoted  for  the  thunderstorm  cases  be  util- 
ized for  frontal  wind  shear. 


Much  more  data  on  fronts  is  needed  before  an  adequate  model  of  a synoptic 
front  wind  field  is  possible.  Such  data  may  be  available  at  the  National 
Severe  Storms  Laboratory  in  Norman,  Oklahoma,  but  require  further  processing. 
Also,  to  utilize  the  10  sec  average  wind  speeds  generated  by  the  NSSL  Program 
effectively,  more  must  be  determined  about  the  turbulence  levels  in  major 
fronts. 

One  additional  set  of  data  that  the  engineer  interested  in  carrying  out 
warm  front  simulation  should  be  aware  of  is  a one-dimensional  profile  compiled 
by  Keenan  [4-2].  This  model  was  developed  from  data  obtained  in  an  accident 
which  occurred  in  Tokyo  in  1966.  Table  4-1  illustrates  the  warm  front 
one-dimensional  model. 


TABLE  4-1 

WARM  FRONT  WIND  SPEED  PROFILE  [4-2];  W^  = 0 


z (m) 

W^  (m  s"^) 

Wy  (m  s'b 

6.1 

-6.6 

3.9 

30.5 

-3.6 

-2.1 

61 .0 

-4.4 

-12.1 

91.4 

0 

-20.6 

121.9 

5.4 

-30.4 

152.4 

5.9 

-33.5 

182.9 

5.8 

-32.9 

213.3 

5.8 

-32.9 

243.8 

5.7 

-32.4 

274.3 

5.6 

-31 .9 

304.8 

5.6 

-31 .4 

457.2 

5.9 

-31 .4 
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SECTION  5.0 


SUMMARY 

Mathematical  models  of  wind  fields  associated  with  thunderstorms, 
neutral  and  stable  boundary  layers,  and  storm  fronts  have  been  developed 
based  on  reported  experimental  data.  These  data  have  been  tabulated  and 
presented  in  Appendices  2A,  3A  and  4A,  respectively.  Computer  programs  have 
been  developed  which  utilize  table  lookup  and  interpolation  routines  to  com- 
pute wind  speeds  at  spatial  positions  and  under  stability  conditions  as 
called  for  by  the  main  calling  computer  program.  These  computer  programs  are 
presented  in  Appendices  2C  and  3C. 

Turbulence  simulation  techniques  have  been  developed  which  permit  the 
super-positioning  of  randoi"  fluctuations  having  the  characteristics  of 
turbulence  associated  with  the  atmospheric  phenomena  of  interest. 

The  wind  speed  models  for  thunderstorms  and  major  fronts  represent  the 
three  components  of  wind  velocity  in  a vertical  plane  500  m in  extent  and  of 
variable  lengths  in  the  horizontal  direction.  In  simulating  an  approach, 
application  of  these  models  assumes  that  the  flight  path  lies  in  the  plane  of 
the  wind  speeds  with  that  plane  being  centered  over  the  runway.  These  models, 
therefore,  give  realistic  simulation  of  wind  shear  due  to  a thunderstorm 
moving  parallel  to  the  runway  and  allow  the  flight  mechanics  computation  to 
have  departures  from  the  flight  path  in  the  vertical  direction.  Any 
departure  of  the  aircraft  out  of  the  plane  in  the  lateral  direction  is  not 
accurately  modelled.  However,  evidence  is  shown  that  the  correlation  between 
the  wind  components  in  the  x-  and  z-direction  is  large  and,  therefore,  it  is 
anticipated  that  the  same  is  true  in  the  y-direction.  Therefore  the  wind 
disturbances  encountered  by  an  aircraft  displaced  by  the  wind  a reasonable 
distance  out  of  the  plane  of  the  wind  field  will  be  nearly  similar  to  the  in- 
plane winds. 

To  simulate  a flight  path  which  is  diagonal  to  the  direction  of  the 
motion  of  the  front,  one  must  assume  that  the  wind  fields  modelled  are  two- 
dimensional  over  a very  large  lateral  scale.  This,  however,  is  not  likely 
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and  to  provide  wind  fields  which  permit  this  type  of  simulation  will  require 
considerably  more  research  into  the  three-dimensional  nature  of  thunderstorms 
and  major  fronts. 

The  turbulence  model  developed  to  accompany  the  mean  wind  field  data  are 
believed  to  be  quite  realistic  with  the  exception  that  the  value  of  the 
turbulence  parameters  in  the  frontal  case  are  vertically  unknown.  The  turbu- 
lence model  for  the  neutral  and  stable  boundary  layers  is  based  on  extensive 
data  sets  and  theoretical  analysis  and  is  believed  to  be  very  representative 
of  what  happens  in  the  actual  atmosphere.  The  thunderstorm  turbulence  model 
is  based  on  advanced  turbulence  concepts  and  does  include  large-scale  velocity 
fluctuations  of  a non-Gaussian  nature,  however,  coherence  matching  between 
large  gusts  which  are  in  the  development  stage  [5-1]  should  possibly  be  added 
to  the  simulation  technique.  Finally,  the  parametric  inputs  to  the  model 
require  better  definition  which  can  only  be  achieved  with  more  experimental 
effort.  Currently,  it  is  proposed  to  use  the  values  of  turbulent  intensity 
measured  at  elevations  of  several  kilometers  because  these  are  the  better 
defined  values.  Moreover,  comparing  the  thunderstorm  data  measured  by  ground 
based  towers  to  those  backed  out  of  on-flight  data  recorders  shows  that  there 
is  quite  a large  difference  between  10  sec  averaged  data  and  much  shorter 
averaging  period  data.  This  suggests  the  available  thunderstorm  data  be  re- 
analyzed using  a smaller  time  increment. 

There  is  considerable  need  to  measure  turbulence  in  the  vicinity  of  the 
ground,  particularly  in  the  downdraft  center  of  thunderstorms  and  also  to 
measure  the  gust  gradient  which  could  occur  in  these  severe  storms  phenomena. 
The  same  conclusion  holds  true  for  large-scale  frontal  motion  where  very 
limited,  if  any,  wind  velocity  measurements  as  well  as  turbulence  measure- 
ments have  been  made.  The  type  of  measurement  which  appears  most  needed 
would  be  from  an  array  of  towers,  preferably  500  m in  height,  and  distributed 
in  a line  at  least  over  1 km  in  extent.  Such  an  array  of  towers  is  obviously 
very  expensive  and  the  data  reduction  associated  with  the  vast  amounts  of 
data  generated  is,  in  turn,  costly.  A tower  array  is  available  at  NASA 
Marshall  Space  Flight  Center  with  the  limitations  that  the  towers  are  only 
24  m in  height  [5-2]. 
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As  techniques  for  remote  sensing  of  wind  speed  become  available,  typical 
examples  being  laser  Doppler,  acoustic  Doppler  or  Doppler  radar  methods,  more 
probing  of  thunderstorms  and  of  major  synoptic  fronts  will  undoubtedly 
increase. 

Aircraft  are  also  being  used  to  study  thunderstorms  and  are  providing 
meaningful  data  at  great  altitudes.  It  is  not  likely  that  much  probing  of 
thunderstorms  or  severe  frontal  weather  will  be  carried  out  with  aircraft 
below  500  m.  However,  data  at  this  level  is  extremely  critical  in  developing 
simulator  wind  field  models  for  flight/hazard  definition. 

A possible  technique  for  obtaining  the  needed  data  would  be  to  couple 
aircraft  measurements  with  those  of  the  tower  array  currently  in  existence  at 
the  Atmospheric  Sciences  Division,  NASA/Marshall  Space  Flight  Center, 
Huntsville,  Alabama.  A coordinated  program  utilizing  aircraft  passes  over 
the  tower  array  with  the  towers  in  operation  would  map  out  a well  defined 
two-dimensional  wind  field  and  could  also  provide  three-dimensional  wind 
fields  when  the  angle  of  the  approaching  storm  is  oblique  to  the  tower  array. 

Mathematical  models  currently  being  developed  have  promise  for  use  in 
wind  shear  simulation  applications,  however,  they  are,  in  turn,  dependent 
upon  input  parameters  which  must  be  experimentally  measured  before  meaningful 
results  can  be  calculated. 

A statistical  model  for  the  stable  and  neutral  boundary  layer  has  been 
developed  which  allows  the  probability  of  a given  wind  speed  and  stability 
condition  occurring  simultaneously  to  be  estimated.  The  model  also  predicts 
the  frequency  of  the  combined  conditions  occurring  within  a specified  period 
of  time.  This  model  assumes  a Rayleigh  distribution  of  wind  speeds  for  any 
given  site  in  the  mainland  U.S.A.  This  is  a justifiable  assumption  based  on 
a study  of  138  sites  by  Justus,  et  al.  [5-3].  The  probability  distribution 
of  Richardson  numbers  is  not  quite  as  well  validated,  being  based  on  only  two 
sets  of  data.  It  is  believed,  however,  that  this  is  a good  estimate  of  the 
distribution  of  stability  conditions  with  wind  speeds. 

The  statistical  model  for  the  thunderstorm  is  rather  weak  and  for  the 
synoptic  fronts  nonexistent.  The  best  estimate  of  the  probability  of  a given 
magnitude  of  wind  shear  in  a thunderstorm  that  is  currently  available  is 


based  on  the  ensemble  average  of  20  thunderstorm  cases  along  a prescribed 
flight  path  and  the  standard  deviation  fron.  the  averaged  value.  These  com- 
puted standard  deviations  are  added  at  each  point  along  the  flight  path  as  if 
they  were  independent  points.  This  is  a questionable  assumption  and  further 
work  is  needed  to  develop  the  statistical  model  of  a thunderstorm. 

Before  a statistical  model  of  fronts  can  be  developed  much  more  data  is 
needed.  Due  to  the  lack  of  measurements  of  detailed  wind  speed  profiles 
through  large  fronts  and  also,  to  an  extent,  through  thunderstorms,  reliable 
simulation  flight/hazard  wind  speed  models  cannot  be  developed. 

The  conclusion  of  this  report  is  that  mathematical  models  of  wind  shear 
in  thunderstorms,  stable  and  neutral  boundary  layers  and  fronts  for  flight/ 
hazard  simulation  studies  have  been  developed  which,  based  on  state  of  the 
art  knowledge,  allow  realistic  simulation  of  flight  through  hazardous 
atmospheric  wind  shear.  These  models,  however,  represent  a very  good 
beginning  but  they  should  be  continuously  updated  as  information  becomes 
available  and  standardized  to  establish  consistent  criteria  for  pilot 
training,  avionics  development  and  aircraft  design  and  certification. 
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APPENDICES 

The  following  appendices  are  numbered  according  to  the  section  of  the 
main  body  of  the  report  to  which  they  pertain.  For  example.  Appendix  2 
followed  by  an  alphabetical  symbol  contains  information  relating  to  Section  2, 
Thunderstorm  Wind  Shear.  Appendices  3 and  4 are  referenced  in  the  same 
format. 
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APPENDIX  lA 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 

Friction  velocity 
Stability  parameter,  icu*/fL' 

Coriolis  parameter 

Monin-Obukhov  stability  length  scale 

Horizontal  extent  of  thunderstorm  data  set 

von  Karman's  constant  taken  as  0.4 

Length  scale  of  longitudinal  velocity  fluctuations 

Length  scale  of  lateral  velocity  fluctuations 

Length  scale  of  vertical  velocity  fluctuations 

Mean  speed  of  storm  front 

Longitudinal  wind  speed 

Lateral  wind  speed 

Vertical  wind  speed 

Fluctuating  longitudinal  wind  speed 

Fluctuating  lateral  wind  speed 

Fluctuating  vertical  wind  speed 

Surface  roughness 

Dimensional  height 

Dimensionless  height  z = zf/u* 

Rossby  number  u^^/fz^ 

Power  density  spectrum  function 
Cyclic  frequency 

Turbulence  intensity  of  a wind  speed  fluctuation 


Statistical  correlation  between  cx  and  6 wind  speed  components 

Ensemble  average 

Denotes  dimensionless  parameter 


APPENDIX  2A 


TABULATED  THUNDERSTORM  DATA 

Appendix  2A  contains  the  tabulated  data  interpolated  from  the  thunder- 
storm wind  speed  data  presented  and  documented  in  Goff  [2-1].  These  data  are 
the  fundamental  data  sets  upon  which  the  wind  shear  model  for  thunderstorms 
is  based.  The  user  can  compute  wind  speed  profiles  either  by  hand  directly 
from  the  tables  or  by  coupling  the  data  with  the  lookup  computer  program 
provided  in  Appendix  2C. 

The  rows  and  columns  of  the  table  are  numbered  from  the  left-hand  bottom 
corner.  The  numbers  increase  from  left  to  right  in  the  positive  x-direction 
of  the  original  data  and  from  bottom  to  top  in  the  positive  z-direction  of 
the  original  data.  The  wind  speeds  are  given  in  units  of  meters  per  sec, 
m s \ with  being  positive  in  the  direction  of  frontal  motion,  being 
positive  upward  and  being  positive  into  the  plane  of  the  paper.  The 
upper  portion  of  the  table  covers  columns  1 through  21  and  the  lower  table 
covers  columns  21  through  41.  (Note  column  21  is  repeated  for  symmetry  and 
clarity  of  presentation.) 

The  thunderstorm  case  numbers'  designation  for  this  report  are  listed 
at  the  top  of  the  table.  The  letters  in  parentheses  and  the  following  series 
number  correspond  to  the  thunderstorm  designation  given  by  Goff  [2-1]. 

Also  listed  at  the  top  of  the  table  are  the  frontal  speed,  W^,  and  the 
horizontal  length  scales  for  the  given  wind  record.  The  horizontal  extent 
of  each  data  set  varies  because  of  the  data  reduction  procedure.  Hence,  the 
length  of  field,  L,  in  meters  and  the  horizontal  grid  spacing.  Ax,  (m)  are 
specified  on  each  table.  The  vertical  extent  of  each  field  is  taken  as 
500  m with  50  m vertical  grid  spacing.  Figure  2A-1  illustrates  the  identi- 
fication format  used  for  all  tables. 
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Table  2A-l(a)  Wind  Speed  in  Direction  of  Frontal  Motion,  W 

Case  1(A)  Series  No,  0446,  ^ 

ms,  -ength  of  Field  L = 4000  m / Horizontal  Increment  Ax  = 100  meters 
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Case  1(A)  Series  No.  0446, 
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Table  2A-1(c)  Vertical  Wind  Speed, 

Case  1(A)  Series  No.  0446, 


99 


101 


Table  2A-2(c)  Vertical  Wind  Speed, 

Case  2(B)  Series  No.  1314, 

Length  of  Field  L = 3278.8  m , Horizontal  Increment  Ax  “ 81.97  meters 
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Table  2A-3(c)  Vertical  Wind  Speed,  W. 
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Table  2A-4(a)  Wind  Speed  in  Direction  of  Frontal  Motion,  w 

X 

Case  4(D)  Series  No.  1459, 

6 m Length  of  Field  L = 7692.4  m,  Horizontal  Increment  x - 192.31  meters 
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Table  2A-5(a)  Wind  Speed  in  Direction  of  Frontal  Motion, 

Case  5(E)  Series  No.  1924, 

=16.7  m s , Length  of  Field  L = 11428.4  m , Horizontal  Increment  Ax  - 285.71  meters 
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Table  2A-5(b)  Wind  Speed  Perpendicular  to  Frontal  Motion,  W 
Case  5(E)  Series  No.  1924 

16.7  n s . unetl,  of  Flold  L - 11428.4  „orl4„„’tal  Incro.ent  Ax  - 285.71  .oterx 
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Table  2A-6(b)  Wind  Speed  Perpendicular  to  Frontal  Motion,  W 
Case  6(F)  Series  No.  1933,  ^ 

Length  of  Field  L = 7272.8  m , Horizontal  Increment  Ax  = 181.82  meters 
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Table  2A-8(b)  Wind  Speed  Perpendicular  to  Frontal  Motion, 
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Table  2A-9(c)  Vertical  Wind  Speed,  W 

Case  9(1)  Series  No.  1507, 
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Table  2A-10(b)  Wind  Speed  Perpendicular  to  Frontal  Motion,  W^ 

Case  10  (J)  Series  No.  2206, 

13.1  m s , Length  of  Field  L = 8510.8  m , Horizontal  Increment  Ax  = 212.77  meters 
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Table  2A-10(c)  Vertical  Wind  Speed, 

Case  10  (J)  Series  No.  2206, 

Length  of  Field  L = 8510.8  m , Horizontal  Increment  Ax  = 212.77  meters 
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Table  2A-ll(c)  Vertical  Wind  Speed,  W^ 

Case  11 (K)  Series  No.  2118, 

12.4  m s”^.  Length  of  Field  L = 8163.2  m,  Horizontal  Increment  Ax  = 204.08  meters 
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Table  2A-13(a)  Wind  Speed  in  Direction  of  Frontal  Motion,  W 

X 

^ Case  13(M)  Series  No.  0211, 

m s , Lenght  of  Field  L = 63A9.2  m , Horizontal  Increment  Ax  = 158 
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Table  2A-14(a)  Wind  Speed  in  Direction  of  Frontal  Motion,  W 
Case  14(N)  Series  No.  0109, 

s , Length  of  Field  L = 5128.4  m , Horizontal  Increment  Ax 
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Table  2A-14(b)  Wind  Speed  Perpendicular  to  Frontal  Motion, 
Case  14 ^N)  Series  No.  0109, 
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Table  2A-16(b)  Wind  Speed  Perpendicular  to  Frontal  Motion, 

Case  16(P)  Series  No.  0220, 

> s , Length  of  Field  L = 4444.4  m , Horizontal  Increment  Ax  = 111.11  meters 
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Table  2A-16(c)  Vertical  Wind  Speed, 

Case  16(P)  Series  No.  0220, 

Length  of  Field  L = 4444.4  m , Horizontal  Increment  Ax  = 111,11  meters 


Table  2A-17(b)  Wind  Speed  Perpendicular  to  Frontal  Motion,  W 

Case  17 (Q)  Series  No.  1837, 

= 17.5  m s , Length  of  Field  L • 11428.4  m',  Horizontal  Increment  Ax  = 285.71  meters 
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Table  2A-18(a)  Wind  Speed  in  Direction  of  Frontal  Motion,  W 
Case  18(R)  Series  No.  1651,  ^ 

25.0  m s , Length  of  Field  L = 16000.0  m.  Horizontal  Increment  Ax  = AOO.O  meters 
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Table  2A-19(b)  Wind  Speed  Perpendicular  to  Frontal  Motion,  W 

Case  19(S)  Series  No.  1904,  ^ 

Length  of  Field  L = 7142.8  m , Horizontal  Increment  Ax  = 178.57  meters 
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Table  2A-20(a)  Wind  Speed  in  Direction  of  Frontal  Motion,  W^ 

^ Case  20 (T)  Series  No.  0029. 

19.9  n>  s , Length  of  Field  L = 13333.2  m , Horizontal  Increment  Ax  = 333.33  meters 
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Table  2A-20(b)  Wind  Speed  Perpendicular  to  Frontal  Motion 
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Case  20(T)  Series  No.  0029, 

19.9m  s , Length  of  Field  L = 13333.2  m.  Horizontal  Increment  Ax  = 333.33  meters 


APPENDIX  2B 


GRAPHICAL  ILLUSTRATION  OF  THUNDERSTORM  WIND  SPEED  PROFILES 

This  appendix  contains  graphical  illustrations  of  the  wind  speed  pro- 
files encountered  along  flight  paths  having  a 3°  glide  slope  through  each  of 
the  20  thunderstorm  cases.  The  position  of  each  flight  path  relative  to  the  I 

streamlines  fixed  in  the  frame  of  reference  moving  with  the  storm  is  illus-  • 

trated  in  each  of  the  figures,  and  is  chosen  to  terminate  in  the  left-hand  I 

corner  of  each  data  set.  The  flight  paths  are,  therefore,  essentially  arbi- 
trary flight  paths  which  could  be  encountered  during  any  routine  landing 
through  a passing  thunderstorm. 

The  longitudinal,  lateral,  and  vertical  wind  speeds  encountered  along 
the  glide  slope  are  shown.  The  ordinate  in  the  wind  speed  profiles  is 
height,  z,  nondimensionalized  with  the  length,  H = L tan  3°,  where  L is  the 
length  of  the  wind  field.  Each  profile  is  the  wind  as  seen  by  an  airplane 
traveling  along  the  flight  path  drawn  across  the  streamlines  as  shown  in  the 
upper  figures.  Note  that  the  streamlines  plotted  are  relative  to  the  speed 
of  the  front  which,  for  reference  purposes,  has  been  indicated  on  the  hori- 
zontal wind  speed  profile  with  a vertical  dashed  line.  The  wind  speed 
profiles  are  relative  to  the  fixed  earth  frame  of  reference. 

The  purpose  of  these  illustrative  profiles  is  to  provide  an  impression 
of  the  diversity  and  severity  of  wind  shear  that  can  be  encountered  during  a 
routine  approach  when  a thunderstorm  is  in  the  vicinity  of  the  terminal  area. 

The  data  sets  do  not  necessarily  represent  a wind  field  through  the  center  of 
the  storm  nor  the  worst  downdraft  conditions,  but  represent  storms  which  may 
have  passed  either  to  one  side  or  the  other  of  the  tower  from  which  the  data 
were  measured.  Therefore,  the  profiles  do  not  illustrate  necessarily  the  20 
most  severe  flight  paths  which  would  be  encountered  through  the  given  thun- 
derstorms but  rather  a collection  of  wind  speed  profiles  that  represent  an 
averaged  situation. 


157 


158 


Figure  2B-1  Typical  Wind  Profiles  along  a 3°  Glide  Slope:  Case  1(A)  Series  No.  0446, 
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Figure  2B-3  Typical  Wind  Profiles  along 
W = 8.6  m s"^  H - 5714  tan 


161 


Figu 


Series 


163 


Figu 


165 


Figure  2B-8 
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APPENDIX  2C 


COMPUTER  PROGRAM  FOR  CALCULATION  OF  TURBULENT 
THUNDERSTORM  WIND  FIELDS 

This  appendix  describes  a computer  program  used  to  obtain  the  longi- 
tudinal, lateral,  and  vertical  wind  speeds  at  a given  point  (x,z)  and  all  six 
velocity  gradients  at  that  point.  The  computer  program  has  the  option  of 
superimposing  a non-Gaussian  turbulent  field  on  the  thunderstorm  wind  speeds. 
When  the  turbulence  option  is  utilized,  an  input  for  time  is  required.  The 
specified  time  designates  to  the  program  the  period  over  which  the  turbulent 
fluctuations  are  computed. 

The  subroutine  VELO  calls  seven  additional  subroutines  as  schematically 
illustrated  in  Figure  2C-1.  These  seven  subroutines  are  entitled  TURB, 
ROMDON,  FILTER,  GAUSS,  FFT,  WIND,  and  BUILD. 

Program  Description 

The  program  is  written  in  FORTRAN  IV  and  generates  the  wind  field  with 
turbulence  superimposed  for  (x,z,T)  inputs  from  the  calling  flight  simulator 
program. 

Subroutine  VELO 

For  the  given  position  and  time  the  VELO  subroutine  furnish  the  wind 
velocity  [U,V,W  output]  and  velocity  gradient  [UDX,UDZ,VDX,VDZ,WDX,WDZ, 
output]  in  a flow  field  which  is  one  of  20  cases  of  cold  air  outflow  from 
thunderstorms  detected  by  National  Severe  Storms  Laboratory  during  1971 
through  1974  [2-1]. 

Also,  a non-Gaussian  Dryden  spectrum  turbulence  fluctuation  is  added  to 
the  wind  velocity. 

VELO  first  calls  the  subroutine  FILTER  which  sets  up  the  transfer  func- 
tion for  tailoring  the  random  signal.  It  then  calls  subroutine  TURB.  The 
subroutine  TURB  processes  the  random  signal  generated  in  subroutine  ROMDON 

through  the  prescribed  filter  system  and  outputs  fluctuations  w , w and  w 

X y L. 
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having  a non-Gaussian  distribution  with  a Dryden  spectrum  form  (Equation  2-4). 
The  transfer  to  VELO  is  achieved  with  a common  statement.  VELO  then  calls 
subroutine  WIND  to  obtain  the  mean  wind  velocities  at  (x,z)  to  which  the 
corresponding  fluctuations  are  added.  The  wind  speed  and  wind  speed 
gradients  are  then  computed  and  the  data  is  output  in  units  of  m s~^  and  s \ 
respectively. 


Nomenclature 

Subroutine  VEL0(X, Z,T,HAM,NBB,U,V,W,UDX,UDZ,VDX,VDZ,WDX,WDZ) 

X Position  vector  in  longitudinal  direction  (user's  units) 

Z Position  vector  in  vertical  direction  (user's  units) 

T Controls  time  period  of  turbulent  signal  (sec) 

HAM  Ratio  of  unit  (meter/user's  units) 

NBB  Integer  number  for  calling  a new  thunderstorm  case. 

NBB=1  calls  a new  thunderstorm  data  set  into  storage 

U Velocity  in  longitudinal  direction  (m  s”^) 

V Velocity  in  vertical  direction  (m  s”^) 

W Velocity  in  lateral  direction  (m  s ^ ) 

UDX  Longitudinal  velocity  gradient  in  x direction  (s  ^) 

UDZ  Longitudinal  velocity  gradient  in  z direction  (s”^) 

VDX  Vertical  velocity  gradient  in  x direction  (s"^) 

VDZ  Vertical  velocity  gradient  in  z direction  (s'^ ) 

WDX  Lateral  velocity  gradient  in  x direction  (s~^) 

WDZ  Lateral  velocity  gradient  in  z direction  (s  ^) 


[INPUT] 

[INPUT] 

[INPUT] 

[INPUT] 

[INPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 


T is  incremented  with  each  integration  time  step  of  the  calling  program 
and  is  a control  variable  for  the  length  of  the  random,  turbulent  signal 
generated.  When  VELO  is  first  called  in  it  generates  a turbulent  signal  at 
height  z for  a specified  time  period  TEND.  Each  time  VELO  is  called  follow- 
ing the  initial  case,  i.e.,  NBB  = 1,  T is  compared  with  TEND,  when  T ^ TEND  a 
new  turbulent  signal  at  the  current  height  z of  time  period  TEND  is  generated. 
Thus  a quasi-turbulence  variation  with  height  is  achieved. 

NBB  is  set  equal  to  unity  for  the  first  time  VELO  is  called.  It  auto- 
matically changes  NBB  = 2 until  another  thunderstorm  is  desired  at  which  time 

the  calling  program  must  reset  NBB  = 1. 
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Listing  of  Subroutine  VELO 


iUbKuOi'iNi:.  VliLiOlX,Z,r,HA'4,f4bb,U,V«w,UijX,uUi6>VLiA,VUL,^DX,ivL)^,j 
COM^ON/Sl'if/IX 

i,  3, 12y  ) 

ClIMMlI.>J/ST  7/l)X 
CUMMUN/Sl’l  2/L)i  , NNN 
CUMMON/bll 3/GG( 3 , 128  ) 

lK(Nbb-n  101,101,10b 

101  lX  = tj‘3b4y 

CAuL.  liL.l't.«(Z) 

CALL,  rUHb(Z) 
rti^oso. 
iMCC  = 0 

103  rtNOsfti^D  + bl'^NNN 

los  If  ( I’-TtiJO)  109,108,108 

108  CAbL  I'UKbCZ) 

'JCC  = MCC-*-l 

GO  10  103 

109  bH=T/UT 

i^r  = fc:R  + l-NiMi3»NCC 
JR=GG( 1 , NX) 

VP=GG(2 ,NX) 

•»P=GG(3,Nr) 

CALL  WiNU(X,Z,HAM,Nl3b,WX,WZ,v«X,iMXX,wXZ,MZx,MZZ,W¥X,NVZ) 

J = l«X  + UP 
V=8Z+VP 
W = N ¥ + W P 
\ 0 = N 1 + 1 

UL)X=wXXi-(LG(  1 ,i>iO)-ilP)/DX 

\/L)X  = wZX+(GG(2,WU)-VP)/0X 

WOX  = V*¥X  + (GG(2,nU)-NP)/OX 

ODZ=WXZ 

t/UZ=MZZ 

imOZ=m  ¥Z 

Rbi'OKN 

bNO 


Subroutine  TURB 

The  subroutine  TURB  is  called  by  the  subroutine  VELO.  This  sub- 
routine carries  out  the  calculation  for  the  turbulent  fluctuations  based 
on  non-Gaussian  turbulence  model  of  Reeves,  et  al.  [2-3],  A Dryden 
spectrum  is  used.  A schematic  of  the  non-Gaussian  turbulence  simulation 
model  is  given  in  Figure  2C-2. 
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wide  band  gaussian  linear  scale 

white  noise  filters  factors 

sources 

Figure  2C-2  Physical  Interpretation  of  a Single  Component  of 
the  Non-Gauss ian  Turbulence  Model 


TURB  first  calls  the  subroutine  ROMDON  to  obtain  three  independent  dis- 
crete Gaussian  white  noise  random  functions  N , N.  , N_.  These  signals  are 

a D C 

processed  through  a common  statement  by  the  subroutine  FILTER  having  the 
filter  function 

After  mathematical  analysis  the  simulated  turbulent  fluctuation 
Wy(t)  and  w^(t)  are  determined  and  with  a common  statement  returned  to  VELO. 

Listing  of  Subroutine  TURB 
iUbKUO  ri"4c,  rUKb(Z) 

OlMtNSiUN  KK  c i , 12«  1 , R1  ( 3 , 1 ZK  J , AH  t ) , Ml.  i i , i-.tU  3 , 3 , 1 / h ) 

C JHMUN/S'i  1 2/0  r , flWi'l 
CUiiMUN/bl'l  3/GG(  3,  12H) 

4=1 . 

i r = i 

j = i 

CALL  KUML)UiM(RR,HI,NT,Nr2) 

100  DU  110  K=i,l'iT 

XH(K)=H( IT,J,K)*hH( J,K) 

XI CK  >=U. 

110  cuivriNut 

CALiLi  FfT(XK,Xl,.’JT,MT2,2) 

OJ  111  K = l,i'(T 
GH(iT, J,K)=Xk(K) 

111  cjNriNUfc; 

J=J  + 1 

lC(J-i)  100,100,130 
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liU  DU  140  K=1,NT 

Xi^CK  J=GH(l'r , l,K)*GH(ir,2,K) 

GG(IT,K)=K/(l.+K*H)**0.b*XHCK)+l./(l.+K»K)*»0.t54GH(ll,3,f.  ) 
140  COwriNUk. 

I  r=iT-t>i 

IF(ir-iJ  ^9,99,150 
150  CCJNl'lNUt 
i*ir(irt.Cb»2} 

OJ  loo  1=1,3 

/(i<irL(b,l)  (GG(  1 , K ) , K = 1 , All  1 
le'll-2n54,  Ibb,  IbO 
lb4  <«iKlTe:(6, 3 1 
GJ  TU  loO 
Ibb  W><iri:;(b,4} 

IbO 

1 f'QRMAr(lo(lX,fctl.4)) 

2 e' JKMAX  Cl  5X  , ' LUNGITO  I iM  A L ' ) 

3 t'JKMATC  IbX,  • VtRTlCAL' ) 

4 FJKMATCIdX, 'LATKHAL' ) 

KtTURN 

£M0 


Subroutine  ROMDON 

The  subroutine  ROMDON  is  called  by  the  subroutine  TURB.  This  subroutine 

calls  the  subroutine  GAUSS  which  generates  a random  Gaussian  white  noise 

signal.  Three  series  of  normalized  zero  mean  random  function  Ng(t),  Njj(t) 

and  N (t)  are  transferred  to  the  subroutine  FFT  which  fast  Fourier  transforms 
c 

N (t),  N. (t)  and  N (t)  and  returns  the  transformed  signal  to  TURB. 

d D C 

Nomenclature 


Subroutine  ROMDON  (RR.RI ,NT,NT2) 


RR(1,J) 

Real  part  of  Fourier  transform  of  Ng(t) 

[OUTPUT] 

RR(2,J) 

Real  part  of  Fourier  transform  of  N^(t) 

[OUTPUT] 

RR(3,J) 

Real  part  of  Fourier  transform  of  N^(t) 

[OUTPUT] 

RI(1,J) 

Imaginary  part  of  Fourier  transform  of  N^(t) 

[OUTPUT] 

RI(2.J) 

Imaginary  part  of  Fourier  transform  of  N|j(t) 

[OUTPUT] 

R1(3,J) 

Imaginary  part  of  Fourier  transform  of  N^(t) 

[OUTPUT] 

NT 

Integer  number  of  discrete  noise  signal, 
usually  NT  = 128 

[INPUT] 

NT2 

(2)NT2  _ example,  if  NT  = 128,  NT2  = 7) 
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L i s ting  of  Subroutine  ROMDON 


SUrtKOUriwt,  HUMU(IN(rt«,Rl  , NT,N12) 

CLJrtMON/Sr2/lX 

UlMbw.SlUf-J  «H(i, 128), Rl(3,12«),Xh(12bJ,Xill2ri) 

OiJ  1 lU  1 = 1 , i 

SUm=0. 

L>iJ  105  J = l,l'll 

CALL  GAUS5{IX,1. ,0. ,R) 

RKtl , J)=K 

SUM  = £iUK  + KK(  I , J ) 

103  CU^^lXNUt, 

AVc,  = SUM/NI 
SIG=0 . 

Uu  10b  J = l,iMl' 

Rk(1 , J)=KK( 1 , J)-AVt 
S1G  = SIG  + Kf<(  I , J ) **2 
10b  co.^riNUt 

S1G=(S1G/NT)**U.5 
OU  107  J=1,NT 
KK(i , J)=KK(I , J)/SIG 
10/  Go  'll  iNOrJ 

00  108  d = i,i^'r 
X«(J)=K8(I,J) 

XlCJ )=0. 

108  cuwrii'tut: 

CALL  KKTCXH  ,XI  ,1'iT,  Mr2, 1 ) 

UU  loy  J=1,NT 
KKd  , J J=XK(J) 

8I(i,J)=Xi(0} 
loy  coNXif'iUfc 
110  CU'MflNOt 
RETURN 
END 


Subroutine  GAUSS 

The  subroutine  GAUSS  is  called  by  the  subroutine  ROMDON.  Each  time 
GAUSS  is  called  it  generates  a new  random  number.  Random  numbers  generated 
by  GAUSS  are  white  noise  having  a Gaussian  distribution. 

Nomenclature 

Subroutine  GAUSS(IX,S,AM,V) 

IX  Arbitrary  starting  number,  IX  = 65549 

S Adjusting  number,  setting  S = 1.0  gives  standard 
deviation  of  unity 


[INPUT] 

[OUTPUT] 

[INPUT] 
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AM  Adjusting  number,  setting  AM  = 0.0  gives  zero  mean 
V Random  number  output 


[INPUT] 

[OUTPUT] 


Listing  of  Subroutine  GAUSS 


SUBKKUTli'.fc  GAUSS(IX,S,AM,  V) 

HCAL^b  y 

A=0.0 

JU  bO  1=1,12 

iy=lX*bbbJg 

IK(iy)  5,b,b 

b iy=iy+21474«ib47+l 

b y=iy 

y=y*o.4bbobi jo-y 
ix=iy 
bO  A=A+y 

V = (A-b.U)»b  + AiVl 

KtrOHN 

tNO 


Subroutine  Filter 

The  subroutine  FILTER  is  called  by  the  subroutine  VELO.  The  input  list 
contains  the  height  z.  This  subroutine  sets  up  nine  filter  functions  shown 
in  Table  2C-1,  for  generating  the  three  direction  fluctuating  signal,  w^,  Wy 
and  w^. 

These  are  returned  to  TURB  through  a common  statement  where  H(IT,J,K) 
equals  H^(K)  if  J = 1,  H|^(K)  if  J = 2 and  H^(K)  if  J = 3.  The  value  of 
IT  = 1 specifies  longitudinal  fluctuations;  IT  = 2 vertical  and  IT  = 3 
lateral . 

Listing  of  Subroutine  FILTER 

SUbHUUTJivb  KILftRCZ) 

CuMWUM/Sl'4/Ht  3 , 3 , 12h) 

COMMON/Sl  12/DT,NNlV 
N 1 = N N N 

Zl=Z*i.281 

UMtAN=b. 

ALU=(Zl/(0.l77  + 0.832*0.001*Zn**1.2)/J.281 
Al,V  = AI,U 
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TABLE  2C-1 


"cr 

X X 

X 


1 + — ^ s 
'^x 


Lateral  Component 

H,(S)  H|,(S)  H^(s) 


Vertical  Component 


Hjs)  Hj,(s)  H^(s) 
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ALiWsALiJ 

AMU  = 0.10b*H./C0.177  + 0.d32*O.U01*Zn**0.4 

AMV=AMU 

ammsamu 

Ul)U  = ALU/UMc.AN 

VDUsAbV/UMeAN 

WOOsALin/ZUMtAN 

UL)=12d**0.5 

L)t=i.**U.b 

e'C=l ./  (2.*UT) 

DO  no  K = 1,NT 
S = KC*(K-n  /NT 

rt(  1 , 1 ,K)=-*.  ♦AMU^UOd/  ( 1 . 2 . *UuU*i> ) 
H(l,2,K)=l./(l.+2.*UJU*S) 
H(l,i,K)=AMlJ*(2.*UJU)*»0.5/(  1 .+UUU»S) 
H(2,lfKJ=AMU»L)n*VOU**2/Cl.+2.*VuU»S) 

M(2,2,KJ=6/(1 .+2.*VUd*S)**2 

h(2,3,K)=AMV*VuU**0.5*(l.+lJt.*VUU*^J/(l.  + VuU*S)*»2 
riC3,l»'<)=Aiviw*L)t)*i«.rjU**2/(  1 . + 2.*wGU*S) 
tif  3, 2,K)=.S/  (l.  + 2.»-^llJU*S)**2 

HC3, 3,K)=A,'1rt*wtJLI**0.b*(  1 .+Oc;*w(JO»S)/(  1 . i-rtUU  + S)  **2 
110  COiNTiNUt 
^£TljKN 
tNL) 


Subroutine  FFT 

The  subroutine  FFT  provides  a discrete  fast  Fourier  or  inverse  Fourier 
transformation  and  it  is  called  by  the  subroutine  ROMDON  and  the  subroutine 
TURB. 

In  the  subroutine  ROMDON,  FFT  is  utilized  to  transform  N^,  and 
from  the  time  domain  to  the  frequency  domain. 

In  subroutine  TURB,  after  the  frequency  function  has  passed  through  the 
filter  function,  FFT  is  used  to  transform  the  frequency  function  back  to  the 
time  functions  a(t),  b(t)  and  c(t). 

Nomenclature 

Subroutine  FFT(XR,XI ,N,NU,NGG) 

XR  Array  of  the  real  part  of  the  function  to  be  transformed 

XI  Array  of  the  imaginary  part  of  the  function  to  be 

transformed 

N Total  number  of  the  Array  (N);  usually  N=128 
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[INPUT' 

[OUTPUT' 

[input' 

[output; 

[INPUT] 


[INPUT] 


NU  (2)NU=N;  usually  NU=7  which  gives  N=128 

N66  If  N6G=1,  the  subroutine  performs  the  Fourier  Transform 

If  NGG=2,  the  subroutine  performs  the  Inverse  Fourier 
Transform 


Listing  of  Subroutine  FFT 

SUttKOUTltJt;  FFT(XR  ,X1,N,NU,NC.G1 
UiMtNSIUN  XK(N),XI(N) 

If (NGG.£v<.  1 ) GO  TO  321 

UO  i20  1 = 1,1'! 

320  Xl(I)=-Xi(I) 

321  N2=N/2 
WU1=NU-1 
K = 0 

00  100  0 = 1, fvU 

102  DU  101  l = l,i'i2 

P = lttlTR(K/2**NUl  ,NU) 

At<G  = b.  2831H!3*F/f  LQArCM) 

C=CUS( ARG) 

S = Sli'J  ( AHG) 

M=K+1 
K1  i<2=Kl  + N2 

Tk=XK(KlN2 J*C+XI CK1N2)*S 
ri=Xl(KlW2)*C-XK(KlN2)+S 
X8(K1N2)=XH(K1 )-TH 
XltK,li'l2)=XHKl)-TI 
XK(K1 )=Xk(Kl  ) + rK 

xi(Kn=xi(Ki  ) + ri 

101  «.  = K-H 

N=K+N2 

IhCK.Lf.'O  GO  rO  102 
K~0 

i^«J1=n01-1 
100  (^2  = n2/2 

DO  103  K=1,N 

l=iai TKCK-l ,NU)+1 

If (l.Lt.KJ  GO  TO  103 

TK=XK(K) 

ri=Xi(K) 

XK(KJ=XK(1J 

X1(KJ=XI(1) 

XH( 11=TK 
X1(1)=T1 

103  CONl'lNJt 

If  (nGG.L^j.  1 ) GO  TO  121 
DO  124  1 = 1, N 
XK(l)=XHli)/N 
X1(I)=X11I)/!J 


[INPUT] 
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124  CJNlIN'Jh: 

121  RiiTUKN 

fu.'iCriON  ihitkoj.iml)) 
Ji=j 

ltiirH  = 0 

UJ  20U  1 = 1, 1^0 
J2=Jl/2 

iellK=lbiXK*2  + (Jl-2<‘J2) 
^00  J1=J2 
KblUHN 
£i>IU 


Subroutine  WIND 

The  WIND  subroutine  is  called  by  the  subroutine  VELO  for  the  given 
position  (x,z)  from  VELO.  Interpolation  by  an  area-weighting  method  for  the 
velocity  at  the  point  x,z  in  the  flow  field  based  on  the  tabulated  data  for 
the  given  thunderstorm  data  set  placed  in  storage  by  subroutine  BUILD.  The 
area-weighting  interpolation  method  is  illustrated  in  Figure  2C-3. 

The  area-weighting  scheme  calculates  the  velocity  and  velocity  gradient 
of  each  point  by  using  the  four  nearest  neighboring  grid  point  velocities. 


Figure  2C-3  Area-Weighting  Technique 
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The  velocity  at  point  P is  given  by: 


where 


A = A"!  + A2  ■*■  A^  + A^ 

The  same  interpolation  method  is  used  for  the  velocity  gradients. 


Nomenclature 


Subroutine  WIND  (X.Z.HAM, KCK,WX.WZ,WY,WXX,WXZ,WZX,WZZ,WYX,WYZ) 

(X.Z) 

Position  vector  (users  units) 

[INPUT] 

HAM 

Ratio  of  units  (users  units) 

[INPUT] 

KCK 

Integer  number  for  storing  a new  thunderstorm  case. 

If  KCK=1 , the  subroutine  BUILD  is  called  and  a new 
thunderstorm  data  set  is  read  in;  otherwise  KCK=2. 

[INPUT] 

WX 

Velocity  in  longitudinal  direction 

WZ 

Velocity  in  vertical  direction 

[OUTPUT] 

WY 

Velocity  in  lateral  direction 

WXX.WXZ 

Velocity  gradient 

WZX.WZZ 

Velocity  gradient 

[OUTPUT] 

WYX,WYZ  Velocity  gradient 
Listing  of  Subroutine  WIND 

SJbrtUUflrte.  WlNDCX,Z,HAM,KCK,<X,w2,i(ili,«AA,«»XZ,»fZX,W^Z,,wxJC»*<YZ} 
L)I.*ie;Niii(JiJ  axX(2,2),0<Z(2,2)»0ZX(2,2J,uZZe2,i),UltXC2,Z), 
iiiUXZl  2 , 2) 

vJ  J-I^UN/l’lV  A(  41  ,ll,3),OX,PZ 

11  f JK  Am  r ( / / , ' »♦»  X-VAbUt  IN  bOhrtUUTlNt  iNlt'K  lc>  UJlS  f't  FbLtw 
s'KtbliJN  X = ',e;i3.b,//J 

21  i-  JKAmT  (// , ' *»♦  Y-VAL|J£  IN  SUbHuUllNi  li<Ie,K  IS  UUrslUt,  ^ LUki 
S ' KtGlUN  » - ',tli.6,//) 

31  r JKmai(2X,  ' ♦♦♦P'JINI  X = Z = ',bl3.0,'  IS  JufSlUt') 

IFCACk.Nc,.  I ) ^,ll  J'j  100 
CAbo  biJlbO 
NCA=2 
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lUU 

IF  CXF.Lr.o.  ) oU  iu  210 

J J 1 vJl  1 = 2 , •*  1 
Al  = (l-1  )»'JX 

t r ( Xt' . Lr.  , X 1 ) o'J  rj  102 

101  JiJi'l  <0e. 

.’ir<ll’t,(,r>,ll)  XP 

102  1:»=X-1 

OJ  10  J .J  = 2,ll 
2.1  = (.J-1  J*0Z 

IF  ( ZH  . L.C. . 4 1 J llu  I'l)  10  4 
10  1 C J < 1 1 0 F 

mUlTtCb, 21 1 ZH 

104  0^=0-! 

ML,KA=(AF-CiF-l)*UX)/OX 
tlei'A=(ZP-(JF-l)»UZ)/UZ 
x=  1 

105  VAL=(1.0-ALFA)*(1.0-bc;rAj*A(lP,JH,Xl+tle.TA*(1.0-ALFA)*A(IP,JP 
S + l,^)•*•Ab^A♦(1.0-BK^A)♦A(IP+l,JP,^)+AbFA*Ht.'lA»A(lP+l,JP+l,K) 

lK(^-21  lOb,107,10H 

1 0 b « X = V A Li 
10/  <1 1-  VhL 

n = K+  I 
;5J  I'Ll  10b 

loo 

UJ  20b  1=1,2 
OU  20b  J=l,2 

lp'(  IP.tJ.  1 . AniO.  1 . Ed.  1 ) r,J  lU  201 
1 p I i P . c.ij  . +0  . Anu  . I . to  , 2 ) oU  lU  20'2 

L)XX(l,J)=0.:j/oX»(A(lP  + l,J?+J-l,l)-AUP  + i-2,JP4j-l,l)l 
DZXC 1 , J)=U.3/0X*(A( lP+1, JP+J-1 ,2J-AiiP+I-2,JP+0-l,21) 
o/XCl,oJ=0.b/OX»(A(IP4l,JP4j-l,3)-A(lt'+l-2,JP4J-l,i)) 
oJ  fU  201 

201  OXXCl,>J)  = lA(2,JH+J-l,n-A(l,JP  + J-l,ll)/0X 
OZXC  1 ,J1  = (A(2,JP-*'J-1,2)-A(1  , JPfO-1  , 21  ) /UX 
Ol(XCl,Jl  = CAt2,JP+J-l,3)-A(l,JP4j-l,lJl/L)X 
OJ  ru  203 

202  OXXC2,J1=Ca(41,JP+J-1,1)-AC40, OP +0-1,111/ OX 
OZX  C 2 , 0 1 = l A( 41 ,0P  + J-1 , 21 -A(40 ,0P  + J-1 , 2 ) 1 /OX 
01/X12, J1  = IA(41,JP  + . 1-1,31-4140, JP  + J -1,311/ OX 

203  ir  (OP.  to.  1 . AND.  J .KO.  1 1 GO  '1  U 204 
lf'(JP.E.3,10.AND.J.EJ.21  GD  lo  20b 

OXZ(I,01=O.b/OZ*(A(lP+l-l,0P40,ll-A(lP+i-l,JP+J-2,ll} 
UZZ(I,J}=O.b/oZ*(A(IP+l-l,0P+J,21-A(lP+l-l,JP+J-2,21) 
01/Ztl,01=0.b/LiZ*(A(IP+l-l  ,0P40 ,31-AUP4l-l,JP4J-/,311 
v>  J 1 iJ  20b 

204  0<Z(i,ll=lA(lP+l-l, 2, 11-A(1 P+1-1,1, 111/OZ 
OZZ(l,ll=(A(IP+i-l,2,2)-A(lP4l-l,l,Zll/UZ 
0 1Z(1,11  = 1a(1P  + 1-1,2,31-aUP'M-1,1»-i1)/DZ 
u J 1 O 2 0 d 
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20b  L)*<iCI,2J  = CA(lP  + l-l,ll,lJ-AtiH  + l-l,lU,l})/iJZ 
JZZ(l,2)  = {A(lF  + I-l,ll,2J-AUt^^i-l,lU,2J)/uZ 
t>i^(i,2)  = (A(iF+I-l,ll,3)-AClF  + l-l,lU,3J}/i)Z 
CJiU'li^Oe. 

-«AA=(1.0-ALi-A)*(1.0-b£rA)*UXXtl,U+bLlM^(]  . O-Abt-  A ) »UXX  ( 1 
*,2)+Abl-A+U.O-AbrA)*OXX(2,U  + ALfA*ric:j.A»OXX(i:,^) 

^XZ=(1.0-ALhA)*(1.0-Oc:rA)+UXZll,n+clt,lA*(l.0-ALbA)*DXZCl 

♦ ,2)+AL.t'A*(l.U-bfc,fA)*OXZ(2,l)  + ALKM*DilA»OXZl2,2) 
AbX=(l.O-ALKA)  + (l.0-tit:rA)*UZX(l,n+bc:iA*(l.O-ALt-'A)*UZX(l 

♦ , 2 J + ALl‘  A*l  1 .O-BKTA  ) ♦DZX  ( 2 , 1 ) + ALh  A*Btl  A*  JZX  (2,2) 

^ZZ=(  1 .0-AuKA  ) ♦ ( 1 .O-bti'A)  *L»ZZ(  1 , 1 J+btl  A*C  1 .U-ALFA)  *UZZ(  1 
♦.2)-*-AH-A*Cl.U-HbrA)*DZZ(2,l)+ALl‘A*tie,l'A*uZZC2,2) 

X = C 1 .O-ALKA ) ♦( 1 .O-bKf A) ♦biXCl , 1 ) +bbl A*C1 .O-ALf  A) ♦OYXl  1 

♦ ,2)+AL.fA*ll.O-dErA)*DKXC2,l)+AliFA*BbrA*OiXl2,2) 
•«iZ=(1.0-Al,bA)*(1.0-bbrA)*DKZ(l,lJ+BblA»(1.0-ALFAJ»U)(Z(l 

♦ ,2J+ALi»-A*(1.0-htTA)*UKZC2,l)i-Al.tA>t‘btTA*0YZC2,2) 

1,J  ru  21*5 

210  ^-<lfbCb,Jl)  XP,2,F 
213  KciJKN 
bAto 

Subroutine  BUILD 


The  subroutine  WIND  is  called  by  the  subroutine  BUILD  which  reads  indi- 
vidual sets  of  wind  speed  data  into  storage.  The  thunderstorm  cases  are  read 
in  sequential  order  based  on  the  numbering  system  utilized  in  this  report; 
that  is,  Case  1 corresponding  to  Case  A,  Serial  No.  0446  of  Reference  [2-1] 
is  read  in  first  and  remains  in  storage  until  KCK  in  list  statement  of  sub- 
routine WIND  is  assigned  the  value  1.  The  subroutine  WIND  then  calls  the 
next  thunderstorm  case  in  numerical  order. 

In  subroutine  BUILD,  the  grid  system  for  the  wind  field  is  shown  in 
Figure  2C-4  where  IT=41  and  JT=11. 


Figure  2C-4  Grid  System  of  Wind  Field 
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IT 


At  point  (I.J)  the  normal  velocity  W^(I,J)  and  vertical  velocity  W^(I,J)  are 
stored  as 

A(I,J,1)  = W^(I,J);  A(I,J,2)  = W^(I,J) 

A common  statement  is  used  to  transfer  the  data  A(I,J,K)  to  the  subroutine 
WIND. 


Listing  of  Subroutine  BUILD 

SUbKdU  riwf:  BUILD 
jiiHcwsuit  iBim) 

Cj"(MtJN  / rT/A(  41 , 1 1 , i)  , DX  , DZ 
1 f JBMAfCZuX,  'HOKiZ^MTAL  ViLQCl'lI  ' ,/,/J 

i C JKMATd 

J 1- Ji<riAi'(4X,  , 2X,  ZIPS.  1 ) 

4 f JK.'>Ai’C/*/,2uX,  ' Vt-KllCAL  VhLUCll)^  ' ,/,/J 

3 K JHMAK  2f  lu.4,e'b.  1 ) 

o r' Jr<i<lAT(3A,  'uX  =' ,f'B.2,3X, 'UZ  = ' , ^ . X , JX , ' *lt,AN  VX  =',^5.1) 

7 KJHMArtJX,'J/I',2X, 21(13, 2X1) 
b r JKrtAT(/,/,20X,  'LAFthiAL  V LLUC 1 1 i ' , / , / J 

Kc,AJ(b,b)  UA,DZ,IJMe,Ai^ 

«Kii't(b,o)  DX,OZ,UrttAN 
UJ  yO  1=1, 41 
90  lt)(l)  = i 
A = 1 

<»Kll't(o,l) 

100  .Vr<llt,(o,7)  (1B(I),  1 = 1,41,21 

DJ  10b  J = 1 , I 1 

KiAJ(b,21  ( A( I , J , K ) , 1 = 1 , 41  1 
lrC^.^,e,.21  OU  i’ll  10b 
UJ  101  1=1,41 

A(l,J,M=A(l,J,Kl+U'1fe,AiV 

101  CJ.<lll40t, 

10  3 C J n 1 1 N U C, 

u J 102  J=1  , 1 1 
JXJ=12-J 

lOX  AHllfcCb, 31  JK J , ( A( 1 , J<J , K1 , 1 = 1 , 41  , 2 1 
^ = ^•^•l 

lp(K-3)  io3,104,10b 

103  ■vKli'L(b,4l 

oj  ru  100 

104  »«Xlft;(b,91 
sjJ  i'U  100 

10b  KtrUKN 
e.  4u 
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APPENDIX  3A 


TABULATED  STABLE  AND  NEUTRAL  BOUNDARY  LAYER  DATA 

This  appendix  contains  the  tabulated  data  interpolated  from  the 
extensive  measurement  of  wind  speeds  under  varying  conditions  of  stability 
reported  by  Clarke  and  Hess  [3-1]. 

The  tabulated  values  are  the  longitudinal  and  lateral  wind  speeds 

^ A ^ A ^ 

W (z)  = W /u*  and  W (z)  = W /u*  for  neutral  and  stable  boundary  layers  for 
^ ^ y y /\ 

y > 0 where  the  dimensionless  height  z = zf/u*.  The  data  actually  stored  on 
computer  cards  and  used  with  the  computer  lookup  routine  given  in  Appendix  3C 
is 

AW^  = W^(z  = 0.15)  - W^(z) 
and 

AW^  = Wy(z)  - Wy(z  = 0.15) 

These  data  have  been  converted  to  W and  W using  the  relationships 

X y 

W (z)  = aQ  (z  = 0.001)  - Aft  (z)  + ftjz  = 0.001) 

XX  XX 

and 

ft^(z)  = Afty(z)  - AW^(z  = 0.001) 
where 


ft^(z  = 0.001)  = K‘’[p,n(0.001  Ro  + 1)  + 0.01125  y] 

To  extend  the  tables  to  values  of  z < 0.001,  wind  speeds  can  be  computed  from 
ft^(z)  = 1 {£n(Ro  z + 1)  + 4.5  zy/ic};  y ? 0,  0 < z 5 0.001 
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and 


Wy(z)  = 0 

The  data  are  presented  in  a right-hand  coordinate  system  with  positive 
in  the  direction  from  left  to  right  and  positive  into  the  plane  of  the 
paper. 

Tabulated  results  are  given  for  values  of  y ranging  from  0 to  200  in 
increments  of  Ay  = 10  listed  across  the  top  of  the  tables  and  for  values  of 
ranging  from  0.001  to  0.15  listed  in  the  vertical  left-hand  column. 

From  these  tables  the  user  can  construct  wind  speed  profiles  for  given 
values  of  y,  u*  and  f.  Note  that  z = zf/u*  and  Ro  = u*/zf. 
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N > 


Table  3A-l(a)  Longitudinal  Wind  Speed,  W^,  Ro  =10^ 
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76 
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77.2 
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7 J.2 
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, 1 

74.4 

75 

. 0 
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7b,  7 

7b.  7 

75.7 

75.8 

75.8 

O.OoO 

70.  7 

7 1.4 

72. 

, b 

72.0 

7 3 

. 2 

74.1 

71.4 

73.8 

71.9 

74.0 

74.2 

0 . 0 b 0 

OH  . 7 

04,4 

70. 

, 6 

7u.  / 

7 1 , 

. J 

72.1 

71.4 

71. b 

71,4 

71.2 

71.1 

0.040 

bb . 2 

b *» . 2 

66. 

9 

bo  . 4 

07. 

. J 

od  . 1 

be.o 
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67.8 

67.7 

67.7 
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bu  . 0 

bO.O 

6 1 . 

5 

b 1 . 2 

bl  . 

. b 

02.4 

62.0 

02.9 
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bo . 0 
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52. 

J 
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14.3 

34.5 

14.8 

35.1 

U.OOl 

ib.d 

2b.  1 

26. 

4 

29.7 

27. 

.0 

27.2 

27.5 

27.8 

28.1 

28.4 

28.7 

196 


Table  3A-l(b)  Lateral  Wind  Speed,  W^,  Ro  =10^ 

0 10  20  JO  40  bO  60  7U  60  90  100 

O.lbO  -IB.O  -20.4  -72.9  -26,6  -2e./  -J2.U  -Jl,4  -JO,  7 -29.8  -28.9  -28.0 

0.140  -17. H -19. b -22.4  -24.9  -28.0  -il.2  -JO. 6 -29. b -2b. 9 -27,9  -27.0 

O.IJO  -17.1  -19,4  -21.7  -^^.2  -27,2  -3o.4  -29.7  -2a. 9 -2b. 0 -27.0  -2<>.l 

0.120  -17. U -19.1  -21.2  -2J.9  -/o.J  -29.6  -2b. b -2b. 1 -27.1  -26.1  -26.2 

0.110  -Ib.J  -18.4  -20.6  -22. D -'26.3  -2b. 3 -27. o -2b. 7 -26,6  -24, e -23.7 

0,100  -16,6  -17. b -19,9  -21. d -24. J -27.2  -26.2  -25.2  -24,0  -23,0  -22.1 

0.090  -14.6  -17.  1 -19.4  -20.9  -23.2  -26.9  -24.7  -23.6  -22.3  -21. J -20.4 

U.ObO  -13.6  -16.6  -17,4  -lb. 9 -20,9  -23.2  -22.3  -21.4  -2o.4  -19, e -lb, 8 

0,070  -12.1  -13.6  -16.0  -10,6  -Id.o  -2l.l  -20.3  -19,3  -18,0  -17,2  -16.6 

0.060  -10.4  -11.3  -12.4  -13.9  -lo.4  -19.4  -lb. 4 -17.2  -16.2  -l4.3  -13.9 

0.060  -8.7  -9.8  -10.9  -12. u -14.2  -17.0  -16.9  -14.7  -13.3  -12.3  -n,'6 

0.040  -7.0  -b.l  -9.2  -10.2  -12.1  -14,6  -13.3  -12.1  -10. V -9.9  -9.0 

0.030  -b.3  -6.2  -7.2  -b.4  -10.0  -11. b -10. o -9.4  -b,2  -7.2  -6.3 

O,o20  -J.7  -3,9  -4,4  -6.4  -7.1  -9.1  -7.'s»  -b.b  -6.b  -4,9  -4.0 

O.olO  -l.b  -l.b  -2.0  -2.6  -3.0  -4.b  -4.1  -3.4  -2.9  -2.4  -1.8 

C.006  -0,8  -O.b  -0.9  -1.1  -1.0  -2.1  -l.b  -1,6  -1.3  -1.0  -0.0 

O.OUl  0.0  0.0  0.0  0.0  U.O  O.o  0.0  0.0  0.0  O.o  0.0 
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Table  3A-2(a)  Longitudinal  Wind  Speed 


. W^,  Ro  -10^ 
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62.9 

U.OHO 

3V.1 

Jb  • 0 

40.1 

44. to 

4b, 3 

51.1 

5 3.2 

55,2 

b6.y 

56.3 

59.4 
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1 10 
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75.8 

76,1 

76.4 
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0 .090 

09. 1 

7 0.0 

7D.  7 

71.1 

M.7 

72.4 

7 2.4 

72,4 

72.4 

7 2.0 
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Table  3A-2(b)  Lateral  Wind  Speed,  W^,  Ro  =10^ 


0 10  20  30  40  bO  00  70  BO  VO  100 

0.130  -18.0  -20.4  -22.0  -2S.6  -2b. 7 -32.0  -31.4  -30.7  -29. B -28.9  -28.0 
O.ltO  -17.4  -19.8  -22.3  -24.9  -2o.o  -31.2  -30.0  -29.8  -28.9  -27.9  -27.0 
0.130  -17.1  -19.4  -21.7  -24.2  -21.2  -30.4  -29.7  -28. V -28.0  -27.0  -2b. 1 
0.120  -17. 0 -19.1  -21.2  -23.4  -20.3  -29.3  -28.8  -28.1  -27.1  -20.1  -2b. 2 
0.110  -10.3  -18.4  -20. b -22.0  -23.3  -28.3  -27.0  -2o.7  -2b. 0 -24.8  -23.7 
0.100  "13. b -17. H -19.9  -21.8  -24.3  -27.2  -20.2  -2b. 2 -24.0  -23.0  -22.1 
0.090  -14.3  -17.1  -19.4  -20.9  -23.2  -23.9  -24.1  -23.3  -22.3  -21.3  -20.4 
0.080  -13.3  b.b  -17.4  -18.9  -20.9  -23.2  -22.3  -21.4  -20.4  -19.0  -18.8 
O.Olo  -12.1  -13.3  -Ib.O  -10. b -18.0  -21.1  -20.3  -19.3  -18.0  -17.2  -16.3 
0.000  -10.4  -11.3  -12.4  -13.9  -10.4  -19.4  -18.4  -17.2  -lb. 2 -14.3  -13.9 
o.obo  -8.7  -9.8  -10.9  -12.0  -14.2  -17.0  -13.9  -14.7  -13.3  -12.3  -11. b 
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Table  3A-3(a)  Longitudinal  Wind  Speed 


, W^,  Ro  =10^ 
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64 . 3 

u.  luo 

JO,b 

3b. 0 

39.7 

44.8 

48 . b 

bl  .9 

b4.b 

b7.3 

bO.O 

62.1 

64,0 

u.u^u 

iO.b 

3b. 0 

39.7 

44.8 

4d  • 0 

bl,9 

b4.b 

b7.3 

b9,9 

01.8 

63,4 

u . OdU 

30. b 

3b. 0 

39.7 

44. d 

4d  . 0 

bl  ,9 

b4.4 

bb.7 

bd . 8 

60. b 

62.0 

0 . U ^0 

30. b 

3b. 0 

39.7 

44.8 

48.4 

bl  .2 

b3,b 

bb  . 0 

b7.b 

b9. 1 

60.6 

U.UeU 

30.4 

34. H 

39. b 

44.0 

4 7.7 

49.9 

bl  .4 

b4.0 

bO.O 

b7,7 

b9.2 

U .0^0 

:^b.4 

32.4 

3b. b 

41 .0 

4b.  1 

47.8 

49.9 

bl.9 

b4.0 

bb,7 

57.2 

0 . 04U 

^0.4 

3C.2 

34.4 

39.1 

42. b 

4b. 4 

4?.b 

49.4 

bl.2 

b2.b 

53.6 

0.  OiO 

34.4 

2f.2 

32.3 

3o-7 

39.8 

42.3 

44.3 

4b. 9 

47.0 

47.9 

48. b 

U.uiu 

22.3 

2e,2 

30.0 

33.4 

in,  1 

34.4 

ib.3 

3b. b 

3b. 0 

3o.8 

3d.b 

U.OlU 

1 7 ,d 

2 : . 0 

22.3 

24.3 

24.7 

24.4 

24.3 

24,2 

24.2 

2b. 0 

26.4 

0.U03 

14.2 

lb. 4 

lb. 6 

1-7.7 

18.0 

18.0 

18.1 

18.3 

18.4 

18.4 

19,7 

O.UOl 

11. b 

i:  .0 

12.1 

12.4 

12.0 

12.V 

13.2 

1 3.b 

13.8 

14.0 

14.3 

lOU 

no 

120 

130 

140 

ibo 

160 

170 

180 

l9o 

200 

U.  1 bO 

04  . J 

65.8 

67. 1 

68 .0 

oe . 9 

69.7 

70,0 

70.  J 

70.6 

70 . 9 

71.1 

0.140 

64 , 3 

65.8 

67.1 

67.9 

68.7 

o9.b 

09 , b 

69.8 

70.0 

70.3 

70.6 

0, 1 JO 

64.3 

o5.7 

06.9 

67.6 

6d  . 3 

69.0 

bV.O 

09.2 

o9.b 

69.7 

70.0 

0. 1 20 

64.3 

65.6 

66.6 

67.1 

0 7 . 6 

68.2 

66 . 4 

ob  • 0 

68.9 

69.1 

69.4 

0.110 

04.3 

65.4 

66.2 

66.  b 

67.0 

e7.8 

o7.9 

86.0 

68.2 

88.4 

68 . 7 

0. 100 

64.0 

64 . 9 

65.6 

6b. 9 

66  . b 

67.2 

o7.3 

87.4 

67 , 4 

67.6 

67.9 

U , O90 

o3,4 

64.3 

6b. 0 

bb.3 

6b. 9 

66.7 

00.7 

06.7 

66.6 

66.8 

67.1 

0 . odO 

e2 . 0 

83, 1 

63.9 

64.3 

64.9 

6b . b 

Ob  . 0 

8b  • 0 

e5  • 6 

65.7 

65.7 

O.O  fsi 

60.6 

61.7 

62.6 

62.9 

o 3 . b 

04.2 

04  • 2 

64 . 2 

6 4.2 

64.2 

64.3 

0.000 

b4.2 

80.2 

61 . 0 

ol  . 1 

61.7 

0/  . o 

64.4 

02.3 

62.4 

62.5 

62.7 

O.ObO 

b7.2 

58.3 

59.1 

b9,2 

b9  . d 

00.6 

00 , 4 

60.1 

59.9 

59.7 

59,6 

0.040 

53.6 

54.7 

55.4 

5b. 3 

bb.  8 

bO.O 

bo  • b 

b6,4 

56,3 

bb,  2 

56.2 

O.OiO 

4b. b 

49.4 

50.0 

49.7 

bU.  1 

bu . 9 

bl  . 1 

bl.3 

51.6 

bl  .8 

52.0 

0.020 

38  . b 

39, b 

4 0.8 

41.4 

42.1 

4i  .0 

43.2 

43.5 

43.8 

44.0 

44.3 

O.OlU 

26.4 

27.4 

28,3 

2d. 9 

29,6 

30.2 

30. b 

30  • d 

31.1 

31.4 

31.0 

0 .005 

19.  / 

20.3 

20.8 

21.3 

21.7 

24.2 

44. b 

22.7 

23.0 

23.3 

23.6 

U.  OUl 

14.3 

14.0 

14.9 

lb. 2 

Ib.b 

lb. 7 

lo.u 

16.3 

16.6 

16.9 

17.1 

200 


Table  3A-3(b)  Lateral  Wind  Speed 


0 

10 

2D  30 

0.150 

-1  d.o 

-20.4 

-22.9  -25,6 

0.140 

-17.4 

-19. B 

-22.3  -24.9 

0. 130 

-17.1 

-19.4 

-21,7  -24.2 

0. 1 20 

-17.0 

-19.1 

-21.2  -23.4 

o.  1 10 

- 1 6 , 3 

-IB. 4 

-20.5  -22.6 

0.100 

-15.5 

-17. H 

-19,9  -21, b 

0.090 

-14,5 

-17.) 

-19.4  -20.9 

0 . 080 

-13.5 

-15.5 

-1 7.4  -lB.9 

0 . 0 7 0 

-12.1 

-13.5 

-15.0  -lo.b 

0 . Oou 

-10.4 

-11.3 

-12.4  - 1 3.9 

0.050 

-b.7 

-9.d 

-10.9  -12.0 

0.U40 

-7.0 

-b . 1 

-9.2  -10,2 

0.030 

-5.3 

-b.2 

-7.2  -B,4 

0.020 

-3.7 

-3.9 

-4.4  -5.4 

0.010 

-1  .B 

-1 .6 

-2.0  -2.5 

0 . 005 

-O.B 

-o.a 

-0.9  -1  . 1 

0 • 00  i 

0,0 

0.0 

0.0  u.o 

100 

1 10 

120  130 

0.150 

-20.0 

- 2 6 . B 

-25, « -25.2 

0.140 

-27.0 

-25.  B 

-24.7  -24.2 

O.liO 

-26.1 

- 2 4 . d 

-23. B -23,3 

0.1  20 

-25.2 

-23,9 

-22.9  -22.4 

0 . 110 

-2  3,7 

-22,5 

-21.5  -21.0 

0 . 1 00 

-22.1 

-21.0 

-20.0  -l9.b 

0 . 090 

-20.4 

-19.4 

-1  b . 5 -1 B . 1 

o • 0 d 0 

- 1 d . d 

-17. d 

-17.1  -16.7 

0.0/0 

-16.5 

-15. d 

-15.3  -15.0 

0 • 060 

•13.9 

-13.5 

-13,2  -13.0 

0 .050 

-11.3 

-11.0 

-10.7  - 10. B 

0.040 

-9.0 

-b.  4 

-8.0  -B . 1 

0 .0  30 

•6,3 

-5.7 

-5.2  -5.0 

U .020 

-4.0 

-3.  1 

-2.7  -3.1 

0. 01  0 

- 1 . e 

-1  . 3 

-1.0  -1.4 

0.U05 

-O.B 

-O.b 

-0,5  -0,6 

0.001 

U.O 

■o.o 

0.0  U.O 

, W , Ro  =10^ 

y 


40 

-i«.7 
-2b. 0 
-27,2 
-2o,  3 
-25. 3 
-24.3 

- 2 0 . y 

- 1 B . O 
-16.4 
-14.2 
-J  2, 1 
-10. 0 

- / . 1 
-3.6 
-1,0 
0.0 

1 40 

-24.2 
-23.3 
-22.4 
-21  ,o 
-20,0 
-Ito.b 

-1  >.l 

- 1 5 . B 
-14.2 
-12.5 
-10. b 

-7.» 
-5.0 
-3.  1 
-1  .*• 
- 0 , o 
0.0 


50 

-32.0 
-31,2 
-30.4 
-2^.5 
-2b. 3 
-27.2 
-25.9 
-23.2 
-21  . I 
-IS. 4 
-17.0 

-14.5 
-11. « 
-S.l 
-4.b 
-2.1 
0.0 

150 

-23,0 
-22.2 
-21.4 
-20.5 
-Ib.b 
-17,3 
-15.  V 
-1*4. b 
-13.3 
-U  . S 
-9.S 
-7  , 7 
-5.3 
-3.0 
-1.2 
-0.5 
0.0 


bO 

-31.4 

-30 . o 
-29.7 
-2b.B 
-27.6 

- 2 b , 2 
-24,  1 
-22.3 
-20.3 

-la, 4 

-15. S 
-13.3 

-10. b 

-7.9 

-4.1 
-1  . d 
0.0 

loO 

-21.6 
-21 , 1 
-20.2 
-19.3 
-17. s 

- 1 b . o 
-15.3 

- 1 4 . U 
-12.7 
-11.3 

-9.7 
-7.B 
-5 . b 
-3.5 
-1.6 
-0.7 
0 . 0 


70 

-30.7 
-29. b 
-2d. 9 
-2H.  1 
-2b. 7 
-25.2 
-23.5 
-21.4 
-19.3 
-17.2 
-14.7 
-12.1 
-9.4 
-b , b 

-3.4 

-1.5 

u.o 

170 

-20. S 
-19, b 
-19.0 
-1  b . 0 
-16.9 
-15. b 
- 1 *4 . 7 

-13.3 

-12.0 

-10.7 

-9.3 

-7.6 

-5  . B 

-3,B 
-1  .B 
-0,  d 

U.O 


BO 

-29, b 
-2B  . 9 
-2b. 0 
-27.1 
-25.6 
-24.0 
-22.3 
-20,4 
-IB.O 
-15.2 
-13.3 
-10.9 
-d.2 
-5.b 
-2.9 
-1.3 
0.0 

1 bO 

-19.1 
-lb. 4 
-17  .b 

- 1 o . b 
-15.6 

- I 4 , b 
-13. b 
-12.5 
-11.3 
-10.2 

-b  . b 
-7.0 
-5.5 
-3.b 
-1.9 
-u.b 
0.0 


90 

-2b. 9 
-27.9 
-27.0 
-26.1 
-24.6 
-23.0 
-21.3 

-1  9.6 

-17.2 

-14.3 

-12.3 

-9.9 

-7.2 

-4.9 

-2.4 

-1.0 

0.0 

190 

-17.2 
-16,5 
-15.7 
-14. b 
-13. « 
-13.0 
-12.3 
-11.2 
-10.1 
•b . 9 

-7.3 

-5,7 

-4.2 
-2.9 
-1 .4 
-O.b 

0.0 


1 00 

-26.0 
-2  7.0 
-2e.l 
-25.2 
-23.7 
-22.1 
-20.4 
-Id.b 

-16,5 

-13.9 

-11.5 

-9,0 

-6.3 

-4.0 

-1,6 

-o.d 

O.D 

200 

-15.0 
-14.3 
-13,6 
-12. b 
•1  1 • B 
-11.1 
-10.5 
-9,  b 
-B.7 
-7,2 

-5.6 

-4,0 

-2.4 

-1.4 

-0.6 

-0.3 

0.0 
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Table  3A 

-4(a)  Longitudinal 

Wind 

Speed,  Vi 

lx.  Ro 

1 =10 

T 

0 

10 

20 

30 

40 

50 

00 

70 

bO 

90 

100 

u.  ISO 

^4.b 

29.2 

3 3.9 

39.0 

42.9 

4o.2 

4 b . b 

51 . 5 

54.2 

bo. 5 

56.6 

U.  MU 

;^4  • 6 

29.2 

33.9 

39.0 

42.9 

4b.  2 

4b. b 

51.5 

54.2 

50,3 

3b. 0 

U.  1 JO 

^4.11 

29.2 

33.9 

39.0 

42.9 

40.2 

4b. 6 

51 .5 

54.2 

56.3 

58.0 

u.  wo 

V4.b 

2 9.2 

33.9 

39.0 

4 2.9 

4 0.2 

4d  . 0 

51  .5 

54.2 

5o.5 

3b.6 

U.  1 10 

24 . b 

29.2 

33.9 

39.0 

4 2.9 

4o.  2 

4b. b 

31.5 

54,2 

30.3 

58.6 

0.100 

24  . b 

29,2 

33.9 

39.0 

42.9 

40.2 

4b. 8 

51.5 

5 4.2 

50.4 

58.2 

O.U'^0 

24. b 

29.2 

33.9 

39.0 

4 2,9 

4o,2 

4b. b 

51.5 

54.2 

5o.  1 

57.0 

U . OtfU 

24. b 

29.2 

33.9 

39.0 

42.9 

40.2 

4b . 0 

51.0 

53.0 

54.7 

50.2 

u.u/u 

24  . b 

29.2 

33.9 

39.0 

42.0 

43.5 

47.7 

49  . b 

51.7 

53.3 

54. b 

0 • OoO 

24.0 

29.1 

33. e 

3d  . b 

«2  • 0 

44 . 1 

40.2 

4b  . 2 

50.3 

32.0 

53.4 

u.udO 

22.0 

20.7 

31.0 

35.9 

49. 3 

42.1 

44. 1 

40.2 

4b. 3 

30.0 

51.4 

U . 0<*U 

20.0 

24.5 

2b. b 

33.4 

3o . b 

39.0 

41.7 

43.7 

43.4 

46.b 

47.9 

O.OJO 

1 b . 0 

22.5 

20.5 

30.9 

34.1 

3o.  0 

3b. 5 

40.2 

41.3 

42.1 

42.8 

U.U^O 

1 o . 0 

20.5 

2 4.2 

27.7 

29.0 

29.2 

29,3 

29.  b 

29,9 

31,0 

32.7 

O.OlO 

1 1 . b 

14.3 

tb,5 

lb. 5 

19.0 

1 b . o 

Ib.o 

lb. 5 

lb. 4 

19.3 

20.0 

0 • 003 

b.4 

9.7 

10.9 

11.9 

12.4 

12.3 

12.4 

12.5 

12,6 

13,2 

13.9 

0.001 

S.b 

6.0 

6.3 

0.0 

b.  9 

7.2 

7.4 

7.7 

b.O 

b,  3 

b.O 

100 

1 lU 

1 20 

130 

140 

150 

100 

170 

1 bO 

190 

2Ut) 

0. 1 so 

se.o 

bO.  0 

61  . 3 

02.2 

03.1 

04.0 

o4,3 

04 . 5 

04  . b 

05.1 

65.4 

0.140 

bb . b 

oO.O 

61.3 

o2. 1 

02.9 

04.7 

o3 « b 

04.0 

04 . 3 

64.0 

64.8 

0 . 1 iO 

bb . b 

oO.O 

61,2 

ol . b 

02.3 

0 4.2 

o3 , 3 

63.4 

03.7 

04 . 0 

64.3 

0.  wo 

bb . 0 

59.9 

o0.9 

01 . 3 

*>1.9 

o2,3 

o2 . 0 

02. b 

03.1 

o3.4 

63.7 

0. 1 10 

bb.b 

59.6 

60.4 

bO.7 

bl  .3 

02. 0 

02.1 

02,2 

02.4 

02.0 

62.9 

U . 1 uu 

be . 2 

59.2 

59.9 

oO . 1 

ou . / 

01.5 

01 .3 

0 1 . 0 

bl.7 

01.9 

62.2 

0 . 0^0 

bl  .o 

5b. 5 

39.2 

39.0 

bo.  2 

00.9 

oO . 9 

o0.9 

0U.9 

01.1 

61.3 

0 . OSO 

30.2 

b7.3 

5e.2 

3 0.0 

59. 2 

59. b 

39.9 

59.9 

59,9 

59.9 

60. 0 

0. 0 70 

34  . b 

50. 0 

5b . 6 

57.2 

57.7 

5b. 4 

3b.4 

5b. 4 

5b. 4 

5b. 3 

5b. 5 

0 . voO 

b3,4 

54.3 

55.2 

35.9 

5b . 0 

5o.9 

30.7 

5o , 5 

50 . 0 

3o  • b 

56.9 

U . 030 

31 . 4 

52. b 

53.4 

53.5 

54 . 0 

34  . b 

34.0 

54 . 4 

54 . 1 

54.0 

53.9 

U .040 

47,0 

49.0 

49.6 

4 9.0 

bo.  1 

30 . a 

30.  b 

50.7 

50.5 

50.4 

50.4 

0 . u Jo 

42. b 

4 J.  7 

44.2 

44.9 

44.4 

43.1 

45.4 

43.0 

43.8 

4o.  0 

40,3 

0.0^0 

J2, 7 

J4.0 

35.1 

33.0 

30.4 

37.2 

37.3 

37.0 

3b. 0 

3b. 3 

38.5 

O.OlO 

20.0 

21.7 

22.5 

23.2 

23. b 

24.3 

24.  b 

25.0 

25.3 

23,0 

25.9 

0.003 

1 J.4 

M.3 

15.1 

15.5 

1 o . u 

lo.4 

10.7 

17.0 

17.3 

17.5 

17.8 

0.001 

b . 0 

b.  H 

9.1 

9.4 

9.7 

lu.O 

10.3 

10.5 

10. b 

11.1 

n .4 
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A 

Table  3A-4(b)  Lateral  Wind  Speed,  W^, 


Ro  =10^ 


0 )U  20 


0. 1 90 

-1«,0  -20.4  -22.9 

0. 1 40 

-17.4  -19.0  -22.3 

0.130 

-17.1  -19.4  -21 . 7 

0.120 

-17. u -19.1  -21.2 

0.110 

-lb. 3 -lb, 4 -20.6 

0.  lUO 

-16.6  -17.0  -19.9 

0.090 

-14.6  -17,1  -19.4 

O.  000 

-13.6  -15.6  -17.4 

0.0/0 

-12.1  -13.6  -16.0 

0 . OOO 

-10.4  -11.3  -12.4 

0.060 

-0.7  -9.0  -10.9 

0.O40 

-7,0  -0,1  -9.2 

o.oJO 

-6.3  -6.2  -7.2 

0 .020 

-3.7  -3.9  -4.4 

0.  OlO 

-1.0  -l.b  -2.0 

U . 009 

-O.b  -O.b  -0.9 

0 . 001 

0.0  0.0  0.0 

100  no  120 

0.190 

-28.0  -26.0  -26. b 

0 « 1 40 

-27.0  -26.0  -24.7 

0. 1 30 

-2o.l  -24, e -23.0 

0. 120 

-26.2  -23.9  -22.9 

0.  no 

-23.7  -22.6  -21.6 

0.  loo 

-22.1  -21.0  -20.0 

0.090 

-^U.4  -19,4  -10.6 

0.000 

-lo.H  -17,0  -17.1 

0 . 0 lo 

-lo.9  -16.0  -15.3 

0 . OoO 

-1 3.9  -13.6  -13.2 

0 .090 

-11.6  -11.0  -10.7 

0.040 

-9,0  -0,4  -8,0 

0.O30 

-b.3  -6.7  -6.2 

0.O20 

-4.0  -3,1  -2.7 

0.010 

-1.0  -1.3  -1.0 

0 . 009 

-0.0  -0.6  -0.6 

0.001 

0.0  O.U  0.0 

30 

-2b. b 
-24. m 
-24.2 
-23.4 
-22.6 
-2  1.0 
-20. b 

-Ib.b 
-13.9 
-12. u 
-1U.2 
-b.4 
-b.4 
-/.b 
-1  .1 
0.0 


40 

-2a.  ■) 
-20.0 
-27,2 
-26.3 
-26.3 
-24.3 

- 2 J . 2 
-20.9 

- 1 6 . b 
-16.4 
-14.2 
-J2.1 
-10.0 

-7  . I 

-3.0 

-1 .0 

u.o 


60 

-32.0 
-31  .2 
-3o.4 
-29.6 
-20.3 
-27.2 

-26.9 
-23.2 
-21.1 
-19.4 
-17.0 
-14.6 
-1 1 . b 

-9.1 

-4.0 

-2.1 

0.0 


bO 

-31.4 
- 3 0 . o 
-29.7 
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Table  3A-5(a)  Longitudinal  Wind  Speed 
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Table  3A-5(b)  Lateral  Wind  Speed,  W^,  Ro  =10 
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APPENDIX  3B 


GRAPHICAL  ILLUSTRATION  OF  STABLE  AND  NEUTRAL  BOUNDARY 
LAYER  WIND  SPEED  PROFILES 

To  provide  an  illustration  of  the  velocity  profiles  which  can  be 
encountered  along  a flight  path  during  approach  through  neutral  and  stable 
boundary  layers,  seven  velocity  profiles  have  been  computed  using  the  model 
developed  in  this  report.  Both  longitudinal  and  lateral  wind  speeds  are 
shown.  The  conditions  for  which  the  velocity  profiles  have  been  computed  are 
shown  in  Figure  3B-1  along  with  wind  speed  profiles  calculated  from  the 
standard  log-linear  equation  for  the  stable  boundary  layer,  see  References 
[3-2  and  3-5].  This  enables  a comparison  of  the  simpler  theoretical  models 
with  the  more  elaborate  simulation  models  developed  in  this  report  to  be 
made.  One  thing  to  note  in  comparing  the  velocity  profiles  is  that  the  log- 
linear,  theoretical  models  give  only  the  longitudinal  component  of  wind 
speed,  whereas,  the  table  lookup  model  provides  both  the  longitudinal  and  the 
lateral  components.  The  lateral  wind  speed  component  is  quite  significant  in 
many  of  the  cases  illustrated  and  represents  an  estimate  of  the  directional 
shear  encountered  during  approach  or  landing.  This  directional  shear  cannot 
be  obtained  from  the  simpler  theoretical  equations. 

Figures  3B-2  and  3B-5  through  3B-7  are  computed  for  a u*  value  of 
0.1  m s”^.  This  value  of  friction  velocity  when  used  to  compute  dimensional 
height  from  the  nondimensional  form  z gives  a value  of  z = 173  m at  z = 0.15 
which  is  the  maximum  vertical  extension  of  the  tabulated  data.  Consequently, 
the  aforementioned  figures  do  not  extend  to  the  500  m level  as  do  Figures 
3B-1,  3B-3  and  3B-4. 

Additional  visualization  of  the  velocity  profiles  which  can  be 
encountered  in  stable  and  neutral  boundary  layers  with  turbulence  super- 
imposed are  given  in  the  text.  Figures  3-8  through  3-13,  Section  3.4. 
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Figure  3B-3  Stable  Boundary  Layer 
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Figure  3B-6  Stable  Boundary  Layer 
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Figure  33-8  Stable  Boundary  Layer 


APPENDIX  3C 


COMPUTER  PROGRAM 

The  subroutine  STB  is  a FORTRAN  computer  program  for  calculating  the 

wind  speeds  W and  W and  the  wind  speed  gradients  9W  /9z  and  9W  /9z.  Inputs 
X y X y 

to  the  program  are  height,  z,  friction  velocity,  u*,  Coriolis  parameter,  f, 
stability  parameter,  p,  and  surface  roughness,  z^.  These  parameters  must  be 
introduced  in  units  of  meters  and  meters  per  second,  respectively.  The 
returned  velocities  and  gradients  have  units  of  meters  per  second  and  inverse 
seconds,  respectively.  Wind  speeds  and  gradients  at  a new  value  of  height, 
z,  are  obtained  by  simply  assigning  a new  value  to  z in  the  calling  program. 

A new  condition  of  stability  is  also  achieved  by  assigning  y a new  value  in 
the  calling  program. 

The  user  has  the  option  of  adding  turbulent  fluctuations  to  the  wind 
speed  by  assigning  the  control  integer  NKK  a value  greater  than  two. 

Subroutine  STB 

The  subroutine  STB  for  a given  height  z and  with  the  parameter,  y, 

specified  provides  the  wind  velocity  W and  W . The  user  has  the  option  of 

X y 

superimposing  turbulence  fluctuations  on  the  mean  wind  speed  if  desired.  The 

turbulence  simulation  routine  uses  the  spectra  developed  in  Section  3.5  and 

the  z-transformation  technique  for  generating  the  fluctuations  w , w and  w . 

X y z 

Subroutine  STB  first  calls  the  subroutine  INIT,  which  reads  the  input 
data  into  storage.  The  control  variable  NST  is  then  set  equal  to  two  and 
further  use  of  STB  does  not  call  the  subroutine  INIT.  STB  then  calls  the 
subroutine  WINDF  which  interpolates  wind  speeds  in  the  z direction  for  a 
prescribed  value  of  y. 

At  the  user's  option,  NKK  (0<NKK)  can  be  set  greater  than  two  and  the 
subroutine  STB  will  call  the  subroutine  DRYDEN  which  processes  the  random 
signal  from  the  subroutine  GAUSS  and  generates  turbulent  fluctuations  having 
the  modified  Dryden  spectra  described  in  Section  3.4.  These  fluctuations 
are  added  to  the  turbulent  wind  field  output. 
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Nomenclature 


Subroutine  STB(Z,AMU,USTR,F,ZO,NST,WX,WY,DWX,DWY,DT,IX,NKK) 

Z Height  (m) 

AML)  Stability  parameter  (-333.33  < y < 216.67) 

USTR  Surface  friction  velocity  u*  (ms 
F Coriolis  parameter,  f (s~^) 

ZO  Surface  roughness,  (m) 

NST  Control  variable  NST  = 1 initiates  the  storage  of  the 
tabulated  wind  speeds 

WX  Wind  speed  in  the  x direction  (m  s"^ ) 

WY  Wind  speed  in  the  y direction  (m  s'^) 

DWX  Gradient  of  W^  in  the  z direction  (s"^) 

DWY  Gradient  of  W^  in  the  z direction  (s"  ) 

DT  Time  increment  At  of  the  turbulent  fluctuation  (s) 

IX  Initiating  integer  for  random  signal  generator,  value 

arbitrary 

NKK  Integer  for  total  terms  of  z transformation 
(if  NKK  s 2 turbulence  is  not  added 
if  NKK  > 2 turbulence  is  superimposed  on  wind 
speed  outputs) 


[INPUT] 

[INPUT] 

[INPUT] 

[INPUT] 

[INPUT] 

[INPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[INPUT] 

[INPUT] 

[INPUT] 


Listing  of  Subroutine  STB 

S JoKUUrii'it,  o'l  -iC  , A viu  , USIX  , t , Z,0  , , >vX  , Ki  , A , Ow  i , Di  , ix  , NKK) 

3wXK(i41,rtl[H(J4),«AMj,Ab*'A,i.lRC.e,Vkt.K,UJi:,K 
CiJ«i-lUN/il/ Ah 

Jc..  1 j liu  TH  100 
i X = obb4'4 
h A b b 1 1'l  L r 
KJofK  = USlK/  IK*Z,0) 

u Si's  2 

iOo  cji’^iiNjfc 
Z = Z H Z 0 

CAbb  /iili<bFlAMO,Z,.VX,rt)(,OK»X,OrtY,USlK,F,KUSiK,lZ) 

IrL^KK-Z)  141,141,151 
141  JJ=0. 

U V S 0 , 

^i)  fo  1 0 1 

151  CAbb  OKYutil  (Z  , Hi'^U  , JSi4  , r , AX  , UO  , L/ V , , Oi  , 1 X , NKK) 

Ibl  aX  = /<X  + OU 
i *0\l 
KKillrtW 
b (iJ 
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Subroutine  DRYDEN 

The  subroutine  DRYDEN  first  calls  the  subroutine  GAUSS  which  generates  a 
Gaussian  white  noise  random  signal.  The  random  signal  is  then  passed  through 
a filter  which  generates  the  fluctuating  output  having  a spectrum  function  of 
the  Dryden  form  which  fits  the  data  of  Kaimal  [3-3].  The  subroutine  DRYDEN 
then  returns  the  fluctuating  wind  components,  DU,  DV  and  DW  to  the  subroutine 
STB. 

The  subroutine  DRYDEN  utilizes  the  stability  parameter,  y,  the  friction 

velocity,  u*,  the  Coriolis  parameter,  f,  and  the  wind  speed  at  height,  z. 

2 

The  scale  frequencies,  Oq,  and  turbulence  intensities,  a , are  computed 
internally  in  the  subroutine. 


Nomenclature 

Subroutine  DRYDEN(Z,SMU,USTR,F,V,DU,DV,DW,DT, IX,NKK) 

Z Height  (m)  [INPUT] 

SMU  Stability  parameter,  y [INPUT] 

USTR  Friction  velocity,  u*  (m  s"^)  [INPUT] 

F Coriolis  parameter,  f (s'^)  [INPUT] 

DU  Wind  fluctuation  in  longitudinal  direction  (m  s"^)  [INPUT] 

DV  Wind  fluctuation  in  lateral  direction  (m  s”^)  [OUTPUT] 

DW  Wind  fluctuation  in  vertical  direction  (m  s"^)  [OUTPUT] 

DT  Time  increment  of  the  turbulent  fluctuations  (s)  [INPUT] 

IX  Starting  integer  for  generating  random  signal  [SPACE] 

NKK  Total  number  of  terms  for  z transform  [INPUT] 


Listing  of  Subroutine  DRYDEN 

SMhKjul  l <e.  ,v,L/y,OY,l)w,  r,iX,NKK) 

1.)  1 c,  .1 S IJ  (I  A I 1 u ) , \ U u J , eo  ( i j , IsloC  J J , MK  c .i  ) , AC  C 3 1 
C l-i'lJ.-i/Sl  / aL 
Z,AI.  = Z/ AL 

kI  = ZAl./(1.-*-4.5  + ZmL/) 

bl  f u(  I )=0.  JlH 
(;ii  r .1  Hi 
bz  u(  I ) =u.  b + y j 
bj  b<4,b4,»ib 

HA  fO(2)=U.026b 
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( j ■ I I J H K 

no  1 r'  ( “^  I ” 0 . U -i  4 J J c / , o / , » o 
o7  (•'0(l)=iJ.O^r)^ 
ijii  I':)  8 4 

88  fi)(l)=^,0*8i 

8':»  l»-(ZA|j-l.)  4^,40,48 
48  lr(/.  A Ij - ]./>/>)  4/, 44, 44 

UM  r J 4 8 

4‘)  srfic  n = i . J-'j.  1 j^iAL  + ^o.s 

(in  I'J  48 

4/  Sic;(  n=&.44-4.3i!*4Ab 

48  si(;(,'j=4iCi(j)/(u.i/;+(j.uOz:7/*4j*»o.4 
oloCl)=Si'vi(8)/(.0.3bi  + U.OOli4»-Z,J^*0.8 
ij  I 4 4 I^(<=  1 , J 

aC(  |<^)=2i.8D♦v♦^■u(^^J/^ 

4 4 C . I I I I '4  J t 

H8=  1 

i ()(J  (M  I lO  1 = 1 , 

(*  #n  I ij  »,  A .J  8 S C 7 A , 1 * , , rt  } 

A ( L ) =8 

I 3 (J  C I 1 1'  1 '1  J r, 

Y(  U=0. 

(i  = ".  A-^C-ACt  ^rtri)4r) 

14  = ( i ) K 1 ( 4 N ( .J  8 8 ) ) ) * t 1 . “ C J / r.  C C ' ( rt  t>  ) 

n I 11]  I = 2 , l''A8 

< ( J ) =C* Y ( i-1  J 4U*AC 1-1  J 

111 

If  1.  ■I88-2  ] 142,148,144 

142  = 

' J n 1’ ' ) 14  3 

14  8 8 \/  = 1 ( J A K ) 

nil  I -J  14  3 
1*4  nArYC'iKKl 
I 4 3 J4  I'  1 .((U  t 
• 18  8=^8441 

Ih  ( •'18  8.  (jlr.  . 3 J (Jj  i J lUO 
8 (■  i'  J 8 'I 


Subroutine  GAUSS 

The  subroutine  GAUSS  is  called  by  the  subroutine  DRYDEN.  Each  time 
GAUSS  is  called  it  generates  a new  random  number.  Random  numbers  generated 
by  GAUSS  are  white  noise  having  a Gaussian  distribution. 
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Nomenclature 


Subroutine  GAUSS  (IX,S,AM,V) 

IX  Arbitrary  starting  number,  IX  = 65549 

S Adjusting  number,  setting  S = 1.0  gives  standard 

deviation  of  unity 

AM  Adjusting  number,  setting  AM  = 0.0  gives  zero  mean 


[INPUT] 

[OUTPUT] 

[INPUT] 


[INPUT] 


V Random  number  output 


[OUTPUT] 


Listing  of  Subroutine  GAUSS 


Id  I 1 Mfc  cAU6,S(  i A , i) , A '1 , V ) 
K c.  A 'j  ♦ I 

A = 0 . 0 

Dll  SO  i = l , Iz 
li(=lA*h3D30 
1 1-  ( i f ) b , b , o 
1 ii=li+21't74'iio4/  + i 

n y = I V 

[ x=  I y 
b 0 4 = A y 

V i-  c A - b . 0 ) ♦ S + A m 
K hj  I'  I '>1 

H) 


Subroutine  INIT 

The  subroutine  INIT  is  called  once  by  the  subroutine  STB  to  initiate 
storage  of  the  data.  Subroutine  INIT  stores  the  data  according  to  the  grid 
system  arrangement  shown  in  Figure  3C-1.  The  data  are  stored  as  WDEL  (I,J,K) 
where  each  row  designated  by  J corresponds  to  the  wind  difference  at  that 
elevation  as  explained  in  Section  3.1.  Each  column  designated  by  I repre- 
sents a prescribed  value  of  y ranging  from  -333.34  to  216.67.  The  index  K=1 
represents  longitudinal  wind  speed  and  K=2  represents  lateral  wind  speed. 

The  stored  data  is  transferred  to  the  subroutine  WINDF  through  a common 
statement. 
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JT 


J 


WOEl. 


2 3 


Figure  3C-1 


Grid  System:  K=1  Corresponds  to  AW^  = W^(2 
and  K=2  Corresponds  to  = Wy(z)  - fly(z  = 0 


Listing  of  Subroutine  INIT 


,3i'iw  j.i  J L Mr,  i M 1 1 


I V f I ,<■  / /;  ( .j  I i.i  c.  H J , 1 1 , / ,1  , 1 i)  . i U ^ IH 1 1 ^ <t  ) » t 
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/‘lo  r 1 1 K >1  A 1 ( H *■  1 0 . 4 ) 

4 0 0 

fc,r.U  lA  DAIa  hO 
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r'  c.  c u i iM  u « 1 « rt  1)  t . u 
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OUc-F 
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(II  3 ( = 1,34 

3 -'cAOC'i,  100)  ( .MUt.L  l 1 , J , ^ J , jsl  , 1 1 ) 

^ ♦**♦♦♦*♦  * r r.l-.O  10  uAIa  i.)«vX/U^  AoU  U^'l/uZ 

1 A X = 0 . ] 3 
4 1 : = . UO  1 
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Ji  2 K=l,2 
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" M ( 1 , 1 1 , A ) = I rtuc.Ll  1 , 1 1 , M - I*  u’M  ( i , 1 0 , M ) / 0<i 
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|"J  2 l = z,10 
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^ CiiiliOc- 
n A I ) = ” 4 u 0 , 


•*I,K  A = 1 . 0 
1 A i’  I A *1 


. 0.15)  - W Jz) 
.15) 


y **  ( J A , 1 I ,/  ) , 
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Subroutine  WINDF 


The  WINDF  subroutine  is  called  by  the  subroutine  STB  for  values  of  u and 
z.  Interpolation  is  performed  by  an  area-weighting  method  for  the  velocity 
at  the  elevation  and  for  the  stability  condition,  p,  from  the  tabulated  data 
placed  in  storage  by  subroutine  INIT.  The  area-weighting  interpolation 
method  is  illustrated  in  Figure  3C-2. 


Figure  3C-2  Area-Weighting  Technique 


The  area-weighting  scheme  calculates  the  velocity  and  velocity  gradient 
at  the  specified  point  by  using  the  four  nearest  neighboring  grid  point 
values.  The  velocity  at  point  P is  given  by: 

'^p  ^ ^ ^ ^3*^3  ^ ^4*^4^ 

where 

A = + A^  + A^  + A^ 

The  same  interpolation  method  is  used  for  the  velocity  gradients. 
Nomenclature 

Subroutine  WINDF(AMU,Z,WX,WY,DWX,DWY,USTR,F,ROSTR) 

AML)  Atmospheric  stability  parameter,  p 
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[INPUT] 


z 

Height  (m) 

Wind  speed  in  the  x direction  (m  s"^) 

Wind  speed  in  the  y direction  (m  s’^) 

Wind  speed  gradient  in  z direction  (s~^) 

Wind  speed  gradient  in  z direction  (s”^) 

Friction  velocity,  (iti  s"^) 

Coriolis  parameter  (s'^) 

[INPUT] 

MX 

[OUTPUT] 

WY 

[OUTPUT] 

DWX 

DWY 

USTR 

[INPUT] 

F 

[INPUT] 

ROSTR 

Rossby  number,  Ro 

[INPUT] 

Listing  of  Subroutine  MINDF 

, r^  I ■ 1 1 1 o * I ' > 1 .i , , A ^ I , i;  I'-  A , ij  A > , 1 1 .s  1 k , r , r(  ij S 1 K , I Z ) 
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S'  *1  k C 1 i ),  I ■<  t J ■*)»''/»  ui , Lir  /^ , u«  !:.»•  , U Kt'.h’ , iiDt  K 

L I '•  ’!  I J / i>  I / 4 Ij 

I ) 'i  1 1 = •'  J - ■<  A U 
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.J=  I 
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.)=.!■►  I 

lt-(J»'jtS.  .Iz)  i>'  f 

If  {/.  . . Cr' ) tjii  III  b 
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1 c.  'I  = J - 1 
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I'  aIjI-  '1  = 1 . - A Lr  A 

k C,  1’  A = ( ^,  - J . 0 0 1 ) / U . H 1 4 y - I 1 6 » - 1 J 

■ I ArJ  r M = 1 . - Hr  1 A 
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V.  = I'UJ'l.  i tVi  , 1 )+hc  A^' 

S ^ :)r.b  ( I '*  J V , 1 I.  , i ) ■*•  Aii“  A I A ♦ i»L  C.L  t 1 .'liJ  , i , I J t-  M i,h  M»i.r  I H ♦ 

S • ! -.  u ( 1 1J  , 1 ^ , 1 J 
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V-  X = •'■'  X ♦ J S r r( 
y = '.  y ♦ 1 ! s I w 

1 ■ ‘ * = E N E ♦ P 
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u 1 II  1 
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I y>  y = 0 . 

A A = J i 1 *<  ♦ ( A bUCi  ( r<(.io  J;h  ♦ Z.+  i.x+l  l.Zb»AiAU*/.  )/b.<* 

. ' i ..I  I 

A E.  K [ 1 r,  (.  6 , 7 1 0 ; 

/lb  H JK  A A r c 3X  , ' Aou  lb  Ovi  > tr  rr,bl0.y‘) 
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I Z = Z,  + 'K‘iI  K/E 


APPENDIX  4A 


TABULATED  FRONTAL  DATA 

Tabulated  data  for  two  cold  air  synoptic  fronts  are  given  in  this 
appendix.  The  tables  have,  due  to  their  length,  been  split  into  six  parts 
covering  columns  1 through  20,  21  through  40,  etc.  The  tabulated  values  of 
W are  values  relative  to  the  storm  motion.  The  frontal  speed  is  given 
at  the  top  of  the  table  for  converting  to  its  value  relative  to  the 
ground. 

The  case  numbers  for  the  two  storm  fronts  studied  are  listed  at  the  top 
of  the  table.  Also  listed  at  the  top  of  the  table  are  the  horizontal  length 
scales  for  the  given  wind  record.  The  horizontal  extent  of  each  data  set 
varies  because  of  the  data  reduction  procedure.  Hence,  the  length  of  field, 
L,  in  kilometers  and  the  horizontal  grid  spacing.  Ax,  are  specified  on  each 
table.  The  vertical  extent  of  each  field  is  taken  as  500  m with  50  m 
vertical  grid  spacing. 
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Table  4A-l(a)  Continued 
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Table  4A-l(b)  Lateral  Wind  Speed,  W , for  Cold  Front;  Case  1,  0307, 


• no  0 003300. 


— -^oooooooo. 


OOr^tOOOO  D 30. >• 
30-«00000000 


oor~o.-.oo»«. 


0>'«0J3<N(N00  O'hO 


oo'«>/)omo  o'^^o 


3 <3  000  3 0 0030 


OOOOi'SfXO  OOO  O 
OOOOOOOOOOO 


oooooo«-oooo 

OOOOOOOO  3-0  0 


0«/^'^«a‘'n.>*^0000 


DOT  .>•  0303 


-oooroooo 


> <3  O D 3«  .X  O O O 3 O 


3 3 3>0  n3  0r»33 


OOO'^-'^OOOOOO 


0'*^/>  3000  V o. 


000000«^D  0C«O 


0000000>DO<0 


3>3«.Dy)00^0000 


O 0/>0>x».%  3033 


•/*»vO  0>T.NOOOO 


^ /«  .*>  .N  -*  O O 3 0 


0300000 


9>0«</>00000000 


■ r-'Toooooooo 


O O.'N  03300  3 00 


^ roooooooo'^ 

03303333333 


i3>  O'*  C O^f 


■o 

O) 

3 

C 


o 

CJ 


a; 


fO 


»00300000^3 

i/1^00000>«0'*>0 

^O*r»0000^0-^0 

A^O'*tOOOr>  OOO 

."N  •■o  »»  .x  ^ o c» 

.XOi  5 3 

OOC.'Tv/> 

O 30J»0«N^0«-^0 

300300-*»-**-'00 

ooo  30»-‘«-»«-'« 

33.'^-0'^or' 

^300^.'N»-*-*»-*03 

t 03300<ro^ 
0000»-*tN  ><.N«30 

300^. •'•.■^*-’-*■-'3^ 

030^  ♦ *j.<»  0 3'^^ 
3033-«  — '^*-  — 


003'^'*'  '■■*••"• 

03030»-<*^*^003 


003uj'n»r>if>«f^r><»« 

ooc^oo^-^^^oc.  o 

OOOO'^O.OO'Of^*-* 

00000»^«^^330 

3000^0J>»/''^0 

03  C*  33^03003 

OO33»W00*^^3 
0303  3 300000 

0303no'»»r^'*3 

03330333030 


v'»^-<33>3i/>r<43on 

mOJ'B  CV3r>«'»3l/> 
y»3^«^^«r»V'*>/NO 

T^  Of'  D 03/)C'«  BO-/> 

3r»  r««rw^(N  BO^ 
w«r  »» 

•tncx  *»  n.>i«»»3>33 


3 3 3- 
'T  ♦»  V « 


I 3 3 ^ 3 -n 

•r  • 3*  ^ p«i  3 


Or-P'OOOOO  300 
“N  ."N  -N  •H  ^ 3 

O'^no  333r>-«0  3 

0 ^ rw  .N  O O O 

.>4  3333  O Oi/)>-<33 

333-^-<^->.-*033 

1 

^330^«-«.'*fHv'>4  03 

3033«-*«-'*^-»000 

I 

.'-  330f-r-  »>.->*0  3 

330  330  30000 

1 

jkOO.^-ni^i  ♦f^-^oo 

33003330030 

I 

-•33»^r»r»»<*>30 

WOOOOSSOOOO 

I 

»*/>  3 — 303 

.^SOO-^^x^OOOO 
I I 

000.>«  O 03<n333 


WtAOCNOOOO^OO 

•^'^oor«-'^''^('4^oo 

« I 

oootnooo<rr‘00 

I I 

0in0r-03OaDO<'<«  3 

330r<««'»«»‘NO*30 

• I 

>wi3D330r^3«»3 
3r4.X«  ^ f V m NO  3 

rny>303303000 

NN»•^»T^f•''N3J 


«rnCN<«0'^r*-'*t30i3 

0333-/>i3«r'*'<r'N3 


<333333«^'r{NO 


i^r*>4-43»nr“J^333 

«33  333rf>V«X^-N  3 

^ ^ o »N  o o 3 

<33  O Oj>3«  r«-NO 

n iN  •-•  3 3 .N  3 3*  3 3 

»3333j13<r-33.‘N3 

.■N^^03n  0<^»‘03 

>••••••••••• 

>033  033««*%>*0 

(V-<M33.N3<n3«33 
» 03  0333‘<r  »'*>  NO 

•-i«>«0  33^  O 33 

>333333^«'*^<^3 

•^^033^3. >4  ^03 
'»••••*•■•••• 

> 333333'r«''*><^0 

■•430r'*M3fNt>  33 
>333333V<'<*>fNO 

OO  »3  ''^.*4  3<N^03 

> 333333»»'*'-"**0 

3»0^^03'N»»33 

>333333'»<f^*'*0 

r‘3^-»033.N303 
>303303  T»^.'>«3 


i O O » 


— or  O-^OO  r^«N< 


228 


B 

CO 


cr» 

II 


o 

CO 

o 


a; 

•/) 


c 

o 


o 

o 


o 


-o 

Q) 

<V 

Q. 

</) 


(J 


QJ 


I 

< 


<0 


ro 

tr> 

It 

X 

< 


£ 


000300  30000 

t t t I I 

0000«^*^'^0'/'  NO 

00300000000 

I I I I I • 

oooo**<«  nx*^  rNO 

OOOOOODOOOO 

• I I I i I 

0000^*^<N'’'»N.-<0 

ooooooooooo 

I t I I I • 

0000000003D 

• I • < • I 

3 0 03^0«-^303 

ooo  oooooooo 

I « I 

OOOOOOO^CNOO 

000003000^0 

OOOOOOOCNOO  3 
OOOOOOOOOOO 

ooooc»o.-«.-»{N.-«o 

OOOOOOOOOOO 


o o 

O O ' 


r000>^0«-<00 
> OOOOOOO  3 


oooo.«- 

o o o o o < 


I iN  M o 

O O o o o 


0000»^'HCN'^.'N.-»0 

OOOOOOOOOOO 

I i I I I I 

o oc>  "TfnrNO 

OO  300000000 

• I I I I I 

00  0 0-*fN'»>/*^.NC» 

> O O O O 3 0 

I I I I I 

ooO^.NiN(N<rmwo 

OOOOOOOOOOO 

I I I « i I I 

ooo  «<N«tNin  »>«^0 

OOOOOOOOOOO 

I • I I i i I 

ooo— •-^•^•■«<NrN^o 

ooo  ^>0000030 

I I I I • • I 

00<^-4-««H^(N'4>«0 

00000003000 

I I • I I I I * 

30'^0300-*‘-»00 

00333303000 

I I I 

3 3-»0O03  3 300 
3 3 3 O O :>  J 3 3 3 3 


o^.-nn^noo  3030 
30000000000 


or-r-j^oino  — ^«3 
-•oooooooooo 


»«0Z'‘^000'»t(N  3 
■ -•-<000  30  0 03 


or->inooooorN<^o 

-•000003  3000 


3*J»T<r'n»%-«0303 

O *... 

O OOOOOOOOOOO 

OO'iJ'  i^'T-NiNOOOO 
a» 

/»  -•03  3 0000303 

0300(N«-H0033 

X, 

-^-«-.00033000 

oxr^mooooooo 

li>  -»3300300000 


P*;r333033-^  •>  : 
003030000  3: 


0 3 3 3 3 


3 O 
O 3 


3 3 3 

0 3 3 3 


•%3-^33333333 

00000303333 

I 

— — 30 

3300000  3033 

til  I It 

^ ^ ^ rN  M >3  o 

303300333  30 


'•I  »>  O 

•3  30! 


► 3 3 

■ OOO  3 


3 0 0 
0 3 0 


O— 

03000300 


O O O O O O O 


rN  -•  O 
0 0 0 3 


330000  30000 

• I I « I • I I I I 

«0 

000003033  33 

I I • I I I I I i I 

,s  .N  « .N  'J-  'f  -•  3 

0030303  0 033 

I I I I I I t t I I 

03000  3 33030 


0 0 3 0 0 


fN  (N 
0 3 0 


3 O 
I I 

-•  O 


(N  »N  -H 

3 0 0 


O O O O 


.N  — 3 


3 — 03 
3 O 3 O 


»N  N fN  ; 

O 3 3 V 


0 3 0 0 

I I i 

I — 3 
I 3 O 


3 0 0 3 
3 0 0 3 


3 0 3 3 
3 3 0 3 


• .«  N .N  .N  — . 


/■O  .-N  (N  — — — 


■3330330333  3 


3 0 0 0 


i^...aM-«(N  N30 
» 3 0033003 


0 0 3 0 


I i I I I 


{N  — O 
O 3 3 3 


, ^ ^ rN  3 — . 


33000003303 


<N  CN  — 

OOO 


O O O O 


<r  — 3 
3 0 0 0 


« O O O 3 


I O O r>  O 


V n (N  3 
0 0 0 3 


O O O fN  O • 


o — — 

ooo 


o o o — o — > 


o — 

0 0 3 


o o o o 


-•  o — 


m » JN  o 
o o o o 


o o o o 


rs  m — o 

0 0 0 3 


3 o 3 3 — — ; 


03303000333 


3 3 3 O O • 


— 30. 

3 O 3 ( 


> O 3 

> 3 3 


0—33 
0 0 0 3 


— 3 T*  O • 


229 


oo«tr>io*o30o<.^o 

OOOOOOOOOOO 


00000000900 


.•N0».''*>»-««.-«OC,O09 


>«««^oooooooo 


-«1{N.>trM^««0000  9 


oo-AOO/><r«.>«MO 

••■••••»••• 

r>«<XM.^MOOOOOO 


>««-«<^90000 


>oo/»ofl  nxi  N'^o 

N>*M»«-«900C»00 


j>An9.'>«r>>TOOO 
■n  .^09000 


999'«9d9'V-«99 

f'*.'#,>l^»-000c»00 


00<^Jir“i^^  N9  99 
IN.-'t^xOOOOOOOO 


M.-«£<^r>rv</t<«>ooo 

♦ '^>*.N.-»«00900 


noC'OO 

’X-XOOOOOOOOO 


O 9 o<^s  oo'nooo 

«.>*Mt-40900090 


OOOOOOOOOOO 
• OOOOOOOOOOO 


'VOOOOOOOOO 


'9  0/>0000r«900 


•^'^ooooooooo 


•^ooooooooo 


<>(««•••••••>• 

•-•OOOOOOOOOOO 


■ •-•I^OOOOOOOO 
• -hOOOO  OOOOO 


-»0*»/^-^.^0(NOOO 

rM^xOOOOOOOOO 


c>«oi<Nr>i.‘N(>4f-(Oooo 

•-•OOOOOOOO  90  o 


>^(N«-40009909 


• -4-4.-40000000 


n^n.N  900 



•-•OOOOOOOOOOO 


O 0>«-40000000 


oo«/>ni^'n.«.-*ooo 

ooooooooo  oo 


'*•^009 
90000  9900000 


'4  0000000 


9JBOOOOOO 


/>y>0'r>r/»'^*-*900 

S*>4t44»*>*< 

00900000  9000 


00-4000  9000  9 
90000  oooooo 


>i-40-4»>on.-«ooo 


00^/>^0^-<000 

OOOOOOOOOOOO 


'4-49t9090<9000 

OOOOOOOOOOO 


CN-40««'000<44-40U 


r-r-  oonr-vrx-^oo 

0**44*4»«**4 

OOOOOOOO  9000 


«>\(>t*^.^<—«400000 

ooooo oooooo 


O-f«-4(N^O«<*(^-<0C. 


a>  <Tt  r*  r-  -o  >o  >'n  -4  -<  o 

.n*********** 

C»  OOOOOOOOOOO 


4r«>«^r«-4>400ooo 

OOOOOOOOOOO 


iTt  9^0r^or-i«i-«00 


»9>aDr-r49ki/^o-4<-<0 

«*■••••••••• 

C OOOOOOOOOOO 


/1^r«•'4o*r49oooo 

00009000000 


04tOiUooiA<ri'«oo 

(44-4-4  _4«-«00  0 000 


O0’>«sa«0'0r^«4-40 
0-4-4000-40000  9 


0£  rt  *>.44.44  00009 
OOOOOOOOOOO 


V-49>  0-40  '^•4-4  o 

>(••••••••••• 

0-4-4-409-400000 


»f4or*»ri4  -4  0090 
>9990000999 


-44%0  r 44  44  -4  00  > 


— sikbr^oAr 


3 tt  0**>  0A**>>4« 


230 


Ci^’^CHC^CH.>4C-4r>t^o 


»-•>  o r*i  jy  ,■>*  > • 

o o»«ir^r^rN"^r>i^» 


f 

.y*  A t^  >4  />  O ©>»•*••-• 

I 

.-«<N  Off-  X'^k'UtN 

r»t»-»*>*<MCNr>ir^  >>«»<»* 

I 


<ytpr^>-'^ry 

O O M .X  .N  I 


I » r*  O ♦ 
I -n  .X  .X 


O O ^ iN  .N  .X  »x  iX 


J ^ .X  .X  .X  ^ ^ 

^r~^3'n7«^  O.-*»*-»d0 

o ^^x.x  x»* 

t 

»«r* 

^ ^ X »x  -X  •-• 


or-  oo>A»xj>^x»xj 

00»<X.Xl*^-'»'**'*>(X* 


Or>  4)00'«3'0'X><XX 
o o ^ .X  X *r  ‘*1  .X 


O'O^O  oxoooxr* 
oo^xx/^k^^^X"- 

O^Ot/^OXO'OXA 

oo^xx'n'^  »<^xi-« 

I 

0 o -<  X x .X  *»  ♦ -n  ,x  •-• 

1 

O O 

t»r'«-»xx.x^»*xx-< 


»r» 

■*\ 


^ w^  o OV  SOX>«f  D 

X«-*ClOO»»XX^^O 

I » 

0-^30WOOO»>'^0 

x.-.o-»^xx.x—  -<0 

I I 

0 />«^0000>7k-BO 

X DQ*-XXXX-«OD 

1 I 

-«3A03«^'^'^003 

mOO>mXXXXX«^9 

I 

00»^X  X XX  x.x*^o 
I 

XOO  OitXXX'^  ♦>/> 

oo*-«*^xxxxx»^s 

^or*'OOxxx.-«y>i» 

I 

0 /»  ♦ O x m X .X  X 3 

^D3^.XXXXX>^^ 

1 I 

OOXXOX^*><*»-^X  £ 

X ^ 3 ^ X X .X  X X -*  O 

I I 

o>nooox'*»'n«r*j^ 

X O X X X X X ^ 3 

I I 

ji^x«»»-*x»'n3x 

X <-•  3 «-•  X X X X X 3 

I I 

A3«3>«-|  >«>r>r3  oo 

^ -»  3 « X X X X X -«  3 

I • 

^00>'4X.n«303 
^ O 3 »•  X X X X X »■•  3 


•-•  3 O X X X X IX  •*.  - 

I I 

>xx/)in'^  o 3 

^ o o X X X i>«  X 'X  3 

« I 

a«o  ^ox  -rt  jy  -oo 
oo-«xxxxxx*«o 

I 

^<*>m3X»'n^^»<A«-t 

00-xXXXXXX»-i3 

I 

.A  X tT*  O X '•»  A •»  « </>  «x 
00«XXiXXXX-x3 
I 

30  A3X'^  A •^O^AA 
u ^ »•  X .X  X X X X « 3 


30A3XAJ>  ■AO 

3».4-.xxxXxx-x3 


o>«BX03>'xnooo 

♦ X-xO»x»^XXX.X-* 

I I I 

O Arx—  ^-x^'Aor-  ♦ 
«XmO  -•XXXX»x»^ 

I I I 

0«’fx3>O0X>^X 

♦ X-x33XXX-x^^ 

I • • 

oxo^oaAfxo 

»X^OO«|^XXX^.- 
I I I 

3 OA«3S«’Ar*<-33 

r x^Oxxmxxxxx 

I • I 

3onnA^o3x»x 

V X ^ 3 •-  tX  X A A X X 
» I I 

3‘>«rxO<AXAoey> 
r X *x  3 oX  X A A .A  .X  X 

* I I 

OO3«>O»0OOXX 

»AOOXXA»^A.X 

I I I 

3AA<>«3X00003 

« X 3 .-«  X A ■V*'  » ^ A .A 

I I 1 

3»40<-iO«x<nDOOOX 

»XO*xXAA^»«fA 

I I 

/y  BO^OOO  B O O B 
A ^ 3 .-4  X -A  A A A .A  X 

( I 

•-•A3«^3  OAr“X-x3 

A*«Of^XXA  AAiA.X 

I I 

0X303XAA*^  B 

X.-iO^XXAAX>*^ 
I I 

XX003  O^xAfx  CA 

X -x  o « X X A A X X "* 

I I 

X X 3 O 3 A lA  X O A A 

X »«  3 — X X X A X X -4 

I I 

XX030Ar-3AA0 
X *x  O *-•  X X X A X X *•« 

I I 

4xXOOO>XA  BAX  B 
X4x04xXXXXXX& 

I I 


-4*-«oa‘0— ♦r-X'^A 

A •,»•••••••• 

^ X-403XXXXXXO 

I I 

►x-xO»03XA^3A 

X ••••••••»•• 

I I 

•-•«x3X«X7>  3A3X3 
T X.-.33AX4XXX-43 

I I 


•x  O ^ 0 ' 


.-4  3 »•  o » 


O A » A X -4 


•4  3 94  B < 


O A 0 A X A 


231 


r 


■o 

a; 


o 

o 


CSJ 

I 

<: 


fC 


o'^rxw^oo^-^r-r-^ 

30»»^-^D»^0-nr»- 

ODO*^iN»»^»*r-nw 

a —•  (>4  eu  ■*  »»r'nr^ 

ooo««r.''<or<*«-*o  c 
•-•  rv  .•>*  V .-s  iv  ♦ » ♦ » r>< 

r'4>0001>00«<->0  3 

'*»aD003«^o«-<.'>«o/t 

I 

00000^0-<<N«^0 

(II 

0«7«000«^0««.>43 

I I I I I 

>300  OOOOOf^O  3 

v««  «>.>ior^«  # « >r 

I I I I I 

30r*y^»/»c^OBoor* 

♦ ♦•>  '#•-•o^N•^v,r  *» 

I I I I I 

^ ^ m oo  «ror«o  B*r 

» r n .>«••«  o •>!  fn  V .n  « 

I I I I I 

or- 

^ »kr»>«MOO.>i'*»  •» 

tor^r-  r*  0a)0«^ 

'n<*»r>»ooo.'N 

• III 

•Aor>B>t/iMO^OoO 


•^OBOO^  T r~ 

o 0-0 

I I I I I I I I I I I 

'*f>3‘l>00««B’r*0'«>-4 


r>  o 3 -• 


n^-v/iyi/>/)/>OB'0 

• • i I I I I I I I I 

^•'<r*-B«3«^0'^ooo 

/>  rivvin  n f>  •/i  o-BB 

I « I I I I I I I I I 

n'>*BflBr>r*B'noB 


n B f ♦ F » < 

• I I I I I 

B -*>  T ^ n ^ ' 


T B B B 

• I I I 

<F  ■*■  3 O 


I f I I I I I I I I I 

»»^4»>«n<noc  oooo 

i • • I I t I I I I • 

J»*».''i^Orf>f»*nOOO 

iB  iB  «!•  >r  * >N  .N  n -4  .-^  ^4 

I I i I I I I I 1 I I 

0«BOJsBO  Br>>300 

BtB«<-B.>IM<-i3  0>30 

* I I I I I I I 

0 3J«OOB030iB«F 

BB'B>|«-I00>’4«-<00 

1 I I I I 

0>B3OBOn0'^MC 

B ^ >N  <•«  O >>l  N .'<«  r-«  3 

• ill 

oo.«>o>3Br'0<Br->-t 

B « cn  3 ,>  .>4  v\  *4  •-< 

I I I 

oor-r-»-»»n30.N»*»B 

t I t 

'^oons  y>  o r-  J'OT' 

I I 

.XO-BO  B«^0  T Br-^ 


' 3 .>«  rx  n «r  > 


I <-•  o «n  «r  ( 


I — »»»  « >» 


•B  I*)  m ( 


'^r-OM3'n-*or-«i,>n 

'->XXC>^^.N*»(X<M(X 

till 

Dtb  3 

'^.X«-43*-*>-»N‘X  X(X  N 

I I I 

r»D  ^•BOr>.xr>-4'*i  u 


X 3 X • 


X X X -« 


f*i  >x  o M (X  n ix  >4  o 

I I 

0 X r-  r.  rs  rx  r>  cb-  3 •-«  tB 

«r  '*>«-io<xrx(xnrx«-«i3 

1 I I 

OBOOO>4W)x>«t/«lO 

».xoOfxmi*»<*^<x-43 

I I 

33B**t(33OB0B(3 

XX.-ax'WX^F^^nX^tO 
I I 

3'^BPv.B333»'l>0 

3 3*<X*T«r-Frx-*3 


fxoocxfxBocxcx-ncr 

'*)••••>•••••• 

^BiBB>BiB>B  OB  B CB 

— 

.X  o C X X »n  B 3 .X  -B  «r 

x**>»»***»., 
mBiBBBiB<BJ>BBBB 
I t I I I I I I I I I 

XOOXXXBJVX-^B 
<-*y>y»B»By>y»<B.B  B B B 


X O 3 X X .X  B > . 


•X3>«XXBO..«X-B 



OBBVBBBBiB  O OB 
•-•  I I I ) I I • I I I I 

X3^»-  XXBB-»X-»l 

B*«»*>***»»« 

OBBrBBBBBOOO 

•-•  I i • I I I I I I I I 

.X  O > .X  B B -4  X «) 

OBB^BBB'BBBBB 


X3  >S-«.-4XBBmX-B 
OB<B*^BBBBBBBB 

I I I I I I I I I I I 

O •?•  •«  -4  X tB  JU  -i  X iX 
OiBB'FBBBiBB  OB  O 


» C»  O'  ^ • 


OB«B%FBBBtBiBB*B  B 


I 3 >-«  ^ V r 


OB  BVBtBBBB  B O O 

>’4||I|||||||| 

-^3  7>.-i-4.-«-rr>  oa-tx 

X**«*aaaaaa. 

DBBtftBBBBBOBO 

*x  I I I I I I I I I I I 

-no  B3a-*aoV^O^M 
DBBtrBBB'BBOOO 


OAT  A X ^ 


232 


lA 

E 


rs. 


O 

o 


0) 

o 


o 


o 

o 


o 


T3 

O) 

q; 

Q. 


rO 

OJ 


CVJ 

I 

< 


O) 


AOBV>>4An'nO'»'3 

3 

r** 

I • I « • I t I 

0^40  n/nA^o>o«4>o 

<r>  •■•••.••••• 

I I » I I I I I 

o>-<‘>>o«pr'r~-  oo 

s •••••••••■• 


o (>•  B o o B r 


I O (7>  O 
I (N  O d 


s o r 
r>*  o -H  ■ 


r*  jv  » />  o ->  o 

3 o 


i •>•  />  : 

•N  O -H  ' 


« O 3 'Tt  ^ 
V o»  r«i  D o 


s 

r*^ 


II 

X 

< 


E 

<VJ 

ro 

II 


u^OO 

t ) I I I I t I 

•-•N/nr‘>0.>0»-«03 
D ^ r*  .-•  o 


."N  <n  » -n  .M  sj 


M.'N3  o<BBr*j<*^a«o 
P'«  I*  •-•  .'M  »»  m 


tt  ,t  ^ ^ 

•—•:>*  m <ri  rr%  , ">4  Z>  O 


JN  O <n  -•»/*»  ,s  .•^  O O 


.-•  -M  -•>/>,>•  » 
N O ' 


I «M  o O 
.•Mrs 


1(11111 

r3>*.'NBOOO 

^oo 

I »•  t I t I t 
^ n ^ j*  ^ 03  BO 

riO«-«r«*.N'»>m.->*f'*oo 


o J\  o r'  o 

p>«o«-*(Nr«*fM  ■^.'•r«oo 


• mo>Bs/)3r>-o 

*C4c*e>*f>tt>4rH(^o 


‘ ^ a c*  </>  or'VOr'O 
•0*^P^P^CNr«».N(NOO 


3 o ^ rx  »»%  . 


• ■Xrx->iCX‘N.NOO 


9 3 o ^ r 


-H  3 > o .-«  O 


i»*o  BOiO'^r^r*  O'n. 

'*'Pxo»xo«r>#rx.'N(N»^; 


ooM^^or> 

•■'('^-•OOPXfNrX.-NOiO 


^9>o  BO«nr*  .**  oo  3 

•.■^^O^fNCNriPiCN^^O 

( I I I I I « 

toa>i>«3«~>>r  It  t-  3 oo 

P^«3»N<N">4-^'^*Nci  3 


9^0  7«(>«z>  3>Mr*>^or^ 

r<«MMOor^r^‘>«(Nrso 


B'X>O.NrX»-<OP' 


<r“»r3»-«»B  BO/no 
i'Hooirxrxrxrx.'MOO 


► (N  « 3 r* 


i I I t » 


rx  --«  o •-«  • 


B B 3 •(  3 


tH  -H  o CM  .'x  .-^  3 3 


3 />  B 3 .M  3 ' 


I /)  > B 3 n 3 


CM»x3'MfMMCMCM,M3  3 



^x.np-r«i(/>>iB3-^o 

CM»^0»-»CMfMfMrs.‘M00 

I I I I I I I * 

CMBcn^CM  O 3>(^0*r  3 

.M  « O •><  CM  .->*  M rM  M 3 0 

( I I I I I I I 

CM3*fMMCMOO>3vrC» 

•M  rH  3 •-•  .M  .>*  .M  MOO 

((  I I I I I I 

rM3>OB<M  oojvo^ro 

CM  .-4  3 ^ fM  CM  cn  .M  .'M  3 3 

( I I I t I I 

<Mf-‘.M3nr‘3>3'»3 
CM«^3fM.MM  n<N'M33 


♦ CM^3-x*'4'X  MCMrx  3 

(1(1111 

3f*-.->»  3 f*3**n^3  O 

♦ fM-*3«'M'M'MCM?M3 

((((((( 

30r'CMy>‘3'»»*v*xo  o 

»^-<3-*-MfM-MfMCM3 


o <n  I'  3 • 


V 4M  «-•  O • 


I CM  M CM  iM  3 


'«’«-•  3 B 


> CM  .-•  3 ^ .***  .*>*  .M  CM  iM  O 

(((((((( 

>30  OOOVV>MOO 


^CM^OCMCMCM.” 

((1(1 

'*>3ii^^  BO*< 


.N  .-M  3 

( I ( 

.MOO 


« o »n  f 


O 3 O «B  O 


CM  •-•  3 CM  .M  .•M  3 *3 

( ( ( I ( ( ( I 


■M  -X  O N M ,M 


I •>  rM  3 3 


«/>  '*»M»-*0»xm.M.MCMM3 

( ( I ( I ( ( ( 

03«n'*>^3>rB.Moo 

i>  ••••••,.••• 

B '«>  .M  M O >«  M .-M  .M  .M  .M  o 

B3X>'X.«>0>/>/>.'M<3  O 
B ••••••«•••• 

♦ CM  M O .t  fM  .M  .M  .M  .M  3 


•'M-<3'M.M  M «>■«%  M33 
( ( I I ( I ( ( I 

B M 4M  B 3 O >4 

CN  33«M‘MCM<^  M.-VOO 

l(((l(((l 

O SOCMVOCM  BO-O^ 

<M0»-«CMr>*pX>ix»'M<MC»0 

((((ll(l( 

OBrM>«>y)^;«>*o  ^B<X 
fMO-xCMrM«*>m<NlMOO 


30«303/>/>lM3B 
•>  .M  .-•  3 .M  .N  M <M  .M  3 


.•MOV* 

.'M  ^ ; 


O />  ./>  .M  o o 

« ."M  CM  .M  CM  CM  3 


r m o o.MO.n 

(CM^O^^CMCNMCMCM  3 


O V m c*>  < 

CM  3 .-*  •>*  rM  "n  • 


o « r*  *«x 

CM  .'M  3 3 


(11(1 
oNy>vB'n«f/'*-xr~o 
^o«-*CMrM'n-n.-M  M30 


cn  3 • 
^ 3 . 


V3Vrf>*f—  -«■  3 
MnnnMN30 


3 > » r*.  on 


cnfM-xO^-ifMCMCMCMlMO 

( ( I I ( ( I ( 

>CMxv^»-*pxOB(MO«^n 

^fM»MO»-*CMJMCM.*M«MO 

1(1(11(1 

^•*>CMr*».MC-  ©"MO^i 

V .N  O ?M  .M  -M  CM  CM  3 

(1(1(111 

o SfMr-r*  •-•'“.M 
^M.«3M  M.M  MM-^O 


o n V n M ■ 


233 


03900000000 

0(Br«mo>^o*«'«'4^«o 

•-•oooooooooo 

•xv^mOOOOOOOO 

0>sr>«00«-«<-i'4^0 

'^oooooooooo 


•r*  ooo 

>0000000003 


/»»»000-^0  DOO 
ooooooo  oooo 

(MCViMOOOOOOOO 

OOOOOOOOOOo 

ooooo  ooooo 
OOOOOOOOOOO 

jU'^oo^  ny^'vt.'NOo 


-•  "M  .N  o o o o r*  o 


00000000>'r^ 


oooo^>o>>i>-no 


• O O O O J>  -U  J 


^ o o o o <■ 


cmoooo^o  Oinoio 
ooooooo oooo 

<NOOOO^r«««{NO 

OOOOOOOOOOO 

f*>i000-0^^''*^^'^0 


>'»0  oor^.'NM.'^^D 


/>  o<ooooooooo 
ooooooooooo 

or>r«(Nooooooo 

ooooooooooo 

r*0D-»*r<*t(N^*  oooo 

OOOOOOOOOOO 

J»»«  6/>««0000 
OOOOOOOOOOO 

j*  >j>  xaDr“f«>*n-n  oo 

ooooooooooo 

.>AOOO#>>0/>'««d 

^^*4^,^000000 

o<sr*«r-«'^o>®'*>o 

oo 

o«»r‘*^r«-Of'<o>/>o 
.*N  «N  .>4  «>l  <*4  «-•  O O 

o^r<«oovo^9>«o 
«rm.^.'n<^(>«p4«-«o  oo 

<»«'’^.n.'*<>4«-«-400o 

OOi>  07>iAO^J0'no 
^<»»<N.N'4'<'^0000 

ooi>r*-<*to«^jDr*-'No 

V.'«M'4'4*<00000 

OJ>0>J«l>>0('>0  NO 

ooooooooooo 

or*  O o OJ»r*  o«nriio 
^N  oooooooooo 

oo^o  ^r>o/»#rNo 
'OOOOOOOOOOO 

OOOOOOOOOOO 

oo'4'4'4i/\^mfn'«  3 
C OOOOOOOOOO 

Oi'*oo»>«i/^or**iN'«o 
OOOOOOOOO  Oo 

OOO  OO^'^  '*.■>•'43 
OOOOOOOOOOO 

<^000  O •>*r*-4  3 

003. 3 ooooooo 


OOOOOOOO  30  o 
•-•OOOOOOOOOoO 

0000000-3  000 

(>«••••••••••• 

O«00000000000 

ooooooooooo 
#>«  ooooooooooo 

ooooooooooo 



''OOOOOOOOOOO 

0*4000000000 



000000000030 

O'^OOOO  30000 

oooooooooooo 

0.‘N00O'«'«'*OO0 

r************ 

00000000-0000 

0'n'rt00'«''''000 

o»*»******»« 

oooooooooooo 

omr«oo'«'''«ooo 

>/>••••••••••• 

oooooooooooo 

'«vm0O''''''000 

«••••*•••••• 

oooooooooooo 

IN|^«*>00'«'''4000 

ooooooooc  oo  o 

■^yt«00'4'''400  o 
oooooooooooo 

«0,n00'00''000 

000003030030 


*4  3 > « *4  (t  n < 


o n r •>  N 


234 


lA 

E 


r««. 

II 


O 

o 


0^ 

O 


ooodooo«40^o 

9 

r>«  OOSO^OOOOOO 

I 

ooooooooooo 

r-  •■•••••••■• 

oooorN3ooooo 

-«  0000f-i030003 

till  I 

•-•t'«(>«OISr4«>40000 
r*  ••••••••••• 

03000000000 

• III  II 

o ••••••«•••• 

'I  00000030000 

I I I I I I I I 

n 

ooooooooooo 

I I I I I I I I I 

« ••••••••••• 

^ ooooooooooo 


rx  f>i  <M  .X  . 


• I I I I I 

^ rx  rx  o 


00*x00-400000 

OOOOOOOOOOO 

I 

•-t'x^v/nrx^ooooo 

OOOOOOOOOOO 

I I I I I I 

o^*-<oo»«ooooo 

OOOOOOOOOOO 

I I 

rxfx»nrx"x<N  oooo  o 

OOOOOOOOOOO 

'n»0/>^00i0-x0 

OOOOOOOOOOO 

I 

•N  n00«-«0  OO  0-4  o 

OOO  oooooooo 

I 

rxrxooo-4000-40 

OOOOOOOOOOO 


•ot-iOOO^OOO- 


-4  OOOOOOOOOOO 


o 

o 


o 


•o 

a> 

01 

Q. 

I/) 


<o 

U 


s. 

01 


CSJ 

I 

< 


.o 


6 


II 

X 

< 


E 

CSJ 

CO 

II 


-4  rx  (N  o • 


I I I 1 I 

I ^ ^ -4  O 


o o o o o o ; 


-4  o o o o > 


OOOO 


• M o 


-4  OOOOOOOOOOO 

till  I I I I « 

-^-4-4O<-i00-4  -4  00 

-4  OOOOOOOOOOO 

III  II 

000<-4-400000  3 

o ••••••••••• 

44  OOOOOOOOOOO 


rx>n-4-4-40ooooo 

OOOOOOOOOOO 

-rw>'»»--x'voooooo 

OOOOOOOOOOO 
t'X'41  'X<->«40  0-444-40 

00000000030 

44044.4-4000000 

OOOOOOOOOOO 

-444000000000 

ooooe^oooooo 

I I 

*x<4\4400000000 


10000044-40 


■oooooooo 


O .>OOOOOOJOO 

o^oo-ao '0  0330 


•00000030000 


oooooooc 


3 O O O O o < 


oooo 


I I I I I I 


0-4  >*  •>«  rx  X . 


♦ »».-X  -4000030 
30000000000 

«0</>«<4>-444-4000 
OO  3000000  OO 

xar»r4./>^'nfx-4-4-40 

OOOOOOOOOOO 

Win  noinw  4irx44-«o 
00000000300 

nW‘*^.X  44  -4  44  0 
00030000000 

•x.'x'N  »'n-x-4-40  OO 

000003  oocoo 

-4044m44-4-40ooo 

OOOOOOOOOOO 

-4CX-4(N000-44400 


I x%  fX  44  (N  -4  « 


0000  OOO  oooo 

1 I I I I I I I I I 
wonwww  r4%,x440 
OOOOOOOO  OOO 

rf>nr‘r“0^*».x'x-40 


> O O 3 0 0 0 0. 


-4“X.X  X*XX-4— 1-4-40 


p o o n w ><  ' 


' O O O vO 


-4  Xi*l44-4-4-4-4>4-40 

OOOOOOOOOOO 

I I I I I I I I I I 

-•rx.x  -4-4-400  0 0 o 

OOOOOOOOOOO 

i I I I I I 


• rxoo -4  000  0 0 


0-4-400000C  OO 

OOOOOOOOOOO 

I I 

0-4-4  00000000 
OOOOOOOOOOO 


0000-4000000 

OOOOOOOOOOO 


O 3 O O >4  44  • 


-4  3 > •!>  ■*  o n 


I rx  44  44  -4  o o o 


■ o o o o o o 


00000-400000 

OOOOOOOOOOO 

-444'*)4trx000333 
00303C  30000 

*«4>/t/t^-430000 

30000000000 

0-4rx'xx.M030C3 
33333  3 30303 


^ o »>  » •»•  p n 


?35 


■o 

0) 

3 


o 

o 


CSJ 

I 

< 


r^r>-r''«r>oo»oo  P 

^M^«»aooor>ir*«o 

till 

as»s»«r«oooox» 

•^MMMv^aoo'xota 

(II* 

aaoaA«a9>oao 
f**p*rM: 

»»  o r*  o I 


> P O >4  .>• 


p>a’«t  pr<i>  3 Bai>4: 

I I I I 

p « « ( 

r»  r*«  a ix  ^ ' 


.'x'tav  naooo'x  . 


I I I 

’OP-**rf>a'^OtN> 


^ a r**  ,x  ,x  p j%  a -4  > 


.X  o A a a . 


P •■  ♦ -X  ,X  O O r4  ,X  .X  .X 
I I I I 

PfXPPP-0«OOOr» 


j»  ♦ /•»  .X  . 


I p a P O r 


■nasopofxoaaa 


pr>^^p>4«>«or 


i/>po«taooM<Mor« 


a‘«toar*«/>o<^'4ar> 


p r«  ^ p a • 


Moijk'npooi'x. 


<x.x>4p»«r««xp««PO 
Bppr*>r>  ppip^r^tr* 

rx«x«»o  p-^r-oixo 

ppoDr^r^p  p«pp‘*t(x 

P'S  pr^p  oppp'nix 

•-••«oor“»ar«-«-*pp 



»p  pp4  ppp«r>'«’'*i(x 

•4*^a  o ppp'^cxpa 

PPPOOPOPP'nxt 

r^aa  op**>ap'«>ap 
PPBP  0 P 0/>p4r‘«f 

oo'n^f'tfxo  o ■nao 

ppr-  ppp  OPP«« 

oor-.x(x»4  ap««40 
pr>pp  op-op<r«« 

oaa  oppy^po'^a 
pppT-yvvt^^p^ 

oaoa«^aoaapo 

«“  9 «/n/*>»*»<*>(x.xi4»  » 

aaaPP'np<4>p»a 

>X(X 

a a poaaar«ap«r 

•x  X .4  ,x  a a a a •-■  x .4% 

aap  paa<-40Pa9> 

■xx^^aaaoox.x 

I 

aa  oaa^^^ap«> 

'X'x.^aaaaaO'-'M 

I I I 

aap^ra^.xi^^/noj** 

(x.x^aooooa^«>4 

I I I I 

oa</>pa<xrx^pp«H 

(xrx<-4oaooooo<>4 

lilt 

ooy>x4«M^poop 

(x(x«>i00oooooa 

I I I I 

a DPo^a«HPi>4no 
M.M^^aaoo^oo 

I I I I 

xf**/>M‘xa.-4'o«fa(x 

M^^^apap-«>4a 

I I I I I 

nn/»«r*»a-*7‘r«n«r 


-«-*»«  a 


.X.X^OOPOOOfXP 

m*********** 

•4Ppxar4PPi/>«mN 

(X  .X  »4  a IN  p >n  p p o 
:««••••••••••■ 

^ppppr*  ppp^fvi/n 

/*>m<xa*np»-ir*ar*(x 
•tpp  BPr^ppytpfvtm 

m .«>  ,x  o p a ix  a > o o 

a»»«*«<*«*** 

MPPPPI*'  PPPV  pivt 

n ««  .X  M r*  o >N  o p a a 
>••••••••••• 

a ppp  nr^r'PP^  *i.n 

4rai«i'4p.N(xor*>  pr* 

p*a**»a»»«»* 

OPPPPr^r-ppp'n^* 

4»«onfxaf«-r»p 
r>  ••••••••••• 

aBapppr-pp^vix 

e ^ -a  a N .N  a p e # 

Paaa******** 

OP  BPPr^r*  PiO'pm.N 

•>4«4aPtOr«> 

^••••••••••» 

OP  sppf*‘r>pp«mcx 

i«t'nrxoPOMa4rp«4 

Vavavavaa*** 

CPPPBr4r>p«/>p>«4||x 

fnmixai^Oa^pmf^^O 

0»pP9r4r>py>p«n(x 

■41  .N  N a p > -4  p .N  N a 

■Xa,.**.****. 

opppp***  ppa»4»x 

rXN«4  B.N  Bi^PIXMa 
aoosr**>»POAf*»X 

-4P  ^B•*•  OA  P'AXa- 


236 


APPENDIX  4B 


GRAPHICAL  ILLUSTRATION  OF  FRONTAL  WIND  SPEED  PROFILES 

In  this  appendix,  the  wind  speed  profiles  which  an  aircraft  would 
encounter  approaching  along  a 3°  glide  slope  through  a cold  synoptic  front 
are  illustrated.  Four  flight  paths  having  3°  glide  slopes  and  spaced  at 
equal  increments  along  the  horizontal  extent  of  the  data  set  are  shown  in 
Figure  4B-1.  Figures  4B-2  and  4B-3  illustrate  the  wind  speed  profiles  along 
the  four  flight  paths.  One  observes  that  the  winds  can  be  significantly 
different  for  aircraft  approaching  at  different  times  through  the  same  front. 
Due  to  the  limited  spatial  extent  of  the  data  set  the  latter  flight  paths,  3 
and  4,  do  not  extend  to  the  full  500  m level  and  consequently  are  only 
plotted  to  the  height  at  which  they  pass  out  of  the  data  set  at  the  right- 
hand  margin. 

The  streamlines  on  which  the  flight  paths  are  drawn  are  based  on  wind 
speeds  relative  to  the  motion  of  the  front  which  is  9.3  m s ^ and  7.1  m s ^ 
for  Figures  4B-2  and  4B-3,  respectively.  The  vertical  dashed  line  on  the 
longitudinal  velocity  profiles  illustrate  the  mean  motion  of  the  front, 
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Case 


igure  4B-1  Flight  Paths  through  Streamline  Pattern  Moving  with  the  Front 
Corresponding  to  Wind  Speed  Profiles  in  Figures  4B-2  and  4B-3 


z i 


-15  -10  -5  0 5 10  15 

(ms  ^ ) 


Figure  4B-2  Wind  Speed 
Cold  Front 


Path  4 


ofiles  Along  Four  Different  Flight  Paths  Through 
lase  1 
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Figure  4B-3  Wind  Speed  Profiles  Along  Four  Different  Flight  Paths  Through 
Cold  Front,  Case  2 
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